Seasonal melting and the formation of 
sedimentary rocks on Mars, with predictions 
for the Gale Crater mound 



Edwin S. Kite^, Itay Halevy^, Melinda A. Kahre*^, 
Michael J. Wolff ^, and Michael Manga '^'^ 



^Division of Geological and Planetary Sciences, California Institute of Technology, 

Pasadena, California 91125, USA 

^ Center for Planetary Sciences, Weizmann Institute of Science, P. O. Box 26, 

Rehovot 76100, Israel 

'^NASA Ames Research Center, Mountain View, California 94-035, USA 

Space Science Institute, 4^50 Walnut Street, Suite 205, Boulder, Colorado, USA 

^Department of Earth and Planetary Science, University of California Berkeley, 

Berkeley, California 94720, USA 

^ Center for Integrative Planetary Science, University of California Berkeley, 
Berkeley, California 94720, USA 



Number of pages: [60 



Number of tables: 1 
Number of figures: 19 



Preprint suhmitted to Icarus 



5 March 2013 



4 Proposed Running Head: 

5 Seasonal melting and sedimentary rocks on Mars 



6 Pleaise send Editorial Correspondence to: 

7 

8 Edwin S. Kite 

9 Caltech, MC 150-21 

10 Geological and Planetary Sciences 

11 1200 E California Boulevard 

12 Pasadena, CA 91125, USA. 

13 

14 Email: edwin.kite@gmail.com 

15 Phone: (510) 717-5205 

16 



2 



ABSTRACT 



18 A model for the formation and distribution of sedimentary rocks on Mars 

19 is proposed. The rate-hmiting step is supply of liquid water from seasonal 

20 melting of snow or ice. The model is run for a O(IO^) mbar pure CO2 atmo- 

21 sphere, dusty snow, and solar luminosity reduced by 23%. For these conditions 

22 snow only melts near the equator, and only when obliquity ^40°, eccentricity 

23 >0.12, and perihelion occurs near equinox. These requirements for melting are 

24 satisfied by 0.01-20% of the probabihty distribution of Mars' past spin-orbit 

25 parameters. Total melt production is sufficient to account for aqueous alter- 

26 ation of the sedimentary rocks. The pattern of seasonal snowmelt is integrated 

27 over all spin-orbit parameters and compared to the observed distribution of 

28 sedimentary rocks. The global distribution of snowmelt has maxima in Valles 

29 Marineris, Meridiani Planum and Gale Crater. These correspond to maxima 

30 in the sedimentary-rock distribution. Higher pressures and especially higher 

31 temperatures lead to melting over a broader range of spin-orbit parameters. 

32 The pattern of sedimentary rocks on Mars is most consistent with a Mars pa- 

33 leoclimate that only rarely produced enough meltwater to precipitate aqueous 

34 cements and indurate sediment. The results suggest intermittency of snowmelt 

35 and long globally-dry intervals, unfavorable for past life on Mars. This model 

36 makes testable predictions for the Mars Science Laboratory rover at Gale 

37 Crater. Gale Crater is predicted to be a hemispheric maximum for snowmelt 

38 on Mars. 



39 Keywords: MARS, CLIMATE; MARS, SURFACE; MARS, ATMOSPHERE; 

40 GEOLOGICAL PROCESSES; MARS 
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Climate models struggle to maintain annual mean temperatures T > 273K 
on Early Mars (Haberle, 1998) (Forget et al., submitted manuscript, 2012). 
Seasonal melting can occur for annual maximum temperatures Tmax ^ 273K, 
which is much easier to achieve. Therefore, seasonal melting of snow and ice 
is a candidate water source for surface runoff and aqueous mineralization on 
Mars. Surface temperatures ~300K occur at low latitudes on today's Mars. 
However, seasonal melting of surface-covering, fiat-lying snowpack does not 
occur because of (1) evaporative cooling and (2) cold-trapping of snow and 
ice near the poles or at depth. Reduced solar luminosity for Early Mars makes 



melting more difficult (Squyres and Kasting, 1994). 



Toon et al. ( 1980 ) modelled the control of Milankovitch cycles on Mars ice tem- 



peratures. Melting is favored when snow is darkened by dust, and when evap- 



orative cooling is reduced by increased pressure. Jakosky and Carr (1985) sug- 
gested that equatorial snowpacks would form at high obliquity. They pointed 
out that melt could contribute to the observed low-latitude erosion. IClowl 



(1987) modeled snowmelt due to the solid-state greenhouse effect. He tracked 
meltwater migration to the base of the snowpack. Several authors have mod- 
eled melting on steep slopes as a candidate water source for young midlatitude 



gullies (Costard et al., 2002 Christensen, 2003). Hecht (2002) considered the 



energy balance for water at the melting point in gully alcoves. Williams et al. 



(2008, 2009) modeled melting of relatively clean snow overlain by a thin, dark 
lag deposit. They found melt rates ~1 kg/m^/hr on steep slopes, but argue 
that this is sufficient to form gullies through either fluvial or debris-flow in- 
cision. Morgan et al. (2010) used a ID atmospheric model to examine water 



ice melting and CO2 frost accumulation. No author has studied Early Mars 



cold-traps, although Schorghofer's model (Schorghofer and Aharonson, 2005 



Schorghofer 2007b 



2010) has been used to track cold-traps for subsurface ice 



over the last 5 Ma (Schorghofer, 2007a). 



The flrst purpose of this paper is to extend the global snowmelt models by 
integrating a new thermal model over all spin-orbit parameters, while account- 
ing for cold-traps. Chaotic diffusion in the solar system makes it almost certain 
that Mars' obliquity (0) has ranged twenty times more widely than Earth's 
obliquity over billion-year periods, and that Mars' eccentricity has had a long- 



75 

76 


il^erm variance 


twice that of the Earth ( 


Touma and Wisdom, 


1993 


Laskar and 


Robutel 




1993 




Laskar et al. 




2004 




Laskar 




2008 


). These wide swings cause 



large variations in insolation and propensity to melt (Figure [T]). 

78 The second purpose of this paper is to understand the water source for sed- 

79 imentary rock formation on Mars (Malin and Edgett, 2000). We focus on 



the hypothesis that supply of water from seasonal melting was the limiting 
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time during interval of sedimentary rocl< deposition 



Fig. 1. Motivation for this paper. Mars underwent tens to thousands of spin-or- 
bit oscillations during the interval of sedimentary-rock deposition. Three are shown 
schematically in this sketch. The geologic record of metastable surface liquid water 
is a wet-pass filter of Mars climate history. Mars orbital parameters vary over a 
wide range, resulting in a correspondingly wide range in tendency to melt. There- 
fore, evidence that the sedimentary rocks formed in a small fraction of Mars' history 



(Lewis et al. , 2010) suggests that negligible melting occurred under mean orbital 



forcing. If Mars sedimentary rock record only records orbital conditions that per- 
mitted surface liquid water, modeling average orbital conditions is neither sufficient 
nor appropriate. Instead, it is necessary to calculate snowmelt for the full range of 
orbital elements that Mars likely sampled over the time interval of sedimentary rock 
deposition. These predictions can then be compared to observations. Because of the 



evidence for orbital pacing of sedimentary rock accumulation ( Lewis et al. , 2008 ) , 



transient warming events are not shown, but may have been critical for generating 
geomorphically effective runoff - see §8.4. 



step in the formation of sedimentary rocks on Early Mars (e.g., Wendt et al. 



2011). Existing evidence for snowmelt-limited sedimentary rock formation is 



discussed in 
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If surface liquid water availability was necessary for sedimentary rock forma- 
tion, then the spatial distributions of liquid water availability and sedimen- 
tary rock detections should correspond. ^ analyzes the global sedimentary 
rock distribution. In the only previous global model of sedimentary rock for- 



mation on Mars, Andrews-Hanna et al. (2007) tracked groundwater flow in a 



global aquifer that is recharged by a broad low-latitude belt of precipitation. 
Groundwater upwelling is focussed in low-lying areas, generally consistent with 



the observed distribution of sedimentary rocks (Andrews-Hanna et al. , 2010 



Andrews-Hanna and Lewis 2011). This model assumes T >273K, in order to 



avoid the development of an impermeable cryosphere. Especially in light of 
the Faint Young Sun predicted by standard solar models, temperatures this 
high may be unsustainable for the long periods of time required to form the 
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sedimentary rocks (Haberle, 1998 Tian et al. , 2010) (Forget et al., submitted 
manuscript, 2012). The model described and analyzed in this paper 
assumes that liquid water is supplied from locally-derived snowmelt, rather 
than a deep global aquifer. Groundwater flow is confined to shallow, local 
aquifers perched above the cryosphere. Annually- and planet-averaged tem- 
peratures remain similar to today's, which reduces the required change in cli- 
mate forcing from the present state. If Mars climate once sustained T >273K, 
then it must have passed through climate conditions amenable to snowmelt 



en route to the modern desert (McKay and Davis, 1991). The converse is not 
true. 



106 
107 
108 

109 
110 
111 
112 
113 

114 
115 
116 
117 
118 
119 
120 
121 
122 
123 



Model predictions for different paleoclimate parameters are compared to global 
data in ^ ^ makes testable predictions for Gale Crater, the target for the 
Mars Science Laboratory (MSL) mission (Milliken et al. , 2010). 



The discussion (^ includes comparison of the snowmelt model to alternative 



hypotheses such as global groundwater (Andrews-Hanna et al. , 2007, 2010 



Andrews- Hanna and Lewis 


2011 


) and cryogenic 


(Catling et al. 


2006 


Niles and Michalski 


2009 


)• 



are stated in ^ 

This scope of this paper is forward modeling of snowmelt production as a 
function of (unknown) Early Mars climate parameters. Beyond a qualitative 
discussion in ^ and ^ there is no attempt to physically model the pro- 
cesses running from snowmelt production to sedimentary rock formation. A 
computationally inexpensive ID model allows us to sweep over a large pa- 
rameter space. The trade-off is that ID models cannot track the effect of 
topographically-forced planetary waves on the atmospheric transport of water 



vapor, which controls snow precipitation ( Vincendon et al. 2010 Forget et al. 



2006 Colaprete et al. , 2005). Any ID snow location prescription is therefore 



an idealization. 



124 2 Snowmelt hypothesis 



125 Liquid water is required to explain sedimentary rock texture and bulk geo- 

126 chemistry along the Mars Exploration Rover Opportunity traverse, and there 

127 is strong evidence for extending this conclusion to other light-toned, sulfate- 



128 
129 


J^earing sedimentary rocks on Mars ( 


Bibring et al. 




2007 


McLennan and 


Grotzinger 


2008 


Roach et al. 


2010 


Murchie et al. 


2009a 


Weitz et al. 


2008). 



130 The hypothesis in this paper is that the water source for sedimentary rocks 

131 on Early Mars was seasonal melting, and that liquid water was infrequently 

132 available so that melt availability was the limiting factor in forming sedimen- 

133 tary rocks. "Sedimentary rocks" is used to mean units comprised of chemical 
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precipitates or siliciclastic material cemented by chemical precipitates, usually 
sulfates. These are recognized from orbit as light-toned layered sedimentary 



136 deposits (Malin et al. , 2010) that characteristically show diagnostic sulfate 



137 features in the near-infrared. This definition excludes layered phyllosilicates, 



138 which usually predate sulfates (Bibring et al. 2006; Ehlmann et al. 2011). 



2.1 What is the evidence that sediment lithification on Mars requires liquid 
water? 



141 
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Erosion to form cliffs and boulders (Malin and Edgett , 2000), ejection of meter- 
size boulders from small, fresh craters (Golombek et al. 2010), resistance 



to crushing by rover wheels, and microscopic texture ( [Okubo , 2007) show 



that most light-toned sedimentary deposits (hereafter "sedimentary rocks"' 
are indurated or lithified. Lithification involves compaction and cementation. 
Water is required to form aqueous cements, and for fluvial sediment trans- 
port. At the Opportunity landing site, evaporitic sandstones (60% chemical 
precipitates by weight on an anhydrous basis) record groundwater recharge 



et al. 


2005 


McLennan and Grotzinger, 2008). Aqueous minerals are present 


^n sedimentary rocks throughout Meridiani and the Valles Marineris (Bibring 


et al. 


2007) 


. Murchie et al. (2009a) argue for water-limited lithification of the 



Valles Marineris sedimentary rocks. Subsurface pressure-solution recrystalli- 
sation can occlude porosity and lithify weak evaporites at ~30 bars, without 



aqueous cementation (Warren, 2006). But in this case water is probably still 
needed to form evaporites at 60% by weight. Some layered sedimentary de- 
posits on Mars might not require liquid water to form, but these deposits are 
usually younger or at higher latitudes than the sulfate-bearing layered sedi- 



2010) 



mentary rocks (Hynek et al. 2003 Bridges et al. 2010 Fenton and Hayward 



Outcrop thermal inertia (Tl) is almost always low at 100m scale, so these rocks 
were either never strongly cemented or have been weakened after exhumation 



163 (Edwards et al. , 2009). 



164 2.2 When did sulfate-bearing sedimentary rocks form': 



165 Sulfate-bearing sedimentary rocks occur relatively late in the stratigraphic se- 



Mars ( 


Murchie et al. 


Mango 


d et al. 




2010 



ThoUot et al. , 2012). Formation of sedimentary rocks peaked in the Hespe- 



169 rian (Carr and Head, 2010), well after the peak of phyllosilicate formation on 



170 Mars (Ehlmann et al. 2011 Fassett and Head 2011 Salvatore et al. 2010 
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Bibring et al. 2006). Phyllosilicates are sometimes interbedded with older 



sulfate-bearing sedimentary rocks (Ehlmann et al. , 2011), but those phyllosil 



icates may be reworked (Barnhart and Nimmo, 2011). Sedimentary rocks also 



postdate almost all of the large-scale, regionally integrated highland valley net- 



,works of the Late Noachian/Early Hesperian (Carr and Head, 2010; Fassett 



and Head, 2011), and are spatially separated from these "classic" valley net- 



works (Hynek et al. , 2010). Therefore, the observed sedimentary rocks cannot 
be the terminal deposits of the classic valley networks. The climate that cre- 
ated the classic valley networks could have been different from the climate that 



formed the sedimentary rocks ( Andrews-Hanna and Lewis, 2011). Sedimen- 



tary rocks do contain some channels, often preserved in inverted relief (Edgett 



2005 


Burr et al. , 


2009 


2010 



2010). Finally, many sedimentary rocks postdate the 



large impacts of the Late Heavy Bombardment, and many have quasi-periodic 



bedding suggesting orbitally-paced deposition (Lewis et al. , 2008, 2010). These 
observations are inconsistent with brief bursts of rapid sedimentary rock for- 
mation during impact-induced greenhouse events. 



187 2.3 What existing data supports the snowmelt hypothesis? 
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Liquid water was in short supply even at the time of sedimentary rock for- 
mation at the Opportunity landing site. Mineralogy indicates low water/rock 
ratios during alteration, and that the cumulative duration of water-rock in- 



teraction was < 100 ka (Berger et al. 2009 Hurowitz and McLennan 2007 



Elwood Madden et al. 2009). Weathering at Meridiani Planum was either iso- 



chemical or at low water/rock ratio or both (Ming et al. , 2008), consistent with 



a rare trickle of snowmelt. Low specific grind energy of sandstones indicates 



weak aqueous cementation (Herkenhoff et al. , 2008). The present day extent 



of sedimentary rock outcrops is small, and the persistence of opal, jarosite 
and olivine in rocks (and olivine in soils) indicates minimal water-rock inter- 



^ction sin ce those minerals crystallized (Tosca and Knoll 2009 Olsen and 



Rimstidt 


2007 


Yen et al. 


2005 



2005). Away from the sedimentary rocks them- 



selves, aqueous mineralization was minor or absent elsewhere on the planet at 



the time when most sulfate-bearing sedimentary rocks formed (Murchie et al 



2009bl ISalvatore et all |20T0l [Mustard et "aL| [2009] [Fassett and Headj [20n 



Hausrath et al. , 2008). Globally, soils formed "with little aqueous alteration 



ijunder conditions similar to those of the current Martian climate" (Bandfield 



et al. , 2011), and elemental profiles indicate top-down mobilization of soluble 



elements (Amundson et al. 2008 Arvidson et al. 2010). 



207 Dividing the total thickness of sedimentary rock deposits by the thickness of 

208 quasi-periodic layers and then multiplying by the obliquity periods thought 

209 to pace accumulation suggests that the sedimentary rocks formed in 1-10 



210 Ma (Lewis et al. , 2010). This is a small fraction of Mars' history. Geomor- 
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phic evidence that the Mars surface environment has only marginally sup- 
ported surface liquid water since the Noachian includes mean erosion rates 
~1 atom/year (Golombek et al. , 2006), together with a sharp post-Noachian 



decline in valley network formation and crater infilling (Fassett and Head 



2008 Forsberg- Taylor et al. , 2004). These data argue for a short-lived and 



downward-infiltrating post-Noachian water supply, suggestive of transient liq- 
uid water that is generated only during brief melt events. 



218 2.4 How could brief pulses of snowmelt form kilometer-thick accumulations 

219 of sedimentary rock? 
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Antarctica's McMurdo Dry Valleys are a terrestrial analog for seasonal-melt- 



Jimited fluvial erosion and sedimentary rock formation at T < 273K (Doran 



et al. 2010 Lee and McKay 2003 Marchant and Head, 2007). Weathering 



and mineralization is confined to lakes, hyporheic zones, and a shallow active 
layer. However, seasonal river discharges reach 20 m^ s~^ (McKnight, 2011), 
fluvial valleys incise > 3m deep into granite (Shaw and Healy, 1980[ ), and 



annually-averaged weathering intensity within the hyporheic zone is greater 



than in temperate latitudes (Nezat et al. , 2001). Ions are concentrated within 



ice-covered lakes by sublimation. Outcrops of gypsum, carbonate evaporites, 
and algal limestone sediments show that sediments have accumulated at the 



base of melt-fed perennial lakes for 300,000 years (McKay et al. , 1985 Hendy 



2000). Dry Valley Drilling Project cores show lithification in older horizons 



(McKelvey, 1981) 



Order-of-magnitude energy and mass balance shows that brief, rare pulses of 
snowmelt provide enough water to form the kilometers of sedimentary rock 
observed on Mars. For solar luminosity reduced by 23%, peak noontime in- 
solation at Mars at perihelion on a moderate-eccentricity (e=0.15) orbit is 
~ 630 W/m^. If the snowpack is dusty then its albedo will be that of Mars 



dust, 0.28 (Putzig et al. , 2005). During melting, radiative losses are (jT^^^ ^ 
320 W/m^, and for a 200 mbar atmosphere a reasonable value for wind-speed 
dependent sublimation losses into dry air is ~60 W/m^. Conductive losses will 
be roughly one-half the diurnal temperature range divided by the diurnal skin 



depth, giving 60 W/m for the snowpack material properties in Carr and Head 



(|2003j) and a lOOK diurnal cycle of surface temperature. Greenhouse forcing 
from a 200 mbar CO2 atmosphere equilibrated with a 230K daily mean surface 
temperature is ~ 60 W/m^(from detailed radiative transfer calculations: Ap- 
pendix [B]). Neglecting all other gains and losses, the net energy available for 
melting is therefore 630(1-0.28) - 320 - 60 - 60 + 60 ~ 100 W/m^, equivalent 
to approximately 1 kg/m^/hr snowmelt. The total water required to form the 
800m-thick Meridiani sediments depends on the water/rock mass ratio (W/R) 



during alteration. W/R is given as <1 by Berger et al. (2009) and <300 by 
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251 



Hurowitz and McLennan (2007). This corresponds to a time-integrated melt 

252 column of either ^0.3 km or ^100 km, respectively, for a bulk Meridiani sand- 

253 stone density of 2.3, and ignoring the contribution of water bound in hydrated 

254 minerals to the solid mass of the deposit. Assuming melting occurs for 10% of 

255 each sol, the melt season lasts for 10% of the year, and 10% of melt is avail- 

256 able for alteration, snowmelt production for only 30 Kyr (for W/R = 1) or 10 

257 Myr (for W/R = 300) provides enough snowmelt to reach the upper limit on 

258 W/R for the entire Meridiani sandstone. The period in question lasted O(IO^) 

259 years, so climate and orbital conditions favorable for surface liquid water at 

260 Meridiani are only needed for <1% of the time. This is consistent with the 

261 wet-pass filter sketched in Figure [TJ 



262 3 Distribution of sedimentary rocks on Mars 
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The Mars Orbiter Camera Narrow Angle (MOC NA) team documented ~4,000 



layered rock outcrops of probable or likely sedimentary origin" (Malin and 



1 

Edgett 


2000 


2001 


|Malin et al. 


2010 



are give in Appendix |Xj The resulting distribution of sedimentary rocks on 
Mars (Figure |2]) suggests that surface water availability was narrowly concen- 
trated near the equator and at low elevations. 64% of sedimentary rocks are 
within 10° of the equator, 60% when the Valles Marineris region is excluded 
(Figure [2^). Blanketing by young mantling deposits may contribute to the 
paucity of sedimentary rocks poleward of 35°, but cannot explain the rarity of 
sedimentary rocks at 10-35° latitude relative to the equatorial belt. The ±10° 
band is not unusual in thermal inertia, dust cover index, albedo, or surface age 
distribution, so a dependence of sedimentary rock on these parameters could 
not explain the latitudinal distribution. 

On average, sedimentary rocks are lower than ancient terrain by 2km (Figure 
[2]d). On Earth sedimentary rocks are low-lying because of sediment transport 
by regional-integrated channel networks. Evidence for regionally-integrated 



jghannel networks on Mars mostly predates the sedimentary rock era ( Carr and 



Head 2010 Fassett and Head 2011). The low-elevation bias is independent of 
the equatorial concentration. Therefore, the low-elevation bias is reflective of 
a planetwide, non-fiuvial process that occurs preferentially at low elevations. 

Sedimentary rock abundance does not decline monotonically away from the 
equator. Abundance away from the equator is much less than in the equatorial 
sedimentary-rock belt, but "wings" of increased sedimentary rock abundance 
are found nearly symmetric about the equator at 25-30S and 20-30N. The 
sedimentary rocks in the southern wing are regionally associated with clusters 
of large alluvial fans (Figure [2): 
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Fig. 2. (a) Latitudinal dependence of sedimentary rocks, masking out recently resur- 
faced terrain (Appendix [A|. Latitude bin size is 10°. Histogram corresponds to 
number of images (left axis). Dark gray bars are the contribution from the Valles 
Marineris region, and black bars represent the rest of the planet. Lines correspond to 
the percentage of images showing sedimentary rock (right axis). Dashed line is the 
percentage of images showing sedimentary rocks once the Valles Marineris region 
is excluded, (b) Elevation dependence of sedimentary rocks, masking out recently 
resurfaced terrain. Elevation bin size is 500m. Gray line is normalized histogram of 
terrain with sedimentary rocks, and black line is histogram of all ancient terrain. 
Dotted gray line is the normalized histogram of terrain with sedimentary rocks, 
after masking out Valles Marineris. Median sedimentary rock elevation is ~2km 
lower than median ancient terrain, (c) Distribution of sedimentary rocks on Mars 



(black dots, from Malin et al. (2010)). Alluvial fans are also shown (purple dots, 
from Kraal et al. (2008)). Blue horizontal lines highlight the ±10° latitude band. 



Dark gray shading corresponds to recently resurfaced terrain. Light gray contours 
show topography, with the +10km contour shown in red. 



289 4 Model 



290 This section describes the seasonal melt model. §4.1 describes the model frame- 

291 work and assumptions. §4.2 describes the ID snowpack thermal model. §4.3 

292 describes the potential-well approximation for warm-season snow locations, 

293 and §4.4 explains how results from ID models are combined to produce pre- 

294 dictions. 
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Assemble Snowpack Energy Balance Look-Up Table 
Using a 1 D Model of the Diurnal Thermal Cycle 



calculate diurnal Tcycle for each season (Ls) 



Step 1: Model annual 
cycle of snowpack 
potential temperature 
on a latitude-season grid 



Step 2: Repeat for all 
orbital conditions 



Step 3: Repeat, altering 
ATand pressure 



Postprocessing: Project onto IVlars Topography 
and Sum Over All Orbital States 



Step 4: Interpolate onto 

topography ^^arylngP ' 
3f datum 




Step 5: Calculate snow 
& melt distribution 



Step 6: For each candidate 
paleoclimate, integrate melt 
predictions over distribution 
of orbital states 







PREDICTIONS: 

Maps of peak melt rate and 
total melt for each 
candidate paleoclimate 




DATA: 

Sedimentary rock 
distribution, geochemistry, 
and paleohydrology. 







Fig. 3. Workflow of the Early Mars seasonal melting model. See text for details. 

295 ^.i Overview of model framework 

296 Controls on Mars snowmelt include:- 
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Spin-orbit properties O = {(f), e, Lp, L^, latitude} control the distribution of 
sunlight at a given location. These include obliquity 0, eccentricity e, solar 
longitude of perihelion Lp, solar longitude Ls, and latitude. Milankovitch 
parameters O' = {0, e, Lp} oscillate or circulate on 10"^"^ yr frequencies, 
and shows chaotic shifts at ~250 Myr intervals (Head, 2011). We iterate 
over all the spin-orbit properties that have probably been encountered by 
Mars over the last 3.5 Ga (Steps 1-2 in Figure |3|. 

Climate parameters C = {P, AT, fsnow} include atmospheric pressure P 
(assumed to be mostly CO2), freezing-point depression/non-C02 greenhouse 
forcing AT, and snow coverage fraction fsnow These are iterated over a 
large range (Step 3 in Figure [3]). They are assumed to vary slowly relative 
to changes in spin-orbit properties. 

Surface material properties, insolation. These are held fixed for the ensemble 
of model runs. Results are sensitive to snowpack TI and albedo. Parameter 



choices and sensitivity tests are discussed in Appendix D 



The predictions of the snowmelt hypothesis (^6) are evaluated as follows 
(Figure [3]). First, for a given trial set of past climate parameters C and orbital 
parameters O, snow temperature for all seasons and latitudes is calculated 
using a ID surface energy balance model (Step 1 in Figure |3]). Using the ID 
model output for a range of P (Step 3), the potential annual-average snow 
sublimation rates are mapped onto topography (Step 4). Warm-season snow 
is assigned to locations most favorable for interanually persistent snow (Fig- 
ure 10). Snowmelt occurs when temperatures at these locations exceed freezing 



(Step 5). These results provide latitude-season diagrams that are appropriate 
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321 



for a flat Mars. Results for this fictitious "cueball" Mars are analyzed in ^3. 

322 To map the results onto real Mars topography, a sequence of latitude-season 

323 output for different P but the same O' is stacked in elevation (results given in 

324 ^4). The output at this stage consists of maps of snow stability for the given 

325 C and O', along with time series of snow temperature and melt rates. Next, 

326 the framework loops over all possible Early Mars O (Step 6), convolving the 



327 outputs with the O' probability distribution function (Laskar et al. , 2004). 

328 The output is now a map of predicted snowmelt on Mars integrated over geo- 

329 logic time, and this can be compared to observed sedimentary rock abundance 

330 and thickness data (^5). These maps are computed for many plausible C. 

331 Assuming the snowmelt model is correct, the C that produces the best corre- 

332 lation between model predictions and data is the best-fit Early Mars climate 

333 (^1). Interpolating in the melt rate output gives a predicted time series at 

334 Gale Crater (Q. 



335 4-2 Thermal model 



Surface hquid water is always unstable to evaporation on a desert planet 
(Richardson and Soto 2008a|[b ). However, transient liquid water can occur 
metastably if temperatures exceed the freezing point, and if P exceeds the 



339 triple point (in order to prevent internal boiling) (Hecht, 2002) 



These conditions are modeled using a ID thermal model (Figure |4]). When 
temperature exceeds (273. 15K - AT), melting occurs and buffers the tem- 
perature at the melting point. Within the snowpack, material properties are 
assumed uniform with depth and heat flow is by conduction and solar absorp- 
tion only. When melt is not present, energy balance in the subsurface layer 
adjacent to the surface is given for unit surface area by (Figure |4]) 



pCsAz—^ = ^(^T^ + LWi +(1 - s,)Q{l)SWi 

radiative terms 
— Sfr — Lfr — Sfo — Lfo (1) 
free convection forced convection 



346 and energy balance at depth z within the snowpack is given by (Figure |4]) 

f)T k rP'T 
pCA.^ = ~^ (l-s,.)Qi.)SWl (2) 

solid— state greenhouse 
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Lfo, Lfr / 

(1-Sr)SWi\ I Sfo, SfrlLWI 



4 J 

Ts — — k = 1 



t 



k = 2 



k = nz 



processes that processes that 

heat surface cool surface 

Fig. 4. Vertical discretization and energy flow in the ID model. Solid horizontal 
lines correspond to solid surface layers numbered K = {1... nz}, dashed horizontal 
lines correspond to atmospheric layers. Tr is the effective atmospheric radiative 
temperature, Ta is the atmospheric surface layer temperature, and Tg is the ground 
surface temperature. The diagonal arrows correspond to energy fluxes: LW ^ for 
greenhouse effect, (1 — Sr)SW \. for insolation attenuated by Rayleigh scattering, 
LW\ for backradiation, and {Lfo,Lfr, Sfo and Sfr} for the turbulent fluxes. Some 
insolation penetrates into the snowpack (dashed continuation of insolation arrow). 
kV^T corresponds to conductive diffusion. 

347 Here, p is snow density, Cg is snow specific heat capacity, Az is the thickness 

348 of the subsurface layer whose upper boundary is the surface, T is the tem- 

349 perature at subsurface level K = {1, 2, ... , n^ } (Figure |4]), k is snow thermal 

350 conductivity, e is the longwave emissivity of ice, LWl is downwelling long- 

351 wave radiation, Sr is the Rayleigh-scattering correction factor, Q{i,2,...,nz} is 

352 the fraction of sunlight absorbed at level z (Appendix [C]), SWl is solar flux 

353 per unit surface area, Sfr corresponds to free sensible heat losses driven by 

354 atmosphere-surface temperature differences, Lfr corresponds to evaporative 

355 cooling by free convection when the atmosphere has relative humidity <1, Sfo 

356 corresponds to forced sensible heat losses caused by cool breezes over warm 

357 ground, and Lfo corresponds to additional evaporative cooling when the wind 

359 are for the 3.5 Gya solar luminosity reported by Bahcall et al. (2001), 23% 

360 below present. 
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361 The ID model draws on previous work by Toon et al. 



Richardson and Mischna ( 2005[) ; WiUiams et aL (2008); Dundas and Byrne 



(jTOSOl); powl (119871); 



i2010); and Liston and Winther (2005). A simpler version is discussed in Kite 
et al. ( 2011b[c ). Representative output is shown in Figure [sj Details of the flux 



365 parameterizations, melt handling, and run conditions are given in Appendix 




4 8 12 16 20 
Mars-hours after local noon 
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Fig. 5. Daily cycle of temperature and melting in the upper 15cm of snowpack. 
Grayscale corresponds to melt fraction at a constant temperature of 273. 15K, with 
contours of 0.1 (edge of gray region corresponds to zero melt). Temperature contours 
are shown wherever the snow is not melting, at intervals of lOK. The blocky outline 
of the "melt = 0" contour corresponds to the 1200s resolution at which model output 
is stored for subsequent processing: the underlying timestep is 12s. Ls = 0° , Lp = 
0°, 6 = 50°. 



367 4-3 Snow location prescription: the potential-well approximation 



368 
369 
370 
371 
372 
373 
374 
375 
376 
377 
378 



For ~100 mbar CO2, the model predicts melt only during the warmest season, 
and usually within a diurnal skin depth of the surface. Experiments, observa- 
tion and theory agree that warm-season snow and ice within this depth range 



on Mars is in cold-trap equilibrium with orbital forcing (Mellon and Jakosky 



1995 




Hudson and Aharonson 




2008 




Hudson et al. 


Aharonson , 


2005 


Boynton et al. 


2002). Because we 



freezing T only where snow is present, the quantity of interest is the annual- 
maximum T experienced by the cold traps, whose location depends on orbital 
conditions and topography. For most orbital conditions, this T is below freez- 
ing, so the greatest interest is in the orbital conditions that maximize the 
cold-traps' annual maximum T. 



379 In the case of no topography (cueball planet), the location of the cold traps 
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380 depends on orbital conditions only. For this case we find the location of cold 

381 traps in the following manner: for each O', the output of the thermal model 

382 for all seasons (Lg) and geographic locations x is used to determine the x 

383 where snow is most likely to be present during the melt season. Melt-season 

384 snow is assumed to be only found at locations where the annually- averaged 

385 sublimation is minimized (Figure [6]). Converged thermal-model output for Lg 

386 is linearly interpolated in time using Kepler's equation (Figure |6]). All x are 

387 then assigned an area-weighted rank, /, scaled from 0% (global minimum 

388 in annually-averaged sublimation; most favorable for snow accumulation) to 

389 100% (global maximum in annually-averaged sublimation; least favorable to 

390 snow accumulation). Ice lost by melting is assumed to be recovered by refreez- 

391 ing close (<100km) to source. Warm-season snow is assumed not to occur 

392 above a critical /, termed fsnow (the percentage of the planet's surface area 

393 that has warm-season snow). Using the /(x) and fsnow, warm-season snow is 

394 assigned to favored geographic locations (Figure [?]). In general, the critical 

395 / for warm-season snow will be greater than the critical / for interannually- 

396 persistent snow, so melting does not have to be supraglacial. The most favor- 

397 able circumstances for aqueous alteration may be where melt occurs during a 

398 seasonal accumulation-ablation cycle which leaves bare soil during part of the 

399 year. 



400 Melting is almost certain to occur when orbital forcing leads to annual-maximum 

401 temperatures above freezing at all latitudes (Figure [?]). Thermal barriers >10 

402 cm thick can insulate snow against diurnal melting, but a sublimation lag cov- 

403 ering all ice is logically impossible, and a debris lag covering all ice is unlikely. 

404 The albedo of pure, fresh, fine-grained snow is high enough to prevent melt- 

405 ing, but contamination with dust is very likely (Appendix [C]). Twice- yearly 

406 transfer of the water ice reservoir across the equator to the cold high-obliquity 

407 winter pole would require unreasonably high seasonally reversing mean wind 

408 speeds. 



With MOLA topography, annual-average sublimation rates are calculated as 
for the cueball case, but now for a range of P that spans mountaintop pressures 
and canyon pressures. The latitude-P grid is then interpolated onto latitude 
and longitude using MOLA topography (Figure 12). This assumes that Mars' 
long-wavelength bedrock topography was in place before the sedimentary rock 
era, in agreement with geodynamic analysis (e.g., Phillips et al. , 2001). This 



also neglects the effect of the adiabatic lapse rate on surface T, which is an 
acceptable approximation at the P of interest here 
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Fig. 6. Seasonal cycle of diurnal-peak temperature and diurnal-mean free sublima- 
tion rate for 3.5 Gya insolation, assuming flat topography. 145 mbar pure CO2 at- 
mosphere. (a,b) present-day orbital forcing {(p = 25.2°, e = 0.093, Lp = 251°); (c,d) 
optimal conditions for melting - high-c/;, moderate e, and Lp aligned with equinox. 
Contours of daily maximum surface temperature are drawn at 180K, 200K, and 
210K and then at intervals of 5K up to a maximum of 270K, only reached in (c). 
White shading corresponds to CO2 condensation at the surface. Sublimation-rate 
contours are drawn intervals of 0.025 kg/m^/sol from to 0.1 kg/m^/sol and then 
at intervals of 0.2 kg/m^/sol. At low e and low (p (similar to today, (a,b)), ice is 
stable at the poles, where temperatures never exceed freezing. In (c,d), ice is most 
stable at the equator, and annual peak temperature exceeds freezing everywhere at 
some point during the year. The thick black line in (c) corresponds to subsurface 
melting at some point during the day. (No melting is predicted for modern orbital 
conditions.) The blockiness of this line corresponds to the underlying seasonal res- 
olution (22.5° in Ls). Solid-state greenhouse effect warms the subsurface relative to 
the surface by up to several K. 

417 4-4 Integrating over all orbital states 



418 For each O', / is obtained for each x, together with the potential snowpack 

419 temperatures and potential melt rates. Then, fsnow maps out the snow distri- 
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Fig. 7. To illustrate the potential-well approximation for finding warm-season snow 
locations. Nominal parameters (Table 1), 146 mbar atmosphere, flat topography. 
Curve corresponds to potential sublimation during a year. Temperatures are annual 
maxima. "MELT" denotes melting of snowpack at some point during the year, if it 
exists at that latitude. Horizontal lines correspond to fsnow values, assuming war- 
m-season snow is found at locations that minimize annually-averaged sublimation, 
and colored red for values of fsnow that lead to melting, (a) Current orbital condi- 
tions. Massive ice or buried ice may exist in the southern hemisphere, but snowpack 
that persists through the warm-season is only likely in the far north, where temper- 
atures are always below freezing. (b) Optimal orbital conditions. Melting occurs at 
all latitudes so melting is inevitable, (c) As for optimal orbital conditions, but for 
perihelion aligned with northern summer solstice. The short, intense northern sum- 
mer displaces the potential-sublimation minimum to 20S. fsnow > 23% is needed 
for melting to occur under these circumstances. The latitude of first melting will be 
near the equator. Note that '273K' is just below the melting point. 

bution, and fsnow and AT map out the melt distribution. The resulting maps 
of snow and snowmelt are weighted using the O' probability distribution func- 



tion (Laskar et al. , 2004), and the weighted maps are summed. (Note that this 



is not the same as the median melt column from a large ensemble of solar 
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424 system integrations; Laskar et al. (2004).) Melt likelihood is then given by 



M^=J {T^a.,. > (273.15 - AT)) (/,„,^ > p(O') dO' (3) 

425 where the "greater than" operator returns 1 if true and if false, and p() is 

426 probability. 



19 



Symbol 



Parameter 



Value and units Source / rationale 



Fixed parameters: 

^vonk 

Cp 

Cs 

9 

kb 

h 

"-snow 

nic 

M„ 

Pat'm,0 
Le 
Le 
Th 

Rgas 

Zo 

Zanem 

a 
t 

p 

Po 
a 



von Karman's constant 
Specific heat, Mars air 
Specific heat capacity of snow 
Mechanical difFusivity of air 
Mars surface gravity 
Boltzmann's constant 

Thermal conductivity of snowpack 

Molar mass of CO2 

Molar mass of H2O 

Molecular mass of H2O 

Pressure of atmosphere now 

Latent heat of water ice sublimation 

Latent heat of water ice melting 

Relative humidity 

Gas constant 

Reference near-surface wind speed 
Roughness length 
Anemometer height 
Albedo 

Emissivity of ice at thermal wavelengths 

Kinematic viscosity of air 

Density of snowpack 

Density of atmosphere (now) 

Stefan-Boltzman constant 

Time of interest 

Mars semimajor axis 

Duration of 1 Mars sol 

Dust concentration 

Dust radius 

Ice grain radius 

Solar constant (now) 



0.4 

770 J/kg/K 
1751 J/kg/K 
14 X 10-* m^/s 
3.7 m/s2 



1.38 X 10" 



m kg s 



5.67 X 10 



0.125 W/m/K 
0.044 kg 
0.018 kg 
2.99 X 10^6 kg 
610 Pa 
2.83 X 10^ J/kg 
3.34 X 105 J/kg 
0.25 

8.3144 J/(mol K) 
3.37 m/s 
0.1 mm 
5.53 m 
K,Q.2'& 
0.98 

6.93 X 10-* m^/s 
350 kg m-3 
0.02 kg m-3 
W m-2 K-* 
3.5 Gyr ago 
1.52366 AU 
88775 s 
~2% by volume 
4 /.im 
1 mm 

1.361 X 10=^ W/m^ 



Carr and Head 


(|2003f 


Hecht \mi\ 






Carr and Head 


(|2003| 



NSSDC 



'Hecht 


(■2002 


Hccht 


mm 



'Millour et al.' |'2008 ll "MY24" a verage 
Polar snow. Brock ct al.| l |2006| 
|Millour et ai., ( ,2U08| 
Albedo of dusty snow. Appendix [P] 



Hccht (2002) 
Carr and Head 

'nssdc 



I 2003 I 



'Murchie et al.|l|2009b| 
InSSOC ' ' ' 

Appendix C 
Appendix C 

Appendix C 

Kopp andTiean I2OII') 



Extrapolation from GCM runs 



Selected variables: 

O Spin-orbit properties 

O' Milankovitch parameters 

C Climate parameters 

4> Obliquity 0-80° 

huB Dundas-Byrne "b" f(^') (Appendix [B|| 

e Eccentricity 0.00 - OrTC 

Lp (Solar) longitude of perihelion — 360° 

Ls Solar longitude - 360° 

M Mean anomaly - 360° 

P Atmospheric pressure 24-293 mbar 

Po Atmospheric pressure at zero elevation 24-293 mbar 

AT Non-C02 greenhouse forcing 0-15 K 

fsnom Fraction of planet surface area with warm-season snow 0-50% 

Qfc Fraction of incident sunlight absorbed at level k 0-100% 

LW^ Greenhouse forcing 

LW't Thermal emission by surface 

SVyj, Insolation 

Sr Rayleigh scattering correction factor 

Lfo Latent heat losses by forced convection 

L-fr Latent heat losses by free convection 

Sfo Sensible heat lost by forced convection 

Sfr Sensible heat lost by free convection 



Table 1 

Selected parameters and variables. 
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428 5 Results and analysis 



429 5.1 Controls on the occurrence of near- surface liquid water on Early Mars 



430 Warm-season snow locations depend on sublimation rates. Diurnal-mean sub- 
limation is shown as a function of P in Figure |8| Losses due to free convection 
decrease with increasing P, because the greater atmospheric density dilutes 
the buoyancy of moist air (Appendix |B| . Losses due to forced convection in- 
crease with atmospheric density (Appendix [B]). Surface temperature increases 

435 monotonically with increasing greenhouse forcing, leading to an uptick in sub- 

436 limation rate for P > 100 mbar. Snow is most stable against sublimation when 

437 P ~ 100 mbar. Therefore, when the atmospheric pressure at zero elevation 



431 
432 
433 
434 



439 



440 



438 (Po) <^ 100 mbar, snow is most stable in topographic lows (Fastook et al. 



2008). When P^ ^ 100 mbar, snow is most stable on mountaintops. 



Melting and runoff depend on energy fluxes around the hottest part of the day. 

441 Figure |9] shows the terms in the energy balance for a snow surface artificially 

442 initialized just below the freezing point. At low P, Lfr exceeds insolation and 

443 melting cannot occur. At high P, Lfr is much less important. Instead, absorbed 

444 insolation and greenhouse warming are balanced principally by radiative losses 

445 and melting. Whether or not surface melting can be sustained will depend on 

446 the partitioning of the subsurface absorbed insolation (Figure [9]). 

447 Setting AT = 5K raises peak melt rate from 0.44 kg/m^/hr to 1 kg/m^/hr. 

448 Melt fraction reaches 1 in the shallow subsurface. Melt is produced for 6 (in- 

449 stead of 4) Mars-hours per sol, and there is some melt in the subsurface for 

450 7.5 (instead of 4) Mars-hours per sol. 

451 Both a and O' affect the energy absorbed by the snowpack. This can be shown 

452 by considering the energy absorbed by equatorial snow at equinox: 



/ l-e^ \" 

Eegmno. ~ (1 " a)LJ ^^^^^^^ 1 (4) 



distance from Sun 

where Ef.quinox is the sunlight absorbed at noon at equinox, is the solar 
luminosity at Mars' semi major axis at geological epoch r, ^ is the minimum 



jugular s eparation between Lp and either Lg = or = 180 (Murray and 



Dermott (2000), their Equation 2.20), and the atmosphere is optically thin. If e 
is large then peak insolation need not occur at equinox. Equation |4] shows that 
moving from average orbital conditions (e = 0.06, \E' = 90) to optimal orbital 
conditions (e = 0.15, \E' = 0) has the same effect on Eequinox as darkening from 
albedo 0.28 to albedo zero. 
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Pressure (Pa) 

Fig. 8. Pressure dependence of sublimation rate at an equatorial site. FREE is Ljr, 
FORCED is Lfo- Solid lines with asterisks correspond to a wind speed that declines 
with increasing P, dotted lines correspond to constant near-surface wind speed of 
3.37 m/s. e=0.11, a ^0.28, Lp=0, L,=0, (/>=50. 

461 5.2 Seasonal cycle and snow locations 



462 Figure [6] shows the seasonal cycle of T and sublimation rate on a flat planet. 

463 Annual average sublimation rate controls warm-season snow location, and 

464 annual-peak snow temperature determines whether melting will occur. The 

465 cold trap latitudes indicated correspond to fsnow — ^ 0, i.e. a single thick ice- 

466 sheet. Suppose instead that warm-season snow covers a wider area - that the 

467 "potential well" of Figure [7] fills up with snow. For modern orbital conditions, 

468 the area of snow stability will then extend south from the North Pole. If warm- 

469 season snow covers more than 43% of the planet - if the cold-trap effect is weak 

470 or nonexistent - then melt is possible even under modern orbital conditions. 

471 For optimal orbital conditions, increasing fsnow spreads the melting area to 

472 form a broad band equatorward of 30°. 



473 5.3 Flat Mars - snowmelt on a cueball 



474 Next, maps are summed over orbital states to find the latitudinal distribution 

475 of surface liquid water at a given P. If liquid water supply limits sedimentary 
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Fig. 9. Pressure dependence of the surface energy balance for an equatorial site with 
Ti = 273. 15K imposed. Wind speed declines with increasing P. Note that subsurface 
melting can occur for Tg^rf < 273. 15K. The absorbed insolation terms decrease 
slightly at high P due to Rayleigh scattering. Greenhouse forcing is stronger than 
in the time-dependent case shown in Figure [5] because the atmospheric temperature 
is assumed to be in equilibrium with a surface at the freezing point. If 1/2 of the 
subsurface absorbed insolation returns to the surface through conductive heating, 
then melting will occur (net cooling at the melting point will be zero) for P > 130 
mbar. However, if none of the subsurface absorbed insolation conductively warms 
the surface, then surface melting will not occur even for the highest P shown (~ 
330 mbar). Sjo and Ljo are weak because of the low surface roughness. e=0.15, 
a^0.28,Lp=0° ,Ls=0° , 0=50°. 

476 rock formation, then the sum over orbital states should correspond to geologic 

477 observations of Early Mars deposits. 

478 Suppose that the planet is fiat, with warm-season snow only present in a 

479 narrow ring at the latitude that minimizes the annually-averaged sublimation 

480 rate. Suppose that we impose climate conditions only just allow melting under 

481 the optimal orbital conditions. Then:- 



What is the latitudinal distribution of snow, melt and melt intensity? Obliquity 
is the strongest control on Mars snowpack stability. The ID model predicts 
that snow is most stable near the equator for >40°, near the poles at = 
{0°,10°,20°}, and at intermediate latitudes (±55°) for = 30°, in agreement 



486 _with other studies (Jakosky and Carr, 1985 Mischna et al. , 2003; Levrard 



et al. 


, 2004 


Forget et al. 


, 2006 



to the hemisphere in which aphehon occurs during summer. Holding e fixed, 
the width of the latitudinal belt swept out by warm-season snow during a 
precession cycle decreases with increasing 0, from ±26° at (^=40° to ±6° at 
0=80° (for e=0.09). Holding cj) fixed, the width of the warm-season snow belt 
increases with increasing e, from ±18° at e=0.09 to ±22° at e=0.16 (for cf) = 
40°). 



494 
495 
496 
497 
498 
499 
500 
501 
502 
503 
504 
505 



of melting occurs for latitudes <10°. Annual column snowmelt is sharply 



concentrated at the equator within the thin melt band (Figure 10). Even 
though the probability of a melt year is just ~0.05% , the orbitally integrated 
expectation for the equatorial snowmelt column is 5 km/Gyr, which is the 
global spatial maximum on this fiat planet. A typical Mars-year of melting 
produces 9 kg/m^ melt at the latitude of warm-season snow. Peak instanta- 
neous melt rate is of order 0.1 kg/m^/hr, and is less sharply concentrated at 



the equator than annual column snowmelt (Figure 10). Runoff is unlikely for 
these low melt rates, so any aqueous alteration must result from infiltration or 
from alteration of silicates within the snow. Low rates of infiltration may be 
sufficient to alter aeolian deposits, which volumetrically dominate the Mars 



sedimentary rock record (Grotzinger and Milliken, 2012). 



What is the distribution of melt and melt intensity with orbital conditions? As 
C is moved towards conditions that allow melt, melt first occurs at >40° 
and e>0.13 (Figure 11). Kite et al. ( 2011b|[c ) explain that high is needed to 
drive snow to the equator, and high e is needed to bring perihelion close to the 



Sun (Figure 11). Melting requires that perihelion occurs when the noontime 
sun is high in the sky - for the equator, this requires Lp ~ 0° or Lp ~ 180°). 
Holding e and fixed and moving Mars through a precession cycle, the cueball 
planet is entirely dry between Lp=45° and Lp=135° inclusive, and between 
Lp=225° and Lp=315° inclusive. 



515 What is the seasonal distribution of melt? All melting occurs near perihelion 

516 equinox. The melt season lasts <50 sols (for e < 0.145). 



How does increasing fsnow affect the results? Pinning snow to ±1° of the opti- 
mum latitude corresponds to fsnow ~ 1%, similar to the present-day value for 
warm-season snow. However, midlatitude pedestal crater heights (44±22.5m; 
Kadish et al. (2010)) are geologic evidence for high fsnow in the Late Ama- 



zonian. The modern surface ice reservoir is 2.9±0.3 x lO^km^ (Selvans et al 



2010 Plant et al. 2007). Assuming secular loss of ice is slow, and that the 
midlatitude pedestal craters correspond to the thickness of a single ancient ice 
layer of uniform thickness, the ice accumulation area is ~46% of the planet's 
surface area. This ice accumulation area must be less than fsnow- At fsnow = 
50% in the model, the melt belt thickens to ±33°, with minor melt activity 



527 around ±50° (Figure 10). Maximum annual column melt and melt rates are 
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Fig. 10. Flat-planet results for a climate that only marginally rises above the freezing 
point on Early Mars. The black solid lines correspond to the orbitally-averaged 
melt column per Mars year (relevant for aqueous alteration). The gray solid lines 
correspond to the peak melt rate experienced at any point in the orbits considered 
(the geomorphically relevant melt rate). The thick lines are for a very small value 
of fsnow (4%), and the thin lines are for fsnow = 50%. The melt events poleward of 
35° for fsnow = 50% correspond to extremely improbable orbital conditions. 

at ±22°, where column melt is 5x greater (98 km/Gyr) and peak melt rates 
are 3x greater than at the equator. Peak melt rates now reach 0.7 mm/hr, so 



530 runoff is conceivable (Figure 10). Although there is still a very strong increase 

531 of melting with increasing e, melting can now occur for any value of Lp and 

532 most e (given >40°). Melt is strongest in the Northern Hemisphere at 15° 

533 < Lp < 45° and 135° < Lp < 165°. It is strongest in the Southern Hemisphere 

534 at 195° < Lp < 225° and 315° < Lp < 345°. The explanation for this behaviour 

535 is that the now-broad belt of warm-season snow never entirely shifts out of 

536 the hot zone. The edge of this belt closest to the perihelion summer pole sees 

537 the sun at zenith shortly before (and shortly after) perihelion solstice equinox. 

538 This dramatic increase in melting does not require any change in greenhouse 

539 forcing or paleopressure, just a change in the way the climate system deposits 

540 snow. 

541 How does pressure change affect snowmelt on a fiat planet? At the highest P 
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542 considered, 293 mbar, low values of fsnow produce a broad band of melting 

543 between ±(15-20)°. There is a secondary peak around ±50°. The maximum 

544 in melt rate snaps away from the equator at fsnow > 35%, and above 35% this 

545 maximum moves to gradually higher latitudes. These patterns are similar at 98 

546 mbar and 149 mbar, although the peak melt rates and expected-seasonal-mean 

547 melt rates are both lower because of the reduced greenhouse effect. At 49 mbar, 

548 melting only occurs for fsnow >35%, and never at equatorial latitudes. Lower 

549 P further restricts melting to high fsnow and higher latitudes. Physically, these 

550 trends correspond to the need for a long day-length (or perpetual sunlight) 

551 to warm the snowpack at low P. This is not possible at the equator where 

552 the day is always ~12 hours long. The melt rate under orbital conditions that 

553 are optimal for melting can be thought of as a potential well in latitude, with 

554 maxima at high latitudes (for high-0 polar summer), and a minimum near 

555 the equator. At low P the melt potential is zero at low latitudes, so large 

556 values of fsnow are needed for melting, which will then occur away from the 

557 equator. Increasing P favors melting, so melt rates increase everywhere, and 

558 for 98 mbar and above the melt potential is nonzero even at the equator, 

559 allowing melting as fsnow nears 0. If Mars had a relatively thick 0(100 mbar) 

560 atmosphere during the sedimentary-rock era, then the progressive inhibition 

561 of equatorial melting is a latitudinal tracer of Mars atmospheric escape. 



The "cueball" results presented in the above sections have strong echoes in 
the geologic record of ancient Mars. The distribution of evidence for liquid 
water shows strong latitudinal banding: sedimentary rocks are concentrated 
at latitudes <15° (§3.1) but have 



wings" at 25-30S and 20-30N, and alluvial 



566 fans are most common at 15S-30S (Kraal et al. , 2008 Wilson et al. , 2012). 



567 5.4 With MOLA topography 



568 



570 



576 



The main meridional trends in model output without topography hold when 
MOLA topography is used. This is because nominal model parameters produce 
snow distributions that are more sensitive to latitude than to elevation. 

At Po = 48 mbar and optimal O' (0=50°, e=0.16, Lp = 0°), snow distribution 
shows only a weak preference for low points. However, peak temperature is 
controlled by elevation because Lfr is stronger (and LW I weaker) at high 
elevations. Therefore melting only occurs in planetary topographic lows close 



to the equator at low fsnow and AT=2K (Figure 12). 



For fsnow = 50%, melting at optimal orbital conditions occurs for all low- 
lying locations equatorward of 30°. When perihelion is aligned with solstice 
at 0=50°, the snow distribution shifts away from the equator, and no melting 
occurs below fsnow ~ 60%. 
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Fig. 11. The sensitivity of annual-peak snowpack temperature to Milankovitch forc- 
ing on an idealized Mars lacking topography. Maximum snowpack temperatures 
over a precession cycle (black contours) are highest for high obliquity and moderate 
eccentricity. The probable range of Mars orbital elements (color ramp, with white 
shading least probable and red shading most probable) is much broader than that 
of Earths orbital elements (Gyr range shown by blue rectangle). Black diamond 
corresponds to Mars' present-day orbital elements. Vertical dashed line divides 
<40° (for which warm-season snow is generally found at high latitude), from (/) >40° 
(for which warm-season snow is generally found at low latitude). AT = 5K, P = 
145 mbar, a wO.28, Faint Young Sun. 



580 Low (j) is much less favorable for snow melting than high (p ( Jakosky and Carr 



1985). For < 30° and perihelion aligned with solstice, snow is most stable 
poleward of 60°, but these most favored locations never reach the freezing 
point. As fsnow is raised, melting will first occur at lower latitudes because 
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Fig. 12. A snapshot of snowmelt distribution for a single example of orbital forcing, 
showing role of snow stability and melting potential. (p=50°, e=0.145, Lp = 0°, Pq 
= 49 mbar, AT = 5K. Areas equatorward of the magenta line correspond to fsnow 
< 20% - likely snowpack locations. Areas poleward of the blue line are hot enough 
to see melting at some point during the year, if snow were present. Where the snow 
zone intersects the hot zone, some melting will occur. Melt zones for small values 
of fsnow are shaded in warm colors, fsnow < 10% is shaded yellow, fsnow < 5% 
is shaded red, and fsnow < 2% is shaded maroon. Notice that extensive melting 
in Valles Marineris requires high fsnow or a different phase of the precession cycle. 
Thick black line corresponds to the boundary of terrain resurfaced since sedimentary 
rocks formed. This terrain is not included in the warm-shaded areas. Landing sites 
of long-range rovers are shown by green circles: - Ga = Gale Crater; Gu = Gusev 
Crater; MP = Meridiani Planum. Grayscale contours in background are topographic 
contours at intervals of 1.5 km from -5 km up to +10km. 



these receive more sunlight. The most favored locations are S Hellas and the 
lowest ground around 40N. These are the midlatitude locations where scal- 



586 loped depressions are most prominent ( Soare et al. , 2007 Zanetti et al. , 2010 ) 



587 although these features might not require liquid water to form (Lefort et al. 



2009 2010) and the model is not directly applicable to Upper Amazonian 



features. 



590 Snow and melt distributions on MOLA topography depend on the trade-off 

591 between P and sublimation rate, which controls snow stability (Figure |8]). For 

592 example, suppose wind speed on Early Mars was much higher than modelled. 

593 Then the relative importance of wind-speed-dependent turbulent losses in the 

594 surface energy balance would increase. This would increase the importance of 

595 elevation (~1/P) in setting snow location, relative to latitude which sets SW^. 

596 The snow and melt distributions for this "windy early Mars" (not shown) are 

597 broader in latitude and more concentrated in low areas (especially Northern 

598 Hellas, but also Northern Argyre and the Uzboi-Ladon-Margaritifer corridor). 
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599 5. 5 Effects of climate on snowmelt distribution with MOLA topography 



600 Summing the orbital-snapshot maps of melt likelihood (/p(O)dO) shows the 

601 effect of C = {P, AT, fsnow} on melt likelihood averaged over geological time. 



For Po = 49 mbar and for small values of AT (5K) and fsnow (2%), warm- 
season snow is found primarily in (Figure 13 1) Valles Marineris, the circum- 
Chryse chaos, the Uzboi-Ladon-Margaritifer corridor, craters in W Arabia 
Terra, the Isidis rim, northern Hellas, Gale, Aeolis-Zephyria Planum, and parts 
of the Medusae Fossae Formation, as well as at high (>50°) latitudes. However, 
warm-season snow only melts very close to the equator (Figure [13]d) - in 
Gale, the circum-Chryse chaos, Meridiani Planum, Aeolis-Zephyria Planum, 
the Isidis rim, and the floors of the Valles Marineris canyons. Even in central 
Valles Marineris, among the wettest parts of the planet under this climate, 
melting occurs with probability <0.5% (e.g., 5 Myr of melt years during 1 
Gyr). As fsnow is increased to 5-10% at AT = 5K, melting in Meridiani Planum 
and Valles Marineris becomes more frequent. Melting in Northern Hellas does 
not occur until either fsnow or AT is greatly increased. 



As the atmosphere is lost, melting becomes restricted in space as well as 
time (AT = 5K, Po = 24 mbar, fsnow = 0.1%, Figure [T4)i) . The last hold- 
outs for surface liquid water on Mars are Gale Crater, du Martheray Crater, 
and Nicholson Crater in the west-of-Tharsis hemisphere, and the floors of the 
Valles Marineris canyons in the east-of-Tharsis hemisphere (Figure 14 a,). Gale 
Crater (near 6S, 135E) is usually a hemispheric maximum in snowmelt for 
marginal-melting climates. Melting can only occur for very improbable orbital 
combinations under this climate. If they occurred at all, wet periods would be 
separated by long dry intervals. 



624 At Po = 293 mbar and low fsnow, low-latitude snow is restricted to high ground 

625 and so is melt. Figure [T4]d shows the melt distribution for fsnow = 10% and 

626 AT = 7.5K. For fsnow > 20%, snow is still most likely at high ground, but the 

627 melt pattern flips: melt occurs at all elevations, but it is most common at low 

628 ground as in the low-P case. As fsnow — ^ 100%, melt extent is limited only 

629 by temperature. This is maximized at low elevations, because of the greater 

630 column thickness of greenhouse gas (and because of the adiabatic atmospheric 

631 temperature lapse rate, neglected here). 



Hellas is usually the most favored area for snowmelt within the midlatitude 
ancient terrain. A climate relatively favorable for melting in Northern Hellas 
is shown in Figure [II}:; (AT = 15K, Po = 24 mbar, fsnow = 40%). The contour 
of locally maximal snowmelt extending from deepest Hellas, to the Northern 
Hellas floor, to crater Terby, is intriguing in view of recent descriptions of thick 



ipackages of sedimentary rock in these locations (Wilson et al. 2007 Ansan 



SNOW: maximum 0.33 




MELT: maximum 7.1x10' 




OE 60E 120E 180E 240E 300E OE 

Fig. 13. Probabilities of (upper panel) warm-season snow and (lower panel) melting 
for P = 49 mbar, AT = 5K, and fsnow = 2%. Background is shaded relief MOLA 
topography, illuminated from top left. Maximum probability on the warm-season 
snow map is 0.33, maximum of the melt map is 7.1 x 10"'^ - the location is Hydaspis 
Chaos for both maxima. Contours are at 1%, 5%, 10%, 25%, 50%, 75% and 90% 
of the maximum value. Because melting requires unusual orbital conditions, while 
warm-season low-latitude snow only requires high obliquity, the lowest colored con- 
tour in the snow map is greater than the highest colored contour in the melt map. 
Thick black line corresponds to border of ancient terrain, and grayed-out contours 
are probabilities on recently resurfaced terrain. Long-range rover landing sites are 
shown by red circles:- Ga = Gale Crater; Gu = Gusev Crater; MP = Meridiani 
Planum. 
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is not as favorable for snowmelt. 



641 For the wettest conditions considered (e.g., AT = 15K, Po = 49 mbar, fsnow 

642 = 40%, Figure [l4|i) , melt occurs more than 25% of the time in most places 

643 equatorward of 30°. Such wet global climates grossly overpredict both the 

644 spatial extent of sedimentary rock formation on Early Mars and the extent of 

645 surface aqueous alteration, as discussed in the next section. 
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Fig. 14. Sensitivity of snowmelt maps to extreme variations in model parameters. 
Background is shaded relief MOLA topography, illuminated from top left. Colored 
contours correspond to snowmelt probabilities on ancient terrain. Contours are at 
1%, 5%, 10%, 25%, 50%, 75% and 90% of the maximum melt hkehhood, which is 
given to the top left of each panel. Black line corresponds to border of ancient terrain, 
and grayed-out contours are snowmelt probabilities on recently-resurfaced terrain. 
Long-range rover landing sites are shown by red circles:- Ga = Gale Crater; Gu 
= Gusev Crater; MP = Meridiani Planum, (a) Parameters that only marginally 
allow melting even under optimal orbital conditions: P = 24 mbar, AT = 5K, 
fsnow = 0.1%. (b) High P drives snow (and melt) to high ground: P = 293 mbar, 
AT = 7.5K, fsnow = 10%. This is inconsistent with the observed concentration 
of sedimentary rock at low elevations, but may be relevant to the distribution of 
older valley networks, (c) Parameters that produce snowmelt in Hellas: P = 24 
mbar, AT = 15K, fsnow = 20%. For legibility, the 1% contour is not shown for this 
subfigure. (d) Very high fsnow and AT predict a latitudinally broader distribution 
of sedimentary rocks than observed: P = 49 mbar, AT = 15K, fsnow = 40%. 
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646 6 Snowmelt in space: understanding the distribution of sedimen- 

647 tary rocks on Mars 



648 6.1 Comparison of global data to model output: implications for Early Mars 

649 climate state 



The full climate ensemble consists of 343 orbitally integrated melt-likelihood 



maps similar to those in Figure 14 To reduce this to a manageable number 



for analysis, A;- means clustering was used (Press et al. , 2007). The spatial 
variability of the melt-likelihood maps was normalized by the within-map 
mean and within-map standard deviation, and clustering was carried out on 



655 these self-standardized maps. Representative results are shown in Figure [15 

656 together with the mean melt-likelihood maps for each of the climate clusters 

657 identified. 

658 Dark blue dots correspond to perpetual-global-desert climates. Zero melting 

659 is predicted on horizontal surfaces under all orbital conditions. 

660 At Po > 150 mbar, all AT, and low-to-moderate fsnow, melting occurs at 

661 high elevation (red cluster). Some melt also occurs at mid-southern latitudes. 

662 Neglect of the adiabatic lapse rate will lead to growing inaccuracy at high Pq, 

663 but will not alter the conclusion that warm-season snow will be driven to high 

664 ground at high Pq, far from the places where sedimentary rocks are observed 

665 (Figure llSb). 



For a wide range of Po, all AT, and moderate-to-high fsnow, the framework 

is 



predicts a broad swath of low latitude melting (amber cluster). Figure 15: 
effectively a map of maximum snowpack temperature - as fsnow becomes large, 
warm-season snow is no longer restricted by elevation. Melting is most intense 
at low elevation because of the increased CO2 column, but the overall pattern 
is diffuse in both elevation and latitude. This contrasts with the strongly- 
focussed observed sedimentary rock distribution (Figure |2]). Melt probabilities 
are large over a broad part of the planet, which does not sit easily with thermal 
infrared data indicating that most soil on Mars did not experience volumetri- 



675 cally important aqueous alteration (Bandfield et al. , 2011) 



For Po < 150 mbar and at least one of low AT or low fsnow, the model pre- 
dicts focused, equatorial melting (cyan cluster), in excellent agreement with 



observations (Figure 15 i). The agreement is especially good given the sim- 
plicity of the model physics (§4) and the fact that we are considering one of 
three objectively-defined classes of paleoclimates rather than the optimum C. 
Independently of the snowmelt model, this result suggests that latitude and 
elevation are the main controls on sedimentary rock distribution on Mars, be- 
cause our model physics does not include 3D effects. We highlight seven points 
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of data/model agreement:- 

(1) The thickest sedimentary rock exposures on Mars are in Valles Marineris 



(up to 8km), Gale Crater (5km), and Terby Crater (3km) (Murchie et al 



2009a[ [Anderson and Belli [20lol [Wilson et aLf [20071 |Ansan et aL| |2011D 



Sedimentary layered deposit thicknesses in the chaos source regions are 



689 up to ~lkm (Aram; Glotch and Christensen (2005)). The Medusae Fos- 



690 
691 
692 



sae Formation is a sedimentary accumulation up to 3km thick (Bradley 



et al. , 2002 ) which may also be aqueously cemented sedimentary rock 



in its lower part (Burr et al. 2009, 2010). With the exception of Terby 
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Fig. 15. Effect of climate on Mars sedimentary-rock distribution, assuming a 
snowmelt water source. Top panel shows clustering of climates into three melt-pro- 
ducing classes, plus perpetual-global-desert climates (dark blue dots). Remaining 
panels show maps of the mean of each of the melt-producing climate classes. Border 
colors correspond to the dots in the climate parameter space which contribute to 
that map. "Max" to the top right of each map refers to the spatial maximum in melt 
likelihood, which is the probability that a given location sees some melting during 
the year. The colored contours correspond to melt likelihoods of 5% , 10%, 25%, 
50%, 75%, and 90% of the spatial maximum for that climate class. Thick black line 
shows the boundary between ancient terrain and recently-resurfaced terrain. 



this is the same set of locations where the focussed, equatorial melt- 
ing paleoclimate class predicts maxima in orb itally- integrated snowmelt. 
The Northern Valles Marineris canyons contain thicker sedimentary-rock 
mounds than the southern Valles Marineris canyons, and are correspond- 
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(2) 



(3) 



(4) 



ingly more favored for snowmelt in the model. 

Gale is a hemispheric maximum in ancient-terrain sedimentary rock thick- 
ness, and is a hemispheric maximum in ancient-terrain snowmelt in the 
model. 

Snowmelt is strongly focussed in the Valles Marineris, the chaos source 
regions, and Gale. Predicted deposit thickness dies away quickly from 
these regions. 

Meridiani Planum is correctly predicted to be a local maximum within a 
broader wedge-shaped Sinus Meridiani outcrop narrowing and thinning 



(5) 



(6) 



(7) 



to the East (Edgett, 2005 



Hynek and Phillips 2008; Andrews-Hanna 



et al. 2010; Zabrusky et al. 



2012 Wiseman et al. , 2011). The concen 



tration of Western Arabia sedimentary rock in mound-filled craters (e.g. 
Crommelin, Firsoff, Danielson, Trouvelot, and Becquerel) is reproduced 
by the model. Alignment of Meridiani Planum with snowmelt maximum 



implies net True Polar Wander <10° since sediment deposition (Mat 



suyama et al. 


2006 


Perron et al. 


and Manga 


2010 


)• 



Noachian surface liquid water, consistent with data (eg., Jaumann et al 



( |2010[ )). 

In the Northern Plains, deep equatorial craters are commonly modified by 
sedimentary infill (e.g., Nicholson, Reuyl). This correlation is reproduced 
by the model. 

The focussed equatorial melting climate cluster predicts strong enhance- 
ment of melting in Northern Hellas relative to other locations in the 



same latitude band (similar to Figure 14). However, AT >10K is needed 
for non-negligible melting away from the equator, so this longitudinal 
enhancement is diluted in the class-average map and is not visible. A 
secondary enhancement within this southern latitude belt is the Uzboi- 
Ladon-Margaritifer corridor. These longitudinal enhancements match data 



on the distribution of sedimentary rocks and alluvial fans (Kraal et al. 



2008). However, the model underpredicts the thickness of Terby fill, rel- 



ative to the equatorial belt of sedimentary rocks. 
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The model predicts snowmelt in the circum-Chryse region out of proportion to 
the observed sedimentary rock deposits. However, chaos and outflow channels 



continued to form around Chryse through the Early Amazonian ( Warner et al 



2009 Carr and Head, 2010), and would have destroyed sedimentary rocks 
deposited earlier. If supraglacial snowmelt crevassed to the base of the ice 
mass and inflated subglacial lakes, seasonal melting could have contributed 
to chasm flooding and overflow. Sedimentary rocks overly chaos in Aram and 



lani (Glotch and Christensen, 2005 Warner et al. 2011). The model does 



not predict snowmelt at Mawrth, consistent with Mawrth's interpretation as 



a Noachian deposit formed under a earlier climate (McKeown et al. , 2009) 



nor does it predict snowmelt at Terra Sirenum, consistent with the nominally 
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741 Late Noachian age of the inferred paleolake deposits there (Wray et al. , 2011 ). 

742 Finally, the model predicts snowmelt mounds within deep Northern Lowlands 

743 craters that appear little-modified in CTX, such as the crater at 94. 5E, ION. 

744 If the snowmelt hypothesis is correct, then these craters must postdate the 

745 sedimentary-rock era in order to avoid infilling by sedimentary rock. This 

746 prediction can be tested with crater counts on ejecta blankets. 

747 The colors assigned to the climate classes correspond to a hot-to-cold sequence 

748 in Early Mars climate parameter space. The high elevation (red) and broad- 

749 swath (amber) classes have melt likelihoods as high as 0.17. The focussed 

750 equatorial (cyan) class shows much lower melting probabilities (< 0.054) and 

751 is wrapped around the perpetual-global-desert climates (dark blue). To obtain 

752 the distribution of snowmelt that is in best agreement with sedimentary rock 

753 data, at least one of Pq, AT or fsnow must be small. This assumes that sedimen- 

754 tary rock accumulation is proportional to the number of years with snowmelt. 

755 Sedimentary rock formation involves nonlinear and rectifying processes, which 

756 could allow the snowmelt production predicted by the broad-swath-of-low- 

757 latitude melting (amber) class to yield a sedimentary rock distribution con- 

758 sistent with observations. Alternatively, a broad initial sedimentary-rock dis- 

759 tribution could be focussed by aeolian erosion, which would preferentially ef- 
face thin deposits. Therefore, data-model comparison supports the focussed- 
equatorial- melting climate class (cyan), and rules out the high-elevation melt- 

762 ing paleoclimates (red class) and the perpetual global desert (dark blue). It 

763 disfavors, but does not rule out, the broad-swath-of- melting (amber) climate 
class. 



760 
761 



764 



765 In summary, if the seasonal-melting hypothesis is correct. Mars paleoclimate 

766 has left a fingerprint in the sedimentary rock distribution. Sedimentary rocks 

767 are distributed as expected if Mars only marginally permitted snowmelt, even 

768 under near-optimal orbital conditions. The climates that give the best fit to 

769 data predict planets on which the wettest geographic location would have been 

770 dry for >90% of the time. During the sedimentary-rock era. Mars was a dry 

771 place. 



772 6.2 Possible implications for other geologic data: valley networks, chlorides, 

773 and alluvial fans 



Regionally-integrated valley networks record overland flow prior to the sed- 
imentary rock era. We find that Mars valley-network elevation distribution 
is biased high by 600m relative to ancient terrain, although this may reflect 
the generally higher elevation of mid-Noachian (as opposed to Early Hespe- 
rian) outcrop (Hynek et al. , |2010 ). High elevation is the fingerprint of high 
Po (Figure 14 d; Figure ■|l5p). This suggests a geologic record of progressive 
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atmospheric loss:- Po is >100 mbar at valley-network time (to drive snow to 
high ground as suggested by valley network elevations), falls to ~100 mbar 
by sedimentary-rock time (high enough to suppress evaporative cooling, low 
enough to allow sedimentary rock formation at low elevation), and falls further 
to the current situation (6 mbar: Lfr prevents runoff on horizontal surfaces). 
The melt rates predicted by our model with nominal parameters are < mm/hr, 
probably insufficient to form the classical valley networks. Processes that could 
link runoff from snowmelt to the formation of classical highland valley net- 
works include:- (1) a stronger greenhouse effect than considered here, with 
or without the orbital variability considered in this paper; (2) increasing e 
to ~0.22, as can occur transiently during the restructuring of solar system 



orbital architecture predicted by the Nice model (Agnor and Lin, 2012); (3) 
transient darkenings from impact ejecta and ash and transient heating from 
impact ejecta. 



794 
795 
796 
797 
798 
799 
800 
801 



Chloride deposits {n = 634) are generally older than the sedimentary rocks, 
extremely soluble, rare in the equatorial sedimentary rock bracelet, and region- 



ally anticorrelated with sedimentary rock (Osterloo et al. , 2010). This excludes 



an erosional mechanism for the latitudinal distribution of sedimentary rocks. 
One possibility is that dust obscures chlorides at low latitudes. Another pos- 
sibility is that chlorides were dissolved in the equatorial band during the melt 
events that lithified the sedimentary rocks. This would imply that melt rarely 
occurred far from the equator. 
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Peak runoff production during the sedimentary-rock era is constrained to 
~ 0.3 ± 0.2 mm/hr (Irwin et al. , 2005 Jaumann et al. , 2010). Melt produc- 



tion at these rates is possible in the climate-ensemble shown. However, runoff 
production will be some fraction of the melt rate, because of refreezing and 
infiltration. Similar to the case of the classical highland valley networks, addi- 
tional energy could be supplied by transient darkenings (ash-on-snow, ejecta- 
on-snow) or transient volcanogenic warming. An alternative way to maximize 
runoff at low Po is a phase lag between and the position of cold traps (which 



is set by orbital forcing, e.g. Montmessin et al. (2007)) and the position of ice 
deposits. For example, an ice deposit built up at 20S at high while Lp ~ 90° 
may melt if it is not removed by sublimation before Lp swings back to 270°. 
This phase lag contrasts to the snow considered in this paper, which is always 
in equilibrium with orbital forcing. 
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Evidence that the alluvial fans are younger than most sedimentary rocks 



(Grant and Wilson, 2011) is consistent with loss of CO2 over time, because 

5.3). Figure 10 shows that at high fsnow 
being negligible at higher lati- 



low Po suppresses equatorial melting 
on a cueball planet, melt rates peak at ±22° 
tudes and several-fold lower at the equator. The latitude of the wing "peak" 
increases with fsnow, but wings exist for a broad range of moderate-to-high 
fsnow and Po- Wings are observed in the latitudinal distribution of alluvial 
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Fig. 16. Possible Mars Science Laboratory primary mission exploration targets in 
the foothills of the Gale Crater mound. Finely layered sulfate-bearing and phyllosil- 



icate-bearing sedimentary rocks, locally reworked by a channel. See Anderson and 



Bell (2010) for geologic context. Part of HiRISE PSP_009294_1750. Image is ~4500 



m across, illumination is from top left. 



822 fans on Mars (Kraal et al., 2008 Wilson et al., 2012). 



823 7 Snowmelt in time: predictions for MSL at Gale Crater 



824 7.1 Testing snowmelt at Gale Crater 



The base of the Gale mound is a good place to test the snowmelt hypothe- 
sis, because snowmelt is predicted at Gale for most of the paleoclimates that 
permit surface liquid water anywhere on Mars|3] Model predictions are for 
the lower unit of the Gale Crater mound, which is known to contain aqueous 



829 minerals (Milliken et al. , 2010), not the spectroscopically bland upper unit 



(Thomson et al. , 2011) for which much less snowmelt is predicted. We think 



that it is not a coincidence that Gale is a good place to test the snowmelt 
hypothesis. MSL was sent to Gale because it hosts one of the thickest sed- 



^ Of the subset of climate states considered that predict snowmelt anywhere on the 
planet, 66% predict snowmelt at the base of the Gale Crater mound. If we say that 
Gale Crater has a "robustness" of 66%, then >99% of ancient surface area scores 
lower for robustness. In addition, for 55% of climates modeled, the base of the Gale 
Crater mound is in the top 1% of the planet for melt likelihood. If we say that the 
base of the Gale Crater mound has a "maximality" of 55%, then >99.9% of ancient 
surface area scores lower for maximality. 
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imentary rock packages on Mars, with mineralogic and stratigraphic hints 
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of chmate change (MiUiken et al. , 2010). The snowmelt model predicts rel- 



atively abundant snowmelt at Gale, even in a changing Early Mars climate. 
If snowmelt is the limiting factor in sedimentary-rock production, then nat- 
urally Gale would be a place that would sustain sedimentary-rock formation 
for the widest range of climate conditions. 

Hypothesis: We hypothesize that the Gale Crater mound is an accumulation 
of atmospherically-transported sediments pinned in place and subsequently 
reworked by seasonal-meltwater-limited processes. (This hypothesis is adum- 
brated in a unpublished M.S. thesis by Cadieux (2011), and in conference 



abstracts by Cadieux and Kah (2011) and by Niles and Michalski ( 2012[ ).) 
Tests: The snowmelt model predicts:- 

• Wet-dry cycles on orbital timescales, with dry conditions most of the time. 

• Mound-scale geochemistry records a succession of closed systems, not a flow- 
through geochemical reactor. If the fluids responsible for alteration were 
in contact with the atmosphere, as is true of all the ancient waters yet 
sampled by meteorites and rovers (Halevy et al. , 2011), then T < 273K 



implies restriction of diagenesis to perched aquifers within meters of the 
surface (more beneath lakes). The Gale Crater mound is 5km high, so this 
predicts that the Gale Crater mound is a succession of tens-to-thousands of 
closed systems. If on the other hand the layers near the top of the mountain 
were precipitated from groundwater that had flowed from the bottom of the 
mountain, then the mountain is a fiow-through geochemical reactor. Basal 
layers would then be vulnerable to alteration by subsequent upwelling fluids. 
If smectite layers are found between Mg-sulfate layers, this would place a 



tight upper limit on fiow-through aqueous chemistry (Vaniman, 2011). 



Clay/ sulfate transitions correspond to a change in silicate input, not a change 
in global environmental chemistry. At Gale and many other sites on Mars, 
sedimentary rocks transition upsection from irregular to rhythmic bedding 



(Grotzinger and MiUiken, 2012). This suggests a change over time in the 



relative importance of transient darkenings from volcanism and impacts, 
versus orbital forcing. Early on, large explosive eruptions and large impacts 
were more frequent - so many melt events were assisted by regional-to-global 
albedo reduction. As volcanism and impacts declined, darkening events be- 



came less frequent, so eccentricity change (Figure 11) emerged as the key 



regulator of melt events. Therefore, we predict that the phyllosilicate layers 
in the base of the Gale Crater mound were altered in-situ, and are strati- 



graphically associated with impact ejecta (Barnhart and Nimmo, 2011) or 
volcanic ash layers. 

Generally homogenous chemistry and mineralogy on ascending the mound. 
With the exception of these events, the protolith is globally-averaged atmospherically- 
transported sediment, and most alteration is local. This leaves little scope 
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875 
876 



for unmixing of major-element chemistry. 

No Gale-spanning lakes (except immediately after the Gale-forming impact). 



877 1^ Local perennial lakes are possible, as in the Antarctic Dry Valleys (|Doran 

878 



et al.[|1998D . 

879 • Isotopic gradients. Within a unit representing a single identifiable melt 

880 event, isotopic trends will depend on the water loss mechanism. If the wa- 

881 ter evaporated, earlier deposits will be isotopicallsedry lighter (in H and O 

882 isotopes) and later deposits heavier. This is due to the preferential evapo- 

883 ration of light water and will give an O isotope trend similar to that seen 

884 within ALH84001. If, on the other hand, the water froze rather than evap- 

885 orated, later deposits will be lighter or no time dependent trend in isotopic 

886 composition will be observed. By contrast, in a groundwater model, if the 

887 supply of groundwater is ~constant during mineralization, then the isotopic 

888 composition of the evaporating fluid will be some steady-state value, which 

889 would depend on the isotopic composition of the upwelling fluid and the 

890 evaporation rate. Lesser variability is expected within a single deposition 

891 event. 

892 • No organic carbon. Slow, orbitally-paced sedimentation and oscillation be- 

893 tween reducing and oxidizing conditions would disfavor preservation of or- 

894 ganic carbon. 



895 7.2 From snowmelt time series to the Gale Crater stratigraphic logs 



896 Seasonal cycles and runoff. Early in the melt season, melt will percolate ver- 



897 tically and refreeze (Marsh and Woo, 1984). Vertical infiltration of snowmelt 



898 can indurate and cement aeolian dust and sand. Draining and channeliza- 

899 tion of melt will lengthen the lifetime of subsurface melt, especially late 

900 in the melt season. The impermeable ice table constructed by refreezing of 

901 early-season melt favors late-season runoff. Runoff and ponding of snowmelt 

902 in ice-covered lakes requires that water reaches channels before it refreezes. 

903 Once water reaches channels, ice cover protects against further freezing. Be- 

904 cause the daily average temperature is below freezing (in general this is not 

905 a requirement for seasonal-melting models, but it is a feature of the model 

906 output considered here), this requires that drainage times through firn are 

907 <1 sol, in turn requiring high drainage density. Channel deposits with high 

908 drainage density are sometimes seen within the sedimentary rocks of Mars 

909 (e.g., HiRISE ESP_020602_1755 and PSP_007474_1745), and feed into much 

910 larger (and much more frequently preserved) inverted channels. A possible ter- 

911 restrial analog for these processes is the Coastal Thaw Zone of the Antarctic 

912 Dry Valleys. 

913 Milankovitch cycles. Snowmelt predictions are mapped onto sedimentology 

914 and stratigraphy in Figure 17 (compare Figure [T|. Wet-dry cycles with pe- 
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riod ~ 20 Kyr are inevitable unless AT ~ 15K. Early in the wet phase of a 
wet-dry cycle, infiltration can provide water for diagenesis of layers that were 



deposited under dry conditions (Figure 17). As cementation reduces perme- 
ability, infiltration will decline and runoff will be increasingly favored. The 
primary control on temperature cycles is precession, with secondary control 
by ~100 Kyr eccentricity cycles. The "steady accumulation" column in Fig- 
ure [17] shows sedimentological predictions for the case where atmospherically- 
transported sediment is lithified by infiltration of snowmelt. The "wet-pass 
filter /disconformities" column shows the case where rock formation only oc- 
curs during wet intervals. This produces major disconformities. Quasi-periodic 
liquid water availability at Gale will not necessarily produce quasi-periodic 
sedimentology. On Earth, orbitally-paced climate signals are recorded with 
high fidelity by abyssal sediments, but are shredded by fluvial processes and 



are barely det ectable in fiuviodeltaic sediments (Palike et al. 2006 Jerol- 
mack and Paola 2010). 



Sequence stratigraphy on Earth divides stratigraphic packages into periods 
of relative sea-level rise and fall. Snowmelt model output suggests that the 
equivalent of sequence stratigraphy for Mars will involve mound accumulation 
during rare wet periods, and mound degradation (and apron accumulation?) 
during more common dry periods. 



935 8 Discussion 



936 8. 1 Validity of model assumptions 



The snowmelt model assumes that the gross accumulation rate of atmospher- 
ically transported sediment integrated over 4 Ga is enough to build up thick 
sedimentary rocks everywhere on Mars, but that the supply of liquid water 
is the limiting step for net accumulation. Alternative limiting factors include 
sediment availability and preservation/exposure. 



Mahn and Edgett (2000, 2001) and Edgett and Mahn (2002) have suggested 



sedimentary rocks were once much more widespread. In this case, restricted 
exposure today would correspond to blanketing or erosion since sedimentary- 
rock time. 



Sediment availability is unlikely to have limited sedimentary rock formation on 
Mars, because gross deposition rates for atmospherically-transported sediment 



,pn tod ay's Mars (10^ ^ /xm/yr; Arvidson et al. (1979); Geissler (2005) 



eFal] ( [20l0| ); [Johnson et"aL] ( |2003[ ); [Kinch et al.| ( |2007D ; [Drube et al. 



Geissler 



(2010)) 



950 are not much less than past accumulation rates of sedimentary rocks: 20-50 
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Fig. 17. Snowmelt model predictions for seven hundred thousand years at Gale 
Crater. Left time series: Potential temperature of snowpack at Gale. Because the 
Solar System cannot be deterministically reverse integrated to 3.5 Gya, the orbital 
forcing is necessarily fictitious, but it is realistic (Laskar solution 301003BIN_A . N006 
for 73.05-73.75 Mya, but with 0.02 added to the eccentricity). The color scale cor- 
responds to / at Gale. Red is unfavorable for warm-season snow at Gale, and blue 
is most favorable for warm-season snow at Gale. The vertical blue line corresponds 
to the melting threshold. Black highlights intervals of melting at Gale for fsnow = 
10%. AT = 6K, Pq = 49 mbar. Stratigraphic logs: Two end-member stratigraphic 
responses to orbitally-paced wet-dry cycles over ~ 50 Kyr interval. Orange corre- 
sponds to sediment accumulated during dry intervals, and blue corresponds to sed- 
iment accumulated during wet intervals. In the left column, the Gale Crater mound 
accumulates steadily with time, and layers are cemented by infiltration during wet 
intervals. In the right column, both accumulation and diagenesis are restricted to wet 
intervals. Lower panels: Precession cycles of temperature and /. Imperfect cyclicity 
results from varying eccentricity. Perihelion during northern-hemisphere summer is 
especially favorable for snow accumulation at Gale. Gale is dry when perihelion oc- 
curs during southern-hemisphere summer: snow accumulation is unlikely, and any 
snow that does accumulate fails to reach the melting point. 



951 /iin/yr at the best-measured site (Lewis et al. , 2008). Planet-wide sand motion 



952 occurs on Mars (Bridges et al. , 2012). Sand transport and dust lifting on Mars 



953 is sensitive to small increases in P (Newman et al. , 2005). Because Mars has 
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lost atmospheric CO2 over time (Barabash et al. , 2007), gross accumulation 



rates for atmospherically-transported sediment were probably > 10^~^ /im/yr 
at the O(IO^) mbar level required for snowmelt. Present-day reservoirs of air- 
fall sediment are large. For example, dust deposits at Tharsis and E Arabia 



Terra are O(10)m thick (Bridges et al., 2010 Mangold et al. 2009). On Early 



Mars, background dust supply would be supplemented by sediment produced 
during impacts and volcanic eruptions. We do not have good constraints on 
the current surface dust (and sand) budget and how important finite dust 
reservoirs are for the current dust cycle, let alone on Early Mars. Therefore, 
applying current deposition rates to make the argument that sediment avail- 
ability is not a limiting factor is fairly speculative. Neverthless, if fluxes and 
reservoirs were as large in the past as today, then sedimentary rock formation 
would not have been limited by the availability of atmospherically-transported 
sediment. The difficulty then is to pin the sediment in place for >3.3 Gyr, and 
cementation by snowmelt is one mechanism that can resolve this difficulty. 

Spatially varying precipitation is ignored. On Earth, "[ajdvances and retreats 
of glaciers are broadly synchronous" (Cuffey and Patersonl 2010), because 



small changes in Earth x, O', and AT overwhelm regional variations in pre- 
cipitation. This ablation sensitivity is what makes glaciers good dipsticks for 
Earth's paleoclimate. On Mars, recent glaciations have laid down geomorphic 
strips near-parallel to lines of latitude, suggesting that longitudinally variable 
precipitation is less important than insolation in controlling precipitation and 



snowmelt ( 


Kreslavsky and Head 




2000 


2003 




Heldmann and Mellon 




2004 


Neumann et al. , 


2003 


Kadish et a 


. 2009; 


Hauber et al. 


2008 


Fassett et al. 


2010 


). The most recent equatorial 


ice deposits formed at intermediate eleva- 



tions on Tharsis and Terra Sabaea (Forget et al. , 2006 Shean, 2010), associ- 



ated with 3D effects such as orographic precipitation. Models disagree about 



where ice should precipitate under different orbital conditions (Mischna et al 



20031 


Levrard et al. , 


2004; 


Forget et al. 


2006 Mischna and Richardson 


2006 


Madeleine et al. , 


2009 


). This motivates follow-up GCM work. 



Atmospheric collapse to form permanent CO2 ice caps is more likely for Faint 
Young Sun insolation and for ~100 mbar initial P (Kahre et al. , 2011) (For- 
get et al., 2012, submitted manuscript). However, snowmelt requires high (j), 
which is less favorable for atmospheric collapse. Will a CO2 atmosphere that 
has collapsed at low reinfiate on return to high 0? A straightforward calcu- 
lation suggests that atmospheres do not stay collapsed. Dividing a 100 mbar 
atmosphere by the current seasonal CO2 exchange rate of ~3 mbar/yr gives 
a reinfiation time of 30 yr, much shorter than orbital change timescales of 10^ 
yr. Therefore the atmosphere is relatively unlikely to be collapsed for orbital 
conditions that optimize snowmelt. 



994 Neglecting the lapse rate in surface temperature is a good approximation for 



995 ipurrent Mars, where surface temperature is set by radiative fluxes (Zalucha 
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et al.l |2010p. Results from the LMD GCM (IWordsworth et al.l (|2012p; Forget 



et al., submitted manuscript) show that the adiabatic lapse rate is not large 
at 250 mbar but is important for ~ 500 mbar. To cross-check, the Ames 
Mars GCM was run at 80 mbar for modern orbital conditions, topography, 
and luminosity. Only a weak increase in surface-temperature coupling to the 
adiabatic lapse rate was found relative to the 6 mbar case. Therefore, neglect 
of the adiabatic lapse rate coupling to surface temperature appears to be 
adequate for ~ 100 mbar. 

The model assumes that instantaneous values of e, cf) and Lp are independent. 
Reverse integrations of the Solar System (obtained from http : //www. imcce . 
fr/Equipes/ASD/insola/mars/DATA/index.html) show statistically signifi- 
cant correlation (p < 0.0003) between e and 0, but with a very small corre- 
lation coefficient (|R| <0.08) and a sign that varies between integrations. The 
weakness of these correlations justifies treating each orbital parameter inde- 
pendently. Mean probabilities exceed median probabilities for high e, but the 
exceedance probability for e = 0.15 is ~0.8 over 4 Gya (Laskar, 2008). 



We assume C changes more slowly than O, because post-Noachian rates of 
volcanic degassing, weathering, and loss to space are small compared to the 
atmospheric reservoir of CO2. This assumption does not consider volcanic- or 
impact-driven transients in AT. 

We assume the freezing point depression for melting is not very large, which is 
appropriate for sulfates (e.g., AT < 4K for the magnesium sulfate - H2O eu- 



^ectic brine). Chloride brines allow liquid at much lower temperatures (Pollard 
eFaLl fT999| IFairen et al.l [20091). 



8.2 Comparison with other proposed mechanisms for sedimentary rock for- 
mation 



1022 Mechanisms for sedimentary rock formation on Mars must define sources of 



jWater, sediment, sulfur, and heat. In the ice- weathering model of Niles and 



Michalski (2009), the water source is an ice sheet. Sediment and sulfur is 



sourced from dust and gas trapped within the ice sheet. The heat source for 
weathering is the solid state greenhouse effect at shallow depths, and geother- 



ipial heating as the ice is buried. In the global-groundwater model (Andrews- 



Hanna et al. 2007 2010 Andrews-Hanna and Lewis, 2011 ), strong greenhouse 



forcing warms the low latitudes to > 273K (long term average). The water 
source is a deep, regional-to-global groundwater reservoir, which is recharged 
by precipitation or basal melting. Sulfur can be either from pyrite or from the 
atmosphere. The seasonal melting model implies conditions that are warmer 
and wetter than the ice-weathering model, but much colder and drier than the 
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1034 global-groundwater model. Snowmelt under a moderately thicker atmosphere 

1035 is the water source, and insolation under unusual orbital conditions supplies 

1036 heat. Sediment is atmospherically transported - ice nuclei, dust-storm de- 

1037 posits, saltating sand, ash, and fine-grained impact ejecta - and it is trapped 

1038 in the snowmelt area by aqueous cementation. The sulfur source is the atmo- 

1039 sphere. 



The main strength of the ice-weathering model is that it is (near-)uniformitarian. 
Ice-sheet sulfate weathering is ongoing on Earth, and there is evidence for re- 



1042 cent sulfate formation on Mars (Mangold et al. , 2010 Masse et al. , 2012). Cur 



1043 rent gaps in the ice-weathering model include the difficulty of explaining in- 

1044 terbedded runoff features (Grotzinger et al. , 2006), except as post-sulfate re- 



1045 working, and a lack of a physical model for the proposed weathering mecha- 



1046 nism. 
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Global groundwater models can explain the location of sedimentary rocks 



and the diagenetic stratigraphy at Meridiani (Andrews-Hanna et al. , 2007 



Andrews-Hanna and Lewis 2011 Hurowitz et al. , 2010). The global ground- 



water model is internally self-consistent and complete. Upwelling rates are 
consistent with inferred sediment accumulation rates. The discovery of gyp- 
sum veins in material eroded from the Shoemaker Formation ejecta in Endeav- 
our Crater has been interpreted as evidence for bottom-up groundwater flow 



(Squyres and Athena Science Team, 2012). Chaos terrain strongly suggests 



Mars had cooled enough to form a cryosphere that could modulate ground- 
water release. Therefore, even in the groundwater model, post-chaos interior 
layered deposits must have formed via a mechanism consistent with T <273, 



such as spring flow (Pollard et al. 1999 Grasby et al. 2003) 



The advantages of the snowmelt model over previous models for the sedimentary- 
rock water source are as follows. The snowmelt model arises from a self- 
consistent climate solution (§4 - §5), liquid water production can "start and 
stop" rapidly relative to Milankovitch cycles, and the equatorial concentra- 
tion of sedimentary rocks emerges naturally (§5). The snowmelt model can 
account for the global distribution of sedimentary rocks (§6). In the snowmelt 
model, the sedimentary rocks form more or less in their current locations, 
with their current layer orientations, and in their current shapes. Most sedi- 
mentary rocks are now in moat-bounded mounds, flUing craters and canyons. 
Groundwater models imply removal of 3>10® km^ of siliciclastic rock to an 
unknown sink (Zabrusky et al. , 2012 Andrews-Hanna 2012). This removal 



is mediated by a major phase of aeolian erosion which produces the moats. 
Structural deformation is also required to tilt the near-horizontal primary dips 
expected for playa-like deposition to the observed present-day draping dips. 
There is no need to appeal to large-scale postdepositional modiflcation in ei- 
ther the snowmelt model or the ice-weathering model. Notwithstanding these 
advantages, the snowmelt model assumes that precipitation is uniform, but 
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in reality it must have been spatially variable. The snowmelt model also does 
not include a physical model for any of the steps linking melt generation to 
bedrock formation. 



1079 8.3 Atmospheric evolution and the decline of sedimentary rock formation 
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Few sedimentary rocks form on Mars now, and there is minimal surface liquid 
water. The only evidence for surface liquid water at the Opportunity land- 
ing site since the current deflation surface was established is minor Na/Cl- 
enriched veneers and rinds (Knoll et al. , 2008). The simplest explanation for 



these changes is CO2 escape to space. The 2013 MAVEN mission will constrain 
the present-day rate of escape to space. Supposing a 50-150 mbar atmosphere 



at sedimentary-rock time (Figure 15, marginally consistent with Manga et al. 
(2012)), a modern reservoir of 12 mbar (Phillips et al. , 2011), and that soil 



carbonate formation has been unimportant, a loss to space of ~40-140 mbar 
over 3.5 Gya is predicted. Total loss of ~40-140 mbar is higher than previ- 
ous estimates of 0.8-43 mbar over 3.5 Gya from extrapolation of ASPERA-3 

and 2.6-21.5 mbar from fits to MHD 



measurements (Barabash et al 



2007) 



models by Manning et al. (2011). An alternative loss mechanism for CO2 is 



Jiptake by carbonate weathering (Kahn, 1985 Manning et al. , 2006 Boynton 



et al. , 2009 Kite et al. , 2011a). However, many sedimentary rocks contain sul- 



ates, and small amounts of SO2 prevent carbonate precipitation (Bullock and 



Moore, 2007; Halevy and Schrag, 2009). Another alternative is that orbital 



conditions needed to drive melting were sampled early in Mars history, but 
not subsequently. 



1099 8.4 The Early Mars climate trade space 



1100 This paper has emphasized unusual orbital conditions, but more than one 
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1101 mechanism could lead to melting of snow or ice on Mars (Figure 18) 



Deposition of ash or fine-grained impact ejecta can lower albedo, driving tran- 
sient runoff events on Early Mars (Equation |4]). Unweathered silicates provide 
the trigger for their own alteration by darkening the snow. In addition to 
Gale Crater (§7.1), this a-hypothesis is relevant for phyllosilicate formation 
at Mawrth: a regionally extensive, layered deposit which may be consistent 



with top-down alteration of ash (Noe Dobrea et al. 



2009| [Bishop et al.j [2008| |Wray et al.\ [20081 pichalski and N 



2010 



McKeown et al. 



oe Dobrea, 2007 



Loizeau et al. , 2007). Layered clays generally predate the sulfate rocks, consis- 



ientwith decline of volcanism and impacts in the Early Hesperian (Ehlmann 



et al. , 2011). Albedo effects are the primary regulator of spatial and temporal 
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Impact shock energy (10-^ yr per event). 
Impact or volcanically triggered 
global greenhouse ('"lO^ yr per event). 
Very large maximum discharge. 



Base state: 
No runoff. 
Most of the time, 
even on Early Mars? 



decreasing albedo, OL 
WET 

Darkening by ash or fine impact ejecta. 
Wastiout < 10^ yr per event; 
Max. runoff production 2-3 mm/tir 



Unusual orbital conditions favor melting. 
Melt optima '"ICf yr each; imperfect cyclicity. 
Max. runoff production 2-3 mm/hr 



Alternative: Early Mars base state is wet. 

Test: Predicts >10^ yr of continuous wet conditions. 



Fig. 18. The Early Mars climate trade space. Assuming runoff did not occur during 
most years on Early Mars, nonzero runoff can be produced by perturbing orbital 
conditions, reducing albedo, heating from greenhouse forcing or impact shock en- 
ergy, or some combination (black axes). All mechanisms can produce runoff ~1 
mm/hr, but are distinguishable (gray arrows) by their limiting runoff and by their 
timescale. 



in2 (Hall et al. 2010) melt production in the Antarctic Dry Valleys. Experimen- 
1113 tally dusted (13-100 g/m^ fine sand) patches of an Antarctic Dry Valleys 



glacier surface showed sharp increases in melt rate (Lewis, 2001). Antarctic 



1115 Dry Valleys snowpack melts faster when it is buried beneath sand (Heldmann 



2012), provided the sand cover is mm-dm thick. 



Heating from longwave forcing or conduction (AT or impact ejecta heating) 
could also drive melting. Increased LW I could result from clouds (but see 



Polaprete and Toon 


(2003)) 


et al. 


2007 


Johnson et al. 



Melting by impact ejecta is discussed by Mangold et al. (2012 ) and Noe Dobrea 



et al. (2010). 



2008 Tian et al. 2010 Halevy and Head 



2012)) 



1123 Solar luminosity is the final dimension of the Early Mars climate trade space. 
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This paper uses a standard solar model (Bahcall et al. 2001 ) that is consistent 



with solar neutrinos and helioseismology, but not elemental abundances in 



the photosphere (Asplund et al. , 2005). Enhanced mass loss from the young 



Sun would help resolve this discrepancy, and would make the young Sun more 



luminous (Guzik and Mussack 2010 Turck-Chieze et al. , 2011 ). To change the 



conclusion that the Sun was faint at the time the sedimentary rocks formed, 
the Sun's subsequent mass loss rate must have been 2 orders of magnitude 



Jiigher than inf erred from nearby solar-analog stars (Wood et al. 2005 Minton 



and Malhotra, 2007). 



Future work could determine which mechanism is responsible using geologic 
observations that constrain discharge and timescale. Albedo reduction events 
are short-lived, with runoff production that cannot exceed 2-3 mm/hr. Opti- 
mal orbital conditions are relatively long-lived, but again runoff production is 
limited by sunlight energy and cannot exceed 2-3 mm/hr. Volcanic- or impact- 
driven events are short-lived, but with potentially very large discharge. 

The climates considered in this paper are extremely cold, comparable in melt 
production to the coast of Antarctica (Liston and Winther, 2005). These 



climates can produce enough water for aqueous alteration, but struggle to 
match peak-runoff constraints. Therefore, there is scope for exploring warmer 
snowmelt-producing climates, comparable to the coast of Greenland. One pos- 
sibility is that cementation of the sedimentary rocks is the result of opti- 
mal orbital conditions, but that river deposits interbedded with those rocks 



(Williams, 2007) record additional transient (non-orbital) warming. 



1147 9 Summary and conclusions 



The work presented here has two parts. First, a seasonal melting framework 
for Mars has been developed that relates candidate paleoclimate parameters 
to the production of seasonal meltwater. This framework has potentially broad 
applications. Second, the model has been applied to a specific problem: the 
origin and distribution of sedimentary rocks on Mars. 

Seasonal melting on Mars is the product of tides of light and tides of ice, 
which move around the planet on Milankovitch frequencies. The peaks of 
these tides rarely intersect. When they do, melting occurs. This water source 
may contribute to sedimentary rock formation. The main conclusions from 
this work are as foUows:- 

• Order-of-magnitude calculations indicate that snowmelt is a sufficient water 
source for sedimentary rock formation on a cold Early Mars. 

• The distribution of sedimentary rocks on Mars is narrowly concentrated at 
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1161 equatorial latitudes and at low elevations. 

1162 • The optimal spin-orbital conditions for snowmelt in cold traps on Mars are 

1163 high obliquity, longitude of perihelion aligned with equinox, and eccentric- 

1164 ity as high as possible. Melting then occurs in the early afternoon, at the 

1165 equator, during perihelion equinox season. 

1166 • A model of snowmelt on Early Mars has been presented, which uses a 

1167 potential-well approximation to track cold traps for all orbital conditions. 

1168 Integrated over all orbital conditions on an idealized flat planet, and assum- 

1169 ing snowpack with the albedo of dust and a ~100 mbar pure CO2 atmo- 

1170 sphere, the model predicts a narrow equatorial concentration of snowmelt is 

1171 predicted if warm-season snow is tightly confined to cold traps. A broad low- 

1172 latitude belt of snowmelt is predicted if warm-season snow is more broadly 

1173 dispersed. 

1174 • With MOLA topography, atmospheric pressure > 100 mbar drives snow 

1175 to high ground. High fsnow allows snowmelt on low ground even at high 

1176 P. Sedimentary rocks are not on high ground, so cither fsnow was high, 

1177 snowmelt was not the water source for the sedimentary rocks, or P was 

1178 ^100 mbar at time of sedimentary rocks. 

1179 • With MOLA topography, a large swathe of parameter space produces a 

1180 snowmelt distribution that is a good match to sedimentary rock locations. 

1181 Enough water is produced to satisfy mass balance for aqueous alteration of 

1182 sedimentary rock. 

1183 • Early Mars climate states that produce the best fit to the spatial distribution 

1184 of sedimentary rocks on Mars are cold. Much warmer climates would lead to 

1185 snowmelt over a large swath of the planet, inconsistent with observations. 

1186 • Climates that permit surface liqiiid water on Mars usually predict snowmelt 

1187 at Gale Crater. Therefore, if MSL does not find evidence for a snowmelt 

1188 contribution to sedimentary rock formation at Gale Crater, this would be a 

1189 decisive failure of the model presented here. 

1190 • Specific predictions for MSL at Gale Crater include generally homogenous 

1191 aqueous chemical processing on ascending the mound, with clay layers cor- 

1192 responding to a change in siliciclastic input, rather than a change in global 

1193 environmental chemistry. The Gale Crater mound should have experienced 

1194 wet/dry orbital cycles, with wet events only during optimal conditions. Ev- 

1195 idence for vertical fluid flow over distances comparable to the height of the 

1196 Gale Crater mound would be a major failure of the model presented here. 

1197 • This is the first physical model to identify Gale Crater as a hemispheric 

1198 maximum for sedimentary rock formation on Mars. The model therefore has 

1199 the potential to relate observations at Gale Crater to global habitability. 
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A Data analysis 



1201 
1202 
1203 
1204 
1205 
1206 
1207 
1208 
1209 
1210 
1211 
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1214 
1215 
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1217 
1218 
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Locations in tlie MOC NA sedimentary rock database are likely to be strongly 
correlated with the true distribution of sedimentary rocks on Mars, even 
though MOC NA did not sample the planet uniformly. MOC NA took 97,000 



images of Mars; 4% showed sedimentary rocks (Malin et al. , 2010) (http:// 



marsjournal.org/contents/2010/ 0001/files/f igurel6.txt) . Although 
MOC NA imaged only 5.5% of Mars' surface (iMahn et all 120101), the Mars 



Reconnaissance Orbiter Context Camera (CTX) has surveyed >75% of the 
planet at comparable resolution to MOC NA (April 2012 Malin Space Science 
Systems press release, http: //www. msss . coni/news/index.php?id=43) and 
has not found large areas of sedimentary rock missed by MOC NA. MOC NA 
targets were selected on a 1-month rolling cycle on the basis of Viking imagery. 



previous MOC images, and the demands of other Mars missions (Malin et al 



2010). Sedimentary rocks were among the highest scientific priorities of the 
MOC NA investigation (IMahn and Edgettl 120001 IMahn et al.l 120101). In the 



same way that oil wells are drilled more frequently in productive basins, there 
is a high density of MOC NA images in areas of sedimentary rocks identified 
early in the mission. Maps of the relative abundance of sedimentary rocks 
show only minor changes when defined using the fraction of MOC NA images 
showing sedimentary rocks within a given spatial bin instead of the absolute 
number of sedimentary rock observations in a given spatial bin. 



1221 
1222 
1223 
1224 
1225 
1226 
1227 
1228 
1229 
1230 



The definition of sedimentary rock used by the MOC NA team excludes at 
least two areas that are sedimentary in origin, the Terra Sirenum drape de- 



posit and a large part of the Medusae Fossae Formation ( Grant et al. , 2010 



Bradley et al. , 2002). However, the Terra Sirenum drape deposit has a distinct 



phyllosilicate-rich mineralogy from the sulphate-bearing sedimentary rocks 



that are the focus of this paper (Ehlmann et al. , 2011), and the entire area of 



the Medusae Fossae Formation is consistently predicted to be a global near- 
maximum in sedimentary rock accumulation by our orbitally-integrated model 
output. Therefore, neither of these omissions from the database is important 
to the data-model comparison. 



Because formation of sulfate-bearing sedimentary rocks peaked in the Hes- 
perian, Terminal Hesperian and Amazonian terrain may conceal underlying 
sedimentary rocks and should be excluded from the analysis. The currently 



1234 available global geological map of Mars (Skinner et al. , 2006) is a digital ren- 



ovation of Viking-era hardcopy maps (Scott and Tanaka 1986 Greeley and 



Guest, 1987 Tanaka and Scott, 1987). Instead of using the old map, the edge 



of young materials was traced using the USGS Mars Global CIS as a base. The 
resulting ("K12")mask covers 45% of the planet, but only 3.5% of the images 
of sedimentary rocks (n = 105). These 3.5% are mostly from the Medusae 
Fossae Formation and the plateaux surrounding Valles Marineris. These rocks 
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appear to represent a late tail in sedimentary rock formation, and so are re- 
tained in the plots of elevation and latitude dependence (Figure 2). Omitting 
them does not change our conclusions. To check that our results do not de- 



pend on the idiosyncrasies of the tracing, Figure 2 from Nimmo and Tanaka 



(2005 ) was georeferenced onto MOLA topography and areas that they mapped 



as Late Hesperian or Amazonian materials (34% of the planet) were traced 
(the "NT05 map"). K12 draws the dichotomy boundary close to the highland 
break-in slope, whereas NT05 draws the boundary near the lowland edge of 
the fretted terrain. Unlike NT05, K12 includes the Medusae Fossae Forma- 
tion as potentially dating from the sedimentary rock era, in virtue of recent 



£esults proving a Hesperian age for large parts of the formation ((Kerber and 



Head 2010 Zimbelman and Scheldt 2012)). The elevation and latitude re 



suits shown in Figure |2] show little change between NT05 and K12 masks. 
Earlier versions of the analysis presented in this paper ( [Kite et al.[ |2011b|[cl ) 
used NT05, and reached unchanged conclusions. NT05 reflects the modern 



understanding of Mars geology better than Skinner et al. (2006), but it was 
intended as a low resolution overview. We believe that K12 is better suited 
than either the Viking-era maps or NT05 for the purpose of masking out 
post-sedimentary-rock terrain. 

Image-center coordinates are assumed to be close to the locations of sedimen- 
tary rocks. Image footprints can be large for orbits early in the MOC NA 
mission, so these phases are excluded. 

The Valles Marineris are a unique tectonic feature containing many sedimen- 
tary rocks. To make sure that conclusions are insensitive to this unique tectonic 
feature, all data was excluded within a large "Valles Marineris box" (260E - 
330E, 20S - 20N). This did not significantly change the latitude or elevation 
dependence. 

The equatorial concentration of sedimentary rocks was previously noted in 
uncorrected data by Lewis (2009[) (unpublished PhD thesis). 



1270 B Details of thermal model 



1271 Radiative terms: A line-by-line radiative transfer model of the atmosphere 



1272 (Halevy et al. , 2009) is used to populate two look-up tables:- LWl as a func- 



1273 tion of Ti and P; and SWl as a function of P and solar zenith angle. The 

1274 radiative transfer model, which for simplicity assumes a clear-sky, pure CO2 

1275 atmosphere with no clouds or dust, is not run to radiative-convective equilib- 

1276 rium. Instead, for each combination of surface P, T, a, and solar zenith angle, 

1277 an atmospheric P-T structure is prescribed and the resulting radiative fluxes 

1278 are calculated. Following the approach of Kasting (1991), the tropospheric 
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lapse rate is dry adiabatic and the stratosphere is approximated as isothermal 
with a temperature of 167 K. A two-stream approximation to the equations 
of diffuse radiative transfer (which accounts for multiple scattering) is solved 
over a wavelength grid with a spectral resolution of 1 cm~^ at frequencies lower 
than 10,000 cm~^ and a spectral resolution of 10 cm~^ at higher frequencies. 
The error induced by this spectral resolution relative to high resolution calcu- 
lations is small compared to the uncertainties in the other model parameters 
2009). The parameterisation of collision-induced absorption is 
Wordsworth et al. (2010), and is based on measurements by 



(Halevy et al. 



the same as in 



Baranov et al. (2004) and calculations by Gruszka and Borysow (1997, 1998). 



The atmospheric temperature profile corresponding to LWl is pinned to the 
diurnal average Ti. Mars' bulk atmospheric radiative relaxation time is ~2 



days at 6 mbar surface pressure (Goody and Belton, 1967 Eckermann et al 



2011), and increases in proportion to atmospheric density. It is assumed to be 



large for the P relevant to melting (>50 mbar). 

Free convective terms: The turbulent flux parameterizations closely follow 



Dundas and Byrne (2010). Sensible heat loss by free convection is: 



Sfr = O.U{T-Ta)ka 



CpUaPa 




(1) 



where is the atmospheric temperature, ka is the atmospheric thermal con- 
ductivity, Cp is specific heat capacity of air, Ua is viscosity of air, pa is density 
of air, g is Mars gravity, and Ap/pa is the difference in density between air 
in equilibrium with the ground and air overlying the surface layer. Ap/pa is 
given by 



Ap _ {nic - m^)esaf(l - r^) 
p rricP 



(2) 



Here, is the molar mass of CO2, is the molar mass of H2O, is the 
relative humidity of the overlying atmosphere, and Cgat is the saturation vapor 
pressure over water ice. The expression for A p assumes that water vapor is a 
minor atmospheric constituent. 



1305 Ta is parameterized as (Dundas and Byrne, 2010) 



rji _ rpb]jBrpl — b]:,B 

-'a ^ min ^ 



(3) 



1306 where T^m is the coldest (nighttime) surface temperature experienced by the 

1307 model, and hoB is the Dundas-Byrne 'b', a fitting parameter. This is an empir- 
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1308 ical model motivated by Viking 2 measurements (Dundas and Byrne, 2010). 



Bdb decreases as P increases, because atmosphere-surface turbulent coupling 
strengthens. bDsiP) is obtained by fitting to the output of GCM runs at 7, 
50, and 80 mbar which employed a version of the NASA Ames Mars GCM 



described in Haberle et al. (1993) and Kahre et al. (2006). Specifically, &d_b(-P) 



is fit to the global and annual average of the temperature difference between 
the surface and the near-surface atmosphere for local times from 11:00-13:00. 

We let 



LeO.UAripaDa 




(4) 



1316 where is the latent heat of evaporation, At] is the difference between at- 

1317 mosphere and surface water mass fractions, and Da is the diffusion coefficient 

1318 of H2O in CO2. 

1319 Forced convective terms: Sensible heat lost by forced convection is given by: 



Sfo = PaCpUsA{Ta - T) 



(5) 



where Ug is the near-surface wind speed. Near-surface winds are controlled 
by planetary boundary layer turbulence which serves to mix the atmosphere 
vertically, so Sfo 7^ is consistent with the assumption of no meridional heat 
transport. The drag coefficient A is given by 



A 



A2 

vonk 



(6) 



1324 where A^onk is von Karman's constant, Zanem. 



1326 
1327 
1328 
1329 
1330 
1331 
1332 
1333 
1334 
1335 
1336 



is anemometer height, and Zo is 



1325 surface roughness. 



Near-surface wind speed Us in the NASA Ames Mars GCM decreases with 
increasing P and decreasing solar luminosity. The four-season average of Eu- 
ropean Mars Climate Database ( "MY24" simulation) globally-averaged near- 



surface wind speeds at the present epoch is 3.37 m/s (Millour et al. , 2008) 



This is extrapolated for P <290 mbar using a logarithmic dependence of Ug on 
P fitted to the global and annual average of Ames Mars GCM model surface 
wind speed for initial pressures of 7, 50 and 80 mbar. Ug is lowered by a factor 
of 1.08 for the Faint Young Sun using the ratio of wind speeds for two 50 mbar 
Ames Mars GCM simulations that differ only in solar luminosity. Simulations 
suggest Us increases with (Haberle et al.l 2003), but this is ignored. Figure 



9 shows the sensitivity of results to Ug = f(P) and Ug 7^ f(P). 
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1337 Latent heat losses by forced convection are given by: 



1370 



Lfo = L^^^Us{esat{^-rh)) (7) 

where My^ is the molecular mass of water, and k is Boltzmann's constant. 
Latent heat fluxes for dirty snow are calculated assuming that the entire ex- 
posed surface area is water ice. Dirt concentrations are at the percent level by 
volume, or less, for all results presented here, so this is acceptable. 

The free and forced fluxes are summed together, rather than considering only 
the dominant term. This matches the functional form of Mars-chamber data 



i Chittenden et al. 2008) and is the standard approach in Mars research (Dun- 



das and Byrne 2010 Williams et al. 2008; Toon et al. , 1980). However, sum- 



1346 ming the terms is an idealization that may overestimate cooling. 

1347 Melt handling: Melt occurs when > (273. 15K - AT). AT is a freezing-point 

1348 depression. It can also be interpreted as any non-C02 warming due to water 

1349 vapor, ice clouds, or SO2, stochastic fluctuations in material properties around 

1350 those assumed in Table 1, or a higher solar luminosity. Additional greenhouse 

1351 warming (freezing point at 273. 15K) implies greater turbulent and LW'^ losses 

1352 at melting than freezing-point depression (freezing point at 273. 15K - AT), 

1353 but AT is small so this difference is ignored here. 

1354 Total melt present and total melt produced are tracked during the sol. Melt is 

1355 not permitted to drain, and the melt fraction is not allowed to affect snowpack 



jnaterial properties except to buffer temperature during refreezing (Liston and 



Winther, 2005) 



Ablation of the snowpack surface by sublimation is not directly tracked. The 
effect on sublimation on snowpack survival is treated indirectly, through the 
potential-well approximation (^3). However, ablation also affects snowpack 
temperature. Implied sublimation rates are ~0.5 mm/sol for conditions favor- 
able to melting. Movement of the snow surface down into the cold snowpack 
corresponds to advection of cold snow upwards (relative to the surface). Snow- 
pack thermal diffusivity is ~2 x 10~^ m^/s. Melting at depths greater than 
~ n/usubi ~ 4 cm may be suppressed by this advective effect. 

Run conditions. Conductive cooling is found by matrix inversion. Vertical res- 
olution is i^2.5mm for nominal parameters, which is 0.033 x the analytic di- 
urnal skin depth. Time resolution is 12s, and the lower boundary condition is 
insulating. 

The initial condition at the surface is slightly cooler than radiative equilib- 
rium, decaying to the energy-weighted diurnal average temperature at depth 
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1372 with an e-folding depth equal to the diurnal skin depth. The model is inte- 

1373 grated forwards in time for several sols using constant seasonal forcing un- 

1374 til the maximum Ti on successive sols has converged (to <0.01K) and the 

1375 diurnal— maximum melt column (if any) has converg ed to <0.018 kg/m^. For 

1376 polar summers, convergence can take an extremely long time as the melt zone 

1377 spreads to cover the entire snowpack, so the integration stops after ~8 sols 

1378 even if the convergence criteria are not met. 

1379 There is no meridional heat transport, seasonal thermal inertia, or CO2 cycle. 

1380 Temperatures are not allowed to fall below the CO2 condensation point. For 

1381 each spatial location, the model is run for many seasons {Lg). The converged 

1382 output is then interpolated on a grid equally spaced in time to recover annual 

1383 means. 

1384 Details about melt- likelihood map construction. Results in this paper are based 

1385 on grids of runs at = {0°, 10°, 20°, 80°}, e = {0, 0.03, 0.06, 0.09, 0.115, 

1386 0.13, 0.145, 0.16}, Lp = {0°, 15°, 30°, 90°} (with mirroring to build up a 

1387 full precession cycle), = {0°, 22.5°, 45°, 337.5°}, and latitude {-90°, -80°, 

1388 90°}, giving 1.5 x 10^ snowpack thermal model runs for each C. Quoted 

1389 results at intermediate values result from interpolation. Statements about C 

1390 are based on interpolation in a grid of runs at P = {4, 8, 16, 24, 48} x 610 

1391 Pa = {24,49,98,146,293} mbar, with AT = OK. AT and fsnow were varied 

1392 in postprocessing. 

1393 To remove longitudinal stripes of high snow probability in the Northern Plains 

1394 that are artifacts of finite model resolution in O' and latitude, the step function 

1395 in {fsnow — f) is replaced by a linear ramp in {fsnow — /)• This is a minor 

1396 adjustment. 



1397 C Snowpack radiative transfer 
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Crystalline water ice is opaque in the thermal infrared, but almost trans- 
parent to visible light. The resulting solid-state greenhouse effect enhances 



snowmelt (Clow, 1987 Brandt and Warren 1993). The purpose of the solid- 



state greenhouse parameterization in this paper is to self-consistently model 
the tradeoff between snowpack broadband albedo {a) and subsurface absorp- 
tion of sunlight. This does not require precisely calculating a as a function of 
dust content, so the model uses simple linear approximations to the radiative 
transfer equations developed for widely-seperated atmospheric aerosols (e.g. 



Kieffer (1990); Calvin et al. (2009)). Although more sophisticated models can 



be employed to take account of aspherical particles, near-field effects, and het- 



erogeneous compositions (e.g. Cull et al. , 2010; Yang et al. , 2002), the lack of 
consensus on their importance leads us to not include them in our algorithm. 
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The solid-state greenhouse parameterization uses the snow radiative transfer 



ipiodel of Brandt and Warren (1993). Ice refractive indices are from Warren 



and Brandt (2008), and are converted to Henyey-Greenstein parameters using 



a standard Mie code following Bohren and Huffman (1983). Mars dust optical 



parameters are calcuted using the refractive indices of Wolff et al. (2006 2009 ) 



An illustration of these parameters for the canonical atmospheric dust sizes is 



shown in Figure 1 of Madeleine et al. (2011 ), but we also employ larger sizes as 



well. The 2000 ASTM Standard Extraterrestrial Spectrum Reference E-490- 
00 is used to describe the wavelength dependence of the direct flux component; 
diffuse flux is neglected as a being a minor perturbation. The young Sun 
was ~100K cooler in the standard solar model. Solar reddening increases a 
by <0.01, so the spectral shift is ignored here. The effect of small amounts 



of meltwater on a. is minor (Warren, 1982) and is also ignored. The effects 



on wavelength-dependent direct-beam semi-infinite albedo (not shown) are 



broadly similar to the idealized "red dust" in Warren and Wiscombe (1980) 



Once optical properties are prescribed, the most important variables are dust 
content, effective dust grain radius, and effective ice grain radius. A given a 
can usually be obtained by several different combinations of these properties. 



The Brandt and Warren (1993) model is used to build a look-up table of 



fractional subsurface absorption as a function of these variables, plus direct- 
beam path length. This length is mapped to depth within soil by multiplying 
by the cosine of the zenith angle. 



The radiative transfer model reproduces the trends found by Clow (1987). The 



Jarger values of the Martian dust single scattering albedo in the optical (Wolff 



and Clancy 2003 Wolff et al. , 2006 2009 ) reduce the amount of melting for 



a given dust concentration. Ice grain size growth is slow in Mars' present day 



polar caps (Kieffer, 1990) but much faster under the near-melting conditions 



that are important for the model presented here. We adopt an effective size 
of 1 mm, corresponding to observed ice-grain radii in hoar layers in Earth 
snowpacks. Not surprisingly, there are no direct measurements of dust content 
in snow on Mars. Dust content in ice has been reported as "a few percent (up to 
at most around 30%)" by volume in the Northern Plains ( Dundas and Byrne[ 
Mol), and ~15% by mass in the SPLD (iZuber et all 120071). We assume ~2% 



dust mass fraction by volume and a dust grain radius of 4/im. 



1444 D Selection of snowpack material properties 



Snow stability and peak temperatures are affected by material properties such 



as Oi and TI. Low (snowlike) TI (Carr and Head, 2003) is used here because 
snow precipitation is also predicted by all General Climate Models (GCMs) 
at high (e.g., iFastook et all (|2008l); iMischna et all (|2003[); iMadeleine et al. 



1449 (2009)). In addition, water ice precipitation was observed on Mars by the 
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Phoenix lander (Whiteway et al. (2011), their Figure 1). High (icehke) TI 
suppresses the diurnal thermal wave and makes melting at the equator much 
more difficult. 
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The impact of parts-per-thousand levels of dust on snowpack albedo and runoff 



is severe (Warren and Wiscombe, 1980; Warren, 1984). Present day observed 



and calculated Mars seasonal H2O snow a is 0.25-0.4 (Vincendon et al. , 2010 



Kereszturi et al. , 2011 ). a. on the South Polar water ice cap is 0.30 (Titus et al. 



2003). To allow melting to create gullies, snowpack surface-layer albedo must 



have been as low as 0.12 (Williams et al. , 2009). Dust storms and dust devils 
occur every year, and caused major changes in regional and global albedo 



Jbetween 1978 and 2000 ( |Geissler| |2005j) and between 2003 and 2007 ( |Putzig 



and Mellon, 2007). Globe-encircling dust storms, which now occur every few 



years, are likely to occur twice every year at high (Haberle et al. , 2003) 



Dust is required to supply ice nuclei for heterogenous nucleation. Therefore, 
it is reasonable to expect snowpack at high to be contaminated with dust. 
Given the likelihood of dust contamination, this paper assumes a = 0.28, the 
same as Mars' light-toned dust continents. This requires a dust concentration 
of 0(1%) (Appendix [C]) . Other likely sources of darkening contaminants on 
Early Mars are volcanic ash and fine-grained impact ejecta. In the words of 



Warren ( 1984 ) "When snow melts, the impurities often tend to collect at the 



surface rather than washing away with the meltwater." Using a low a favors 
melting, which is conservative because the reconstructed paleoclimate will 
involve the smallest change from the current Mars climate that is consistent 
with the geological evidence. 



1474 Figure [T] shows sensitivity to material properties. Increased TI damps the 

1475 diurnal surface temperature cycle. Increasing albedo lowers surface tempera- 

1476 ture for all times of day, especially near noon. The response of the maximum 

1477 temperature within the snowpack (lower loop) is complicated by subsurface 

1478 absorption of sunlight. This occurs at greater depths when dust concentration 

1479 is decreased, even though the snowpack as a whole is more reflective. Because 

1480 buried solar energy cannot easily escape, nighttime subsurface temperatures 

1481 are increased by increasing the albedo. The location of maximum temperature 

1482 moves to steadily greater depths during the night. The energy-burial effect is 

1483 abruptly reversed shortly after dawn (18 hours after noon), when the location 

1484 of maximum temperature returns to the near- surf ace. 
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Fig. 1. Response of snowpack temperatures to adjusting material properties. Color 
corresponds to local hour angle (time after solar noon). The small jumps in the 
temperature loops are interpolation artifacts. e=0.11, 0=50°, and Lp=0°. 



1485 Acknowledgements 



1486 It is a pleasure to thank the following people for their generosity with time, 

1487 ideas and data. We are grateful to Richard Brandt and Steve Warren for shar- 



1488 
1489 



ing the radiation code underlying Brandt and Warren (1993). Aaron Wolf 



provided help with statistics. Joannah Metz supplied discharge calculations 

1490 for the SW Melas Chasma fans. Mikki Osterloo supplied chloride coordinates. 

1491 This work was triggered by discussions with Oded Aharonson, Jeff Andrews- 

1492 Hanna, and Devon Burr. We thank the anonymous reviewers of an earlier 

1493 version of the manuscript, whose comments improved the manuscript. Discus- 

1494 sions with Konstantin Batygin, Bill Cassata, Bill Dietrich, Bethany Ehlmann, 

1495 John Grotzinger, Alex Hayes, Ross Irwin, Vedran Lekic, Alejandro Soto, Ken 

1496 Tanaka and Robin Wordsworth improved the manuscript. E.S.K. is grate- 

1497 ful to Francois Forget and Robin Wordsworth for sharing their Early Mars 

1498 manuscript. We are grateful to the HiRISE team for maintaining a responsive 

1499 public target request program, HiWish, which was useful in this work. 



59 



1500 E.S.K. and M.M. were sTipported by the U.S. taxpayer through NASA Science 

1501 Mission Directorate grants NNX08AN13G, NNX09AN18G, and NNX09AL20G, 

1502 as well as a startup allocation on the NSF Teragrid (TG-EAR100023). Addi- 

1503 tional computing costs were defrayed by a NASA Mars Fundamental Research 

1504 Program grant to Oded Aharonson. 



1505 References 

1506 C. B. Agnor and D. N. C. Lin. On the Migration of Jupiter and Saturn: 

1507 Constraints from Linear Models of Secular Resonant Coupling with the 

1508 Terrestrial Planets. Astrophysical Journal, 745:143, February 2012. doi: 

1509 10.1088/0004-637X/745/2/143. 

1510 R. Amundson, S. Ewing, W. Dietrich, B. Sutter, J. Owen, O. Chadwick, 

1511 K. Nishiizumi, M. Walvoord, and C. McKay. On the in situ aqueous alter- 

1512 ation of soils on Mars. Geochimica et Cosmochimica Acta, 72:3845-3864, 

1513 2008. doi: 10.1016/j.gca.2008.04.038. 

1514 R. B. Anderson and J. F. Bell, IIL Geologic mapping and characterization of 

1515 Gale Crater and implications for its potential as a Mars Science Laboratory 

1516 landing site. International Journal of Mars Science and Exploration, 5: 

1517 76-128, 2010. doi: 10.1555/mars.2010.0004. 

1518 J. C. Andrews-Hanna and K. W. Lewis. Early Mars hydrology: 2. Hydrological 

1519 evolution in the Noachian and Hesperian epochs. Journal of Geophysical 

1520 Research (Planets), 116:E02007, 2011. doi: 10.1029/2010JE003709. 

1521 J. C. Andrews-Hanna, R. J. Phillips, and M. T. Zuber. Mcridiani Planum 

1522 and the global hydrology of Mars. Nature, 446:163-166, March 2007. doi: 

1523 10.1038/nature05594. 

1524 J. C. Andrews-Hanna, M. T. Zuber, R. E. Arvidson, and S. M. Wiseman. 

1525 Early Mars hydrology: Meridiani playa deposits and the sedimentary record 

1526 of Arabia Terra. Journal of Geophysical Research (Planets), 115:E06002, 

1527 2010. doi: 10.1029/2009JE003485. 

1528 J.C. Andrews-Hanna. The formation of Valles Marineris: 3. Trough formation 

1529 through super-isostasy, stress, sedimentation, and subsidence. Journal of 

1530 Geophysical Research (Planets), 2012. doi: 10.1029/2012JE004059. 

1531 V. Ansan, D. Loizeau, N. Mangold, S. Le Mouelic, J. Carter, F. Poulet, G. Dro- 

1532 mart, A. Lucas, J. -P. Bibring, A. Gendrin, B. Gondet, Y. Langevin, P. Mas- 

1533 son, S. Murchie, J. F. Mustard, and G. Neukum. Stratigraphy, mineralogy, 

1534 and origin of layered deposits inside Terby crater. Mars. Icarus, 211:273- 

1535 304, January 2011. doi: 10.1016/j.icarus.2010.09.011. 

1536 R. E. Arvidson, E. Guinness, and S. Lee. Differential aeolian redistribution 

1537 rates on Mars. Nature, 278:533-535, 1979. doi: 10.1038/278533a0. 

1538 R. E. Arvidson, J. F. Bell, P. Bellutta, N. A. Cabrol, J. G. Catalano, J. Cohen, 

1539 L. S. Grumpier, D. J. Des Marais, T. A. Estlin, W. H. Farrand, R. Gellert, 

1540 J. A. Grant, R. N. Greenberger, E. A. Guinness, K. E. Herkenhoff, J. A. 



60 



1541 Herman, K. D. lagnemma, J. R. Johnson, G. Klingelhofer, R. Li, K. A. 

1542 Lichtenberg, S. A. Maxwell, D. W. Ming, R. V. Morris, M. S. Rice, S. W. 

1543 Ruff, A. Shaw, K. L. Siebach, P. A. de Souza, A. W. Stroupe, S. W. Squyres, 

1544 R. J. Sullivan, K. R Talley, J. A. Townsend, A. Wang, J. R. Wright, and 

1545 A. S. Yen. Spirit Mars Rover Mission: Overview and selected results from 

1546 the northern Home Plate Winter Haven to the side of Scamandcr crater. 

1547 Journal of Geophysical Research (Planets), 115:E00F03, 2010. doi: 10.1029/ 

1548 2010JE003633. 

1549 M. Asplund, N. Grevesse, and A. J. Sauval. The Solar Ghemical Gomposition. 

1550 In T. G. Barnes HI & F. N. Bash, editor. Cosmic Abundances as Records of 

1551 Stellar Evolution and Nucleosynthesis, volume 336 of Astronomical Society 

1552 of the Pacific Conference Series, page 25, September 2005. 

1553 J. N. Bahcall, M. H. Pinsonneault, and S. Basu. Solar Models: Current Epoch 

1554 and Time Dependences, Neutrinos, and Helioseismological Properties. As- 

1555 trophysical Journal, 555:990-1012, 2001. doi: 10.1086/321493. 

1556 J. L. Bandficld, A. Deanne Rogers, and C. S. Edwards. The role of aqueous 

1557 alteration in the formation of martian soils. Icarus, 211:157-171, 2011. doi: 

1558 10.1016/j.icarus.2010.08.028. 

1559 S. Barabash, A. Fedorov, R. Lundin, and J. -A. Sauvaud. Martian Atmospheric 

1560 Erosion Rates. Science, 315:501-, 2007. doi: 10.1126/science.ll34358. 

1561 Y. I. Baranov, W. J. Lafferty, and G. T. Eraser. Infrared spectrum of the 

1562 continuum and dimer absorption in the vicinity of the O2 vibrational fun- 

1563 damental in O2/CO2 mixtures. Journal of Molecular Spectroscopy, 228: 

1564 432-440, December 2004. doi: 10.1016/j.jms.2004.04.010. 

1565 G.J. Barnhart and F. Nimmo. Role of impact excavation in distributing clays 

1566 over Noachian surfaces. Journal of Geophysical Research (Planets), 116: 

1567 E01009, January 2011. 

1568 G. Bergcr, M. J. Toplis, E. Trcguicr, C. D'Uston, and P. Pinct. Evidence in 

1569 favor of small amounts of ephemeral and transient water during alteration 

1570 at Meridiani Planum, Mars. American Mineralogist, 94:1279-1282, 2009. 

1571 J.-P. Bibring, R. E. Arvidson, A. Gendrin, B. Gondet, Y. Langevin, S. Le 

1572 Mouelic, N. Mangold, R. V. Morris, J. F. Mustard, F. Poulet, C. Quantin, 

1573 and C. Sotin. Coupled Ferric Oxides and Sulfates on the Martian Surface. 

1574 Science, 317:1206-, 2007. doi: 10.1126/science.ll44174. 

1575 J. P. Bibring, Y. Langevin, J. F. Mustard, F. Poulet, R. Arvidson, A. Gendrin, 

1576 B. Gondet, N. Mangold, P. Pinet, and F. Forget. Global Mineralogical and 

1577 Aqueous Mars History Derived from OMEGA/Mars Express Data. Science, 

1578 312:400-404, 2006. doi: 10.1126/sciencc.ll22659. 

1579 J. L. Bishop, E. Z. N. Dobrea, N. K. McKeown, M. Parente, B. L. Ehlmann, 

1580 J. R. Michalski, R. E. MiUiken, F. Poulet, G. A. Swayze, J. F. Mustard, 

1581 S. L. Murchie, and J.-P. Bibring. Phyllosilicate Diversity and Past Aqueous 

1582 Activity Revealed at Mawrth Vallis, Mars. Science, 321:830-, August 2008. 

1583 doi: 10.1126/science.ll59699. 

1584 C. F. Bohren and D. R. Huffman. Absorption and scattering of light by small 

1585 particles. 1983. 



61 



1586 W. V. Boynton, W. C. Feldman, S. W. Squyres, T. H. Prettyman, J. Bruckner, 

1587 L. G. Evans, R. C. Reedy, R. Starr, J. R. Arnold, D. M. Drake, R A. J. 

1588 Englert, A. E. Metzger, I. Mitrofanov, J. I. Trombka, C. d'Uston, H. Wanke, 

1589 O. Gasnault, D. K. Hamara, D. M. Janes, R. L. Marcialis, S. Maurice, 

1590 I. Mikheeva, G. J. Taylor, R. Tokar, and G. Shinohara. Distribution of 

1591 Hydrogen in the Near Surface of Mars: Evidence for Subsurface Ice Deposits. 

1592 Science, 297:81-85, 2002. doi: 10. 1126/science. 1073722. 

1593 W. V. Boynton, D. W. Ming, S. P. Kounaves, S. M. M. Young, R. E. Arvidson, 

1594 M. H. Hecht, J. Hoffman, P. B. Niles, D. K. Hamara, R. G. Quinn, P. H. 

1595 Smith, B. Sutter, D. G. Gatling, and R. V. Morris. Evidence for Galcium 

1596 Garbonate at the Mars Phoenix Landing Site. Science, 325:61-, July 2009. 

1597 doi: 10.1126/science.ll72768. 

1598 B. A. Bradley, S. E. H. Sakimoto, H. Frey, and J. R. Zimbelman. Medusae 

1599 Fossae Formation: New perspectives from Mars Global Surveyor. Journal 

1600 of Geophysical Research (Planets), 107:5058, August 2002. doi: 10.1029/ 

1601 2001JE001537. 

1602 R. E. Brandt and S. G. Warren. Solar-heating rates and temperature profiles 

1603 in Antarctic snow and ice. Journal of Glaciology, 39:99-110, 1993. 

1604 N. T. Bridges, M. E. Banks, R. A. Beyer, F. G. Ghuang, E. Z. Noe Dobrea, 

1605 K. E. Herkenhoff, L. P. Keszthelyi, K. E. Fishbaugh, A. S. McEwen, T. I. 

1606 Michaels, B. J. Thomson, and J. J. Wray. Aeolian bedforms, yardangs, and 

1607 indurated surfaces in the Tharsis Montes as seen by the HiRISE Camera: 

1608 Evidence for dust aggregates. Icarus, 205:165-182, 2010. doi: 10.1016/j. 

1609 icarus.2009.05.017. 

1610 N.T. Bridges, M.G. Bourke, P.E. Geissler, M.E. Banks, G. Golon, S. Din- 

1611 iega, M.P. Golombek, G.J. Hansen, S. Mattson, A.S. McEwen, M.T. Mel- 

1612 Ion, N. Stantzos, and B.J. Thomson. Planet- wide sand motion on mars. 

1613 Geology, 40(l):31-34, 2012. doi: 10.1130/G32373.1. 

1614 B. W. Brock, I. G. Willis, and M. J. Sharp. Measurement and param- 

1615 eterization of aerodynamic roughness length variations at Haut Glacier 

1616 d'Arolla, Switzerland. Journal of Glaciology, 52:281-297, 2006. doi: 

1617 10.3189/172756506781828746. 

1618 M. A. Bullock and J. M. Moore. Atmospheric conditions on early Mars and 

1619 the missing layered carbonates. Geophysical Research Letters, 34:L19201, 

1620 October 2007. doi: 10.1029/2007GL030688. 

1621 D. M. Burr, M.-T. Enga, R. M. E. WiUiams, J. R. Zimbelman, A. D. 

1622 Howard, and T. A. Brennand. Pervasive aqueous paleoflow features in 

1623 the Aeolis/Zephyria Plana region. Mars. Icarus, 200:52-76, 2009. doi: 

1624 10. 1016/j.icarus.2008. 10.014. 

1625 D. M. Burr, R. M. E. WiUiams, K. D. Wendell, M. Ghojnacki, and J. P. Emery. 

1626 Inverted fluvial features in the Aeolis/Zephyria Plana region. Mars: Forma- 

1627 tion mechanism and initial paleodischarge estimates. Journal of Geophysical 

1628 Research (Planets), 115:E07011, 2010. doi: 10.1029/2009JE003496. 

1629 S. B. Gadieux. Gonstraining Martian sedimentation via analysis of 

1630 stratal packaging, intracrater layered deposits, Arabia Terra, Mars, 



62 



1631 http://trace.tennessee.edu/utk_gradthes/860 . Master's thesis, Uni- 

1632 versity of Tennessee at Knoxville, 2011. 

1633 S. B. Cadieux and L. C. Kah. Intracrater Layered Deposits in Arabia Terra, 

1634 Mars Indicate Potential Wet, Cold, Conditions in Late Noachian-Early Hes- 

1635 perian. In Lunar and Planetary Institute Science Conference Abstracts, vol- 

1636 ume 42 of Lunar and Planetary Institute Science Conference Abstracts, page 

1637 1265, March 2011. 

1638 W. M. Calvin, L. H. Roach, F. P. Seelos, K. D. Seelos, R. O. Green, S. L. 

1639 Murchie, and J. F. Mustard. Compact Reconnaissance Imaging Spectrom- 

1640 eter for Mars observations of northern Martian latitudes in summer. Jour- 

1641 nal of Geophysical Research (Planets), 114:E00D11, December 2009. doi: 

1642 10. 1029/2009 JE003348. 

1643 M. H. Carr and J. W. Head. Basal melting of snow on early Mars: A possible 

1644 origin of some valley networks. Geophysical Research Letters, 30:2245, 2003. 

1645 doi: 10.1029/2003GL018575. 

1646 M. H. Carr and J. W. Head. Geologic history of Mars. Earth and Planetary 

1647 Science Letters, 294:185-203, 2010. doi: 10.1016/j.epsl.2009.06.042. 

1648 D. C. Cathng, S. E. Wood, C. Leovy, D. R. Montgomery, H. M. Greenberg, 

1649 C. R. Glein, and J. M. Moore. Light-toned layered deposits in Juventae 

1650 Chasma, Mars. Icarus, 181:26-51, 2006. doi: 10.1016/j.icarus.2005.10.020. 

1651 J. D. Chittenden, V. Chevrier, L. A. Roe, K. Bryson, R. Pilgrim, and D. W. G. 

1652 Sears. Experimental study of the effect of wind on the stability of water ice 

1653 on Mars. Icarus, 196:477-487, 2008. doi: 10.1016/j.icarus.2008.01.016. 

1654 P. R. Christensen. Formation of recent martian gullies through melting of 

1655 extensive water-rich snow deposits. Nature, 422:45-48, March 2003. 

1656 G. D. Clow. Generation of liquid water on Mars through the melting of a dusty 

1657 snowpack. Icarus, 72:95-127, 1987. doi: 10.1016/0019-1035(87)90123-0. 

1658 A. Colaprete and O. B. Toon. Carbon dioxide clouds in an early dense Martian 

1659 atmosphere. Journal of Geophysical Research (Planets), 108:5025, 2003. doi: 

1660 10.1029/2002JE001967. 

1661 A. Colaprete, J. R. Barnes, R. M. Haberle, J. L. HoUingsworth, H. H. Kieffer, 

1662 and T. N. Titus. Albedo of the south pole on Mars determined by topo- 

1663 graphic forcing of atmosphere dynamics. Nature, 435:184-188, May 2005. 

1664 doi: 10.1038/nature03561. 

1665 F. Costard, F. Forget, N. Mangold, and J. P. Peulvast. Formation of Recent 

1666 Martian Debris Flows by Melting of Near-Surface Ground Ice at High Obliq- 

1667 uity. Science, 295:110-113, January 2002. doi: 10.1126/science.295.5552.110. 

1668 K. Cuffey and W. S. B. Paterson. The Physics of Glaciers, 4th edition: Aca- 

1669 demic Press. 2010. 

1670 S. Cull, R. E. Arvidson, M. Mellon, S. Wiseman, R. Clark, T. Titus, R. V. 

1671 Morris, and P. McGuire. Seasonal H2O and CO2 ice cycles at the Mars 

1672 Phoenix landing site: 1. Prelanding CRISM and HiRlSE observations. 

1673 Journal of Geophysical Research (Planets), 115:E00D16, April 2010. doi: 

1674 10.1029/2009JE003340. 

1675 P. T. Doran, R. A. Wharton, D. J. Des Marais, and C. P. McKay. Antarc- 



63 



1676 tic palcolake sediments and the search for extinct life on Mars. Journal 

1677 of Geophysical Research, 103:28481-28494, November 1998. doi: 10.1029/ 

1678 98JE01713. 

1679 Peter T. Doran, W. Berry Lyons, and Diane M. McKnight. Life in Antarctic 

1680 deserts and other cold dry environments: astrobiological analogs. Cambridge 

1681 University Press, 2010. 

1682 L. Drube, K. Leer, W. Goetz, H. P. Gunnlaugsson, M. P. Haspang, N. Laurit- 

1683 sen, M. B. Madsen, L. K. D. S0rensen, M. D. EUehoj, M. T. Lemmon, R. V. 

1684 Morris, D. Blaney, R. O. Reynolds, and P. H. Smith. Magnetic and optical 

1685 properties of airborne dust and settling rates of dust at the Phoenix landing 

1686 site. Journal of Geophysical Research (Planets), 115:E00E23, August 2010. 

1687 doi: 10. 1029/2009 JE003419. 

1688 C. M. Dundas and S. Byrne. Modeling sublimation of ice exposed by new 

1689 impacts in the Martian mid-latitudes. Icarus, 206:716-728, 2010. doi: 10. 

1690 1016/j.icarus.2009.09.007. 

1691 S. D. Eckermann, J. Ma, and X. Zhu. Scale-dependent infrared radiative 

1692 damping rates on Mars and their role in the deposition of gravity-wave 

1693 momentum flux. Icarus, 211:429-442, January 2011. doi: 10.1016/j.icarus. 

1694 2010.10.029. 

1695 K. S. Edgett. The sedimentary rocks of Sinus Meridiani: Five key observations 

1696 from data acquired by the Mars Global Surveyor and Mars Odyssey orbitcrs. 

1697 International Journal of Mars Science and Exploration, 1:5-58, 2005. doi: 

1698 10.1555/mars.2005.0002. 

1699 K. S. Edgett and M. G. Malin. Martian sedimentary rock stratigraphy: Out- 

1700 crops and interbedded craters of northwest Sinus Meridiani and southwest 

1701 Arabia Terra. Geophysical Research Letters, 29(24) :2179, December 2002. 

1702 doi: 10.1029/2002GL016515. 

1703 G. S. Edwards, J. L. Bandficld, P. R. Ghristensen, and R. L. Fergason. Global 

1704 distribution of bedrock exposures on Mars using THEMIS high-resolution 

1705 thermal inertia. Journal of Geophysical Research (Planets), 114:E11001, 

1706 November 2009. doi: 10.1029/2009JE003363. 

1707 B. L. Ehlmann, J. F. Mustard, S. L. Murchie, J. -P. Bibring, A. Meunier, A. A. 

1708 Fraeman, and Y. Langevin. Subsurface water and clay mineral formation 

1709 during the early history of Mars. Nature, 479:53-60, November 2011. doi: 

1710 10.1038/naturel0582. 

1711 M.E. Elwood Madden, A.S. Madden, and J.D Rimstidt. How long was Merid- 

1712 iani Planum wet? Applying a jarosite stopwatch to constrain the duration 

1713 of diagenesis. Geology, 37:635, 2009. 

1714 A. G. Fairen, A. F. Davila, L. Gago-Duport, R. Amils, and G. P. McKay. 

1715 Stability against freezing of aqueous solutions on early Mars. Nature, 459: 

1716 401-404, May 2009. doi: 10.1038/nature07978. 

1717 C. I. Fassett and J. W. Head. The timing of martian valley network activity: 

1718 Gonstraints from buffered crater counting. Icarus, 195:61-89, 2008. doi: 

1719 10.1016/j.icarus.2007.12.009. 

1720 C. L Fassett and J. W. Head. Sequence and timing of conditions on early 



64 



1721 Mars. Icarus, 211:1204-1214, 2011. doi: 10.1016/j.icarus.2010.11.014. 

1722 C. I. Fassett, J. L. Dickson, J. W. Head, J. S. Levy, and D. R. Marchant. 

1723 Supraglacial and proglacial valleys on Amazonian Mars. Icarus, 208:86- 

1724 100, 2010. doi: 10.1016/j.icarus.2010.02.021. 

1725 J. L. Fastook, J. W. Head, D. R. Marchant, and F. Forget. Tropical 

1726 mountain glaciers on Mars: AltitTide-dependence of ice accumulation, ac- 

1727 cumulation conditions, formation times, glacier dynamics, and implications 

1728 for planetary spin-axis/orbital history. Icarus, 198:305-317, 2008. doi: 

1729 10.1016/j.icarus.2008.08.008. 

1730 L. K. Fenton and R. K. Hayward. Southern high latitude dune fields on Mars: 

1731 Morphology, aeolian inactivity, and climate change. Geomorphology, 121: 

1732 98-121, 2010. doi: 10.1016/j.geomorph.2009.11.006. 

1733 F. Forget, R. M. Haberle, F. Montmessin, B. Lcvrard, and J. W. Head. For- 

1734 mation of glaciers on Mars by atmospheric precipitation at high obliquity. 

1735 Science, 311:368-371, 2006. doi: 10. 1126/science. 1120335. 

1736 N. K. Forsberg- Taylor, A. D. Howard, and R. A. Craddock. Crater degrada- 

1737 tion in the Martian highlands: Morphometric analysis of the Sinus Sabaeus 

1738 region and simulation modeling suggest fluvial processes. Journal of Geo- 

1739 physical Research (Planets), 109:E05002, 2004. doi: 10.1029/2004JE002242. 

1740 P. E. Geissler. Three decades of Martian surface changes. Journal of Geo- 

1741 physical Research (Planets), 110:E02001, 2005. doi: 10.1029/2004JE002345. 

1742 P. E. Geissler, R. Sullivan, M. Golombek, J. R. Johnson, K. Herkenhoff, 

1743 N. Bridges, A. Vaughan, J. Maki, T. Parker, and J. Bell. Gone with the 

1744 wind: Eolian erasure of the Mars Rover tracks. Journal of Geophysical Re- 

1745 search (Planets), 115:E00F11, 2010. doi: 10.1029/2010JE003674. 

1746 T. D. Glotch and P. R. Christensen. Geologic and mineralogic mapping of 

1747 Aram Chaos: Evidence for a water-rich history. Journal of Geophysical Re- 

1748 search (Planets), 110:E09006, September 2005. doi: 10.1029/2004JE002389. 

1749 M. Golombek, K. Robinson, A. McEwen, N. Bridges, B. Ivanov, L. Tornabene, 

1750 and R. Sullivan. Constraints on ripple migration at Meridiani Planum from 

1751 Opportunity and HiRISE observations of fresh craters. Journal of Geophys- 

1752 ical Research (Planets), 115:E00F08, 2010. doi: 10.1029/2010JE003628. 

1753 M. P. Golombek, J. A. Grant, L. S. Grumpier, R. Greeley, R. E. Arvidson, 

1754 J. F. Bell, C. M. Weitz, R. Sullivan, P. R. Christensen, L. A. Soderblom, and 

1755 S. W. Squyres. Erosion rates at the Mars Exploration Rover landing sites 

1756 and long-term climate change on Mars. Journal of Geophysical Research 

1757 (Planets), 111:E12S10, 2006. doi: 10.1029/2006JE002754. 

1758 R. Goody and M. J. S. Belton. Radiative relaxation times for Mars. A discus- 

1759 sion of martian atmospheric dynamics. Planetary & Space Science, 15:247, 

1760 February 1967. doi: 10.1016/0032-0633(67)90193-6. 

1761 J. A. Grant and T. J. Parker. Drainage evolution in the Margaritifer Sinus re- 

1762 gion. Mars. Journal of Geophysical Research (Planets), 107:5066, September 

1763 2002. doi: 10.1029/2001JE001678. 

1764 J. A. Grant and S. A. Wilson. Late alluvial fan formation in southern Mar- 

1765 garitifer Terra, Mars. Geophysical Research Letters, 380:L08201, 2011. doi: 



65 



1766 10. 1029/201 1GL046844. 

1767 J. A. Grant, R. P. Irwin, III, J. P. Grotzinger, R. E. Milliken, L. L. Tornabene, 

1768 A. S. McEwen, C. M. Wcitz, S. W. Squyres, T. D. Glotch, and B. J. Thom- 

1769 son. HiRISE imaging of impact mcgabreccia and sub-meter aqueous strata 

1770 in Holden Crater, Mars. Geology, 36:195-198, March 2008. 

1771 J. A. Grant, S. A. Wilson, E. Noc Dobrea, R. L. Fergason, J. L. Griffes, 

1772 J. M. Moore, and A. D. Howard. HiRISE views enigmatic deposits in the 

1773 Sirenum Fossae region of Mars. Icarus, 205:53-63, January 2010. doi: 10. 

1774 1016/j.icarus.2009.04.009. 

1775 S. E. Grasby, C. C. Allen, T. G. Longazo, J. T. Lisle, D. W. Griffin, and 

1776 B. Beauchamp. Supraglacial Sulfur Springs and Associated Biological Ac- 

1777 tivity in the Canadian High Arctic-Signs of Life Beneath the Ice. Astrobi- 

1778 ology, 3:583-596, November 2003. doi: 10.1089/153110703322610672. 

1779 R. Greeley and J. Guest. I-1802B: Geological map of the eastern equatorial 

1780 regions of Mars. Technical report, U.S. Geological Survey, 1987. 

1781 J. Grotzinger, J. Bell, HI, K. Herkenhoff, J. Johnson, A. Knoll, E. McCartney, 

1782 S. McLennan, J. Metz, J. Moore, S. Squyres, R. Sullivan, O. Ahronson, 

1783 R. Arvidson, B. Joliff, M. Golombek, K. Lewis, T. Parker, and J. Soderblom. 

1784 Sedimentary textures formed by aqueous processes, Erebus crater, Meridiani 

1785 Planum, Mars. Geology, 34:1085, 2006. doi: 10.1130/G22985A.1. 

1786 J. P. Grotzinger and R. E. Milliken. The sedimentary rock record of Mars: 

1787 Distribution, origins, and global stratigraphy. In Grotzinger, J. P., editor, 

1788 Sedimentary Geology of Mars, volume 11 of Special Publications, page 1. 

1789 SEPM (Society for Sedimentary Geology), 2012. 

1790 M. Gruszka and A. Borysow. Roto-Translational Colhsion-Induced Absorption 

1791 of C0_2 for the Atmosphere of Venus at Frequencies from to 250 cm"-l, 

1792 at Temperatures from 200 to 800 K. Icarus, 129:172-177, September 1997. 

1793 doi: 10. 1006/icar. 1997.5773. 

1794 M. Gruszka and A. Borysow. Computer simulation of the far infrared collision 

1795 induced absorption spectra of gaseous C02. Molecular Physics, 93:1007- 

1796 1016, 1998. doi: 10.1080/00268979809482287. 

1797 J. A. Guzik and K. Mussack. Exploring Mass Loss, Low-Z Accretion, and 

1798 Convective Overshoot in Solar Models to Mitigate the Solar Abundance 

1799 Problem. Astrophyszcal Journal, 713:1108-1119, April 2010. doi: 10.1088/ 

1800 0004-637X/713/2/1108. 

1801 R. M. Haberle. Early Mars Climate Models. Journal of Geophysical Research, 

1802 1032:28467-28480, 1998. doi: 10.1029/98JE01396. 

1803 R. M. Haberle, J. B. Pollack, J. R. Barnes, R. W. Zurck, C. B. Leovy, J. R. 

1804 Murphy, H. Lee, and J. Schaeffer. Mars atmospheric dynamics as simu- 

1805 lated by the NASA Ames General Circulation Model. I - The zonal-mean 

1806 circulation. Journal of Geophysical Research, 98:3093-3123, 1993. doi: 

1807 10. 1029/92 JE02946. 

1808 R. M. Haberle, J. R. Murphy, and J. Schaeffer. Orbital change experiments 

1809 with a Mars general circulation model. Icarus, 161:66-89, 2003. doi: 10. 

1810 1016/80019-1035(02)00017-9. 



66 



1811 I. Halevy and J. W. Head. Punctuated Volcanism, Transient Warming and 

1812 Global Change in the Late Noachian-Early Hesperian. In Lunar and Plane- 

1813 tary Institute Science Conference Abstracts, volume 43 of Lunar and Pian- 
ist etary Institute Science Conference Abstracts, page 1908, March 2012. 

1815 I. Halevy and D. P. Schrag. Sulfur dioxide inhibits calcium carbonate precipi- 

1816 tation: Implications for early Mars and Earth. Ceophysical Research Letters, 
1317 36:123201, December 2009. doi: 10.1029/2009GL040792. 

1818 I. Halevy, M. T. Zuber, and D. P. Schrag. A Sulfur Dioxide Climate Feedback 

1819 on Early Mars. Science, 318:1903, 2007. doi: 10.1126/science.ll47039. 

1820 I. Halevy, R. T. Pierrehumbert, and D. P. Schrag. Radiative transfer in CO2- 

1821 rich palcoatmosphcres. Journal of Geophysical Research (Atmospheres), 
1322 114:D18112, 2009. doi: 10.1029/2009JD011915. 

1823 I. Halevy, W. W. Fischer, and J. M. Filer. Carbonates in the Martian meteorite 

1824 Allan Hills 84001 formed at 18 4 C in a near-surface aqueous environment. 

1825 Proceedings of the National Academy of Science, 108:16895-16899, October 

1826 2011. doi: 10.1073/pnas.ll09444108. 

1827 B. L. Hall, G. H. Denton, A. G. Fountain, C. H. Hendy, and G. M. Hender- 

1828 son. Antarctic lakes suggest millennial reorganizations of Southern Hemi- 

1829 sphere atmospheric and oceanic circulation. Proceedings of the National 

1830 Academy of Sciences, 107:21355-21359, December 2010. doi: 10.1073/pnas. 

1831 1007250107. 

1832 E. Hauber, S. van Gasselt, M. G. Chapman, and G. Neukum. Geomorphic 
1333 evidence for former lobate debris aprons at low latitudes on Mars: Indicators 

1834 of the Martian paleoclimate. Journal of Geophysical Research (Planets), 

1835 113:E02007, 2008. doi: 10.1029/2007JE002897. 

1836 E. M. Hausrath, A. K. Navarre-Sitchler, P. B. Sak, C. I. Stccfcl, and S. L. 

1837 Brantley. Basalt weathering rates on Earth and the duration of liquid water 

1338 on the plains of Guscv Crater, Mars. Geology, 36:67-70, 2008. 

1339 J. W. Head. Geologic Evidence for Latitude-Dependent Water-Related De- 

1840 posits on Mars: Implications for Climate History and the Hydrological Cycle 

1841 on Mars. In F. Forget & E. Millour, editor. Mars Atmosphere: Modelling 

1842 and observation, pages 427-430, February 2011. 

1343 M. H. Hccht. Metastability of Liquid Water on Mars. Icarus, 156:373-386, 

1844 2002. doi: 10.1006/icar.2001.6794. 

1845 J. L. Heldmann. Formation and evolution of buried snowpack deposits in 

1846 Pearse Valley, Antarctica, and implications for Mars. Antarctic Science, 

1847 page TBD, 2012. 

1348 J. L. Heldmann and M. T. Mellon. Observations of martian gullies and con- 

1849 straints on potential formation mechanisms. Icarus, 168:285-304, April 

1850 2004. doi: 10.1016/j.icarus.2003.11.024. 

1851 C.H. Hendy. Late Quaternary lakes in the McMurdo Sound region of Antarc- 

1852 tica. Geografiska Annaler Series A - Physical Gegraphy, 82A(2-3):411-432, 

1353 2000. 

1854 K. E. Herkenhoff, M. P. Golombek, E. A. Guinness, J. B. Johnson, A. Kusack, 

1855 L. Richter, R. J. Sullivan, and S. Gorevan. In situ observations of the 



67 



1856 physical properties of the Martian surface. In Bell, J., Ill, editor, The 

1857 Martian Surface - Composition, Mineralogy, and Physical Properties, page 

1858 451. 2008. 

1859 T. L. Hudson and O. Aharonson. Diffusion barriers at Mars surface conditions: 

1860 Salt crusts, particle size mixtures, and dust. Journal of Geophysical Research 

1861 (Planets), 113:E09008, 2008. doi: 10.1029/2007JE003026. 

1862 T. L. Hudson, O. Aharonson, and N. Schorghofer. Laboratory experiments 

1863 and models of diffusive emplacement of ground ice on Mars. Journal of Geo- 

1864 physical Research (Planets), 114:E01002, 2009. doi: 10.1029/2008JE003149. 

1865 J. A. Hurowitz and S. M. McLennan. A ~3.5 Ga record of water-hmited, 

1866 acidic weathering conditions on Mars. Earth and Planetary Science Letters, 

1867 260:432-443, 2007. doi: 10.1016/j.epsl.2007.05.043. 

1868 J. A. Hurowitz, W. W. Fischer, N. J. Tosca, and R. E. Milhken. Origin of 

1869 acidic surface waters and the evolution of atmospheric chemistry on early 

1870 Mars. Nature Geoscience, 3:323-326, 2010. doi: 10.1038/ngeo831. 

1871 B. M. Hynek and R. J. Phillips. The stratigraphy of Meridiani Planum, Mars, 

1872 and implications for the layered deposits' origin. Earth and Planetary Sci- 

1873 ence Letters, 274:214-220, September 2008. doi: 10.1016/j.epsl.2008.07.025. 

1874 B. M. Hynek, R. J. Phillips, and R. E. Arvidson. Explosive volcanism in the 

1875 Tharsis region: Global evidence in the Martian geologic record. Journal of 

1876 Geophysical Research (Planets), 108:5111, September 2003. doi: 10.1029/ 

1877 2003JE002062. 

1878 B. M. Hynek, M. Beach, and M. R. T. Hoke. Updated global map of 

1879 Martian valley networks and implications for climate and hydrologic pro- 

1880 cesses. Journal of Geophysical Research (Planets), 115:E09008, 2010. doi: 

1881 10.1029/2009JE003548. 

1882 R. P. Irwin, HI, R. A. Craddock, and A. D. Howard. Interior channels in 

1883 Martian valley networks: Discharge and runoff production. Geology, 33:489, 

1884 2005. doi: 10.1130/G21333.1. 

1885 B. M. Jakosky and M. H. Carr. Possible precipitation of ice at low latitudes 

1886 of Mars during periods of high obliquity. Nature, 315:559-561, 1985. doi: 

1887 10.1038/315559a0. 

1888 R. Jaumann, A. Nass, D. Tirsch, D. Reiss, and G. Neukum. The Western 

1889 Libya Montes Valley System on Mars: Evidence for episodic and multi- 

1890 genetic erosion events during the Martian history. Earth and Planetary 

1891 Science Letters, 294:272-290, June 2010. doi: 10.1016/j.epsl.2009.09.026. 

1892 D. J. Jerolmack and C. Paola. Shredding of environmental signals by sediment 

1893 transport. Geophysical Research Letters, 37:L19401, October 2010. doi: 

1894 10.1029/2010GL044638. 

1895 J. R. Johnson, W. M. Grundy, and M. T. Lemmon. Dust deposi- 

1896 tion at the Mars Pathfinder landing site: observations and modeling of 

1897 visible/near-infrared spectra. Icarus, 163:330-346, June 2003. doi: 10.1016/ 

1898 S0019-1035(03)00084-8. 

1899 S. S. Johnson, M. A. Mischna, T. L. Grove, and M. T. Zuber. Sulfur-induced 

1900 greenhouse warming on early Mars. Journal of Geophysical Research (Plan- 



68 



1901 ets), 113:E08005, 2008. doi: 10.1029/2007JE002962. 

1902 S. J. Kadish, N. G. Barlow, and J. W. Head. Latitude dependence of Martian 

1903 pedestal craters: Evidence for a sublimation-driven formation mechanism. 

1904 Journal of Geophysical Research (Planets), 114:E10001, October 2009. doi: 

1905 10.1029/2008JE003318. 

1906 S. J. Kadish, J. W. Head, and N. G. Barlow. Pedestal crater heights on Mars: 

1907 A proxy for the thicknesses of past, ice-rich, Amazonian deposits. Icarus, 

1908 210:92-101, 2010. doi: 10.1016/j.icarus.2010.06.021. 

1909 R. Kahn. The evolution of CO2 on Mars. Icarus, 62:175-190, 1985. doi: 

1910 10.1016/0019-1035(85)90116-2. 

1911 M. A. Kahre, J. R. Murphy, and R. M. Haberle. Modeling the Martian dust 

1912 cycle and surface dust reservoirs with the NASA Ames general circulation 

1913 model. Journal of Geophysical Research (Planets), 111:E06008, 2006. doi: 

1914 10.1029/2005JE002588. 

1915 M. A. Kahre, S. Vines, R. M. Haberle, J. HoUingsworth, and C.B.. Leovy. 

1916 Stabilizing the Early Martian Glimate: Effects of Airborne Dust, GO2 Ice 

1917 Gap Albedo, and Orbital Obliquity on Atmospheric GoUapse. AGU Fall 

1918 Meeting Abstracts, 2011. 

1919 J. F. Kasting. GO2 condensation and the climate of early Mars. Icarus, 94: 

1920 1-13, 1991. doi: 10.1016/0019-1035(91)90137-1. 

1921 L. Kerber and J. W. Head. The age of the Medusae Fossae Formation: 

1922 Evidence of Hesperian emplacement from crater morphology, stratigra- 

1923 phy, and ancient lava contacts. Icarus, 206:669-684, April 2010. doi: 

1924 10. 1016/j.icarus.2009. 10.001. 

1925 A. Kereszturi, M. Vincendon, and F. Schmidt. Water ice in the dark dune 

1926 spots of Richardson crater on Mars. Planetary and Space Science, 59:26-42, 

1927 2011. doi: 10.1016/j.pss.2010.10.015. 

1928 H. H. Kicffcr. H2O grain size and the amount of dust in Mars' residual north 

1929 polar cap. Journal of Geophysical Research, 95:1481-1493, February 1990. 

1930 doi: 10.1029/JB095iB02p01481. 

1931 K. M. Kinch, J. Sohl-Dickstein, J. F. Bell, J. R. Johnson, W. Goctz, and G. A. 

1932 Landis. Dust deposition on the Mars Exploration Rover Panoramic Gamera 

1933 (Pancam) calibration targets. Journal of Geophysical Research (Planets), 

1934 112:E06S03, April 2007. doi: 10.1029/2006JE002807. 

1935 E. S. Kite, I. Matsuyama, M. Manga, J. T. Perron, and J. X. Mitrovica. 

1936 True Polar Wander driven by latc-stage volcanism and the distribution of 

1937 paleopolar deposits on Mars. Earth and Planetary Science Letters, 280: 

1938 254-267, April 2009. doi: 10.1016/j.epsl.2009.01.040. 

1939 E. S. Kite, E. Gaidos, and M. Manga. Glimate Instability on Tidally Locked 

1940 Exoplanets. Astrophysical Journal, 743:41, December 2011a. doi: 10.1088/ 

1941 0004-637X/743/1/41. 

1942 E. S. Kite, T. I. Michaels, S. Rafkin, M. Manga, and W. E. Dietrich. Lo- 

1943 calized precipitation and runoff on Mars. Journal of Geophysical Research 

1944 (Planets), 116:E07002, July 2011b. doi: 10.1029/2010JE003783. 

1945 E. S. Kite, S. Rafkin, T. I. Michaels, W. E. Dietrich, and M. Manga. Ghaos 



69 



1946 terrain, storms, and past climate on Mars. Journal of Geophysical Research 

1947 (Planets), 116:E10002, October 2011c. doi: 10.1029/2010JE003792. 

1948 A. H. Knoll, B. L. JoUiff, W. H. Farrand, J. F. Bell, III, B. C. Clark, R. Gellert, 

1949 M. P. Golombek, J. P. Grotzinger, K. E. Herkenhoff, J. R. Johnson, S. M. 

1950 McLennan, R. Morris, S. W. Squyres, R. Sullivan, N. J. Tosca, A. Yen, 

1951 and Z. Learner. Vcnccrs, rinds, and fracture fills: Relatively late alteration 

1952 of sedimentary rocks at Mcridiani Planum, Mars. Journal of Geophysical 

1953 Research (Planets), 113:E06S16, 2008. doi: 10.1029/2007JE002949. 

1954 G. Kopp and J. L. Lean. A new, lower value of total solar irradiance: Evidence 

1955 and climate significance. Geophysical Research Letters, 38:L01706, January 

1956 2011. doi: 10.1029/2010GL045777. 

1957 E. R. Kraal, M. van Dijk, G. Postma, and M. G. Kleinhans. Martian stepped- 

1958 delta formation by rapid water release. Nature, 451:973-976, 2008. doi: 

1959 10.1038/Nature06615. 

1960 M. A. Kreslavsky and J. W. Head. Kilometer-scale roughness of Mars: Results 

1961 from MOLA data analysis. Journal of Geophysical Research, 105:26695- 

1962 26712, November 2000. doi: 10. 1029/2000 JE001259. 

1963 M. A. Kreslavsky and J. W. Head. North-south topographic slope asymmetry 

1964 on Mars: Evidence for insolation-related erosion at high obliquity. Geophys- 

1965 ical Research Letters, 30(15):1815, 2003. doi: 10.1029/2003GL017795. 

1966 J. Laskar. Chaotic diffusion in the Solar System. Icarus, 196:1-15, 2008. doi: 

1967 10.1016/j.icarus.2008.02.017. 

1968 J. Laskar and P. Robutel. The chaotic obliquity of the planets. Nature, 361: 

1969 608-612, 1993. doi: 10.1038/361608a0. 

1970 J. Laskar, A. C. M. Correia, M. Gastineau, F. Joutel, B. Levrard, and P. Robu- 

1971 tel. Long term evolution and chaotic diffusion of the insolation quantities 

1972 of Mars. Icarus, 170:343-364, 2004. doi: 10.1016/j.icarus.2004.04.005. 

1973 P. Lee and C. P. McKay. Mars: Always Cold, Sometimes Wet? In S. Mackwell 

1974 & E. Stansbery, editor. Lunar and Planetary Institute Science Gonference 

1975 Abstracts, volume 34 of Lunar and Planetary Institute Science Gonference 

1976 Abstracts, page 2127, March 2003. 

1977 A. Lefort, P. S. Russell, N. Thomas, A. S. McEwen, C. M. Dundas, and R. L. 

1978 Kirk. Observations of periglacial landforms in Utopia Planitia with the High 

1979 Resolution Imaging Science Experiment (HiRISE). Journal of Geophysical 

1980 Research (Planets), 114:E04005, April 2009. doi: 10.1029/2008JE003264. 

1981 A. Lefort, P. S. Russell, and N. Thomas. Scalloped terrains in the Peneus and 

1982 Amphitrites Paterae region of Mars as observed by HiRISE. Icarus, 205: 

1983 259-268, January 2010. doi: 10.1016/j.icarus.2009.06.005. 

1984 B. Levrard, F. Forget, F. Montmessin, and J. Laskar. Recent ice-rich deposits 

1985 formed at high latitudes on Mars by sublimation of unstable equatorial ice 

1986 during low obhquity. Nature, 431:1072-1075, October 2004. doi: 10.1038/ 

1987 nature03055. 

1988 K. J. Lewis. Solar-forced roughening of Antarctic glaciers and the Martian 

1989 icecaps: How surficial debris and roughness affect glacial melting in Taylor 

1990 Valley, Antarctica and how this can be applied to the Martian icecaps. PhD 



70 



1991 thesis, University of Colorado at Boulder, 2001. 

1992 K. W. Lewis. The rock record of Mars : structure, sedimentology and stratig- 

1993 raphy. PhD thesis, Caltech, 2009. 

1994 K. W. Lewis, O. Aharonson, J. P. Grotzinger, R. L. Kirk, A. S. McEwen, and 

1995 T.-A. Suer. Quasi-Periodic Bedding in the Sedimentary Rock Record of 

1996 Mars. Science, 322:1532-, 2008. doi: 10.1126/science.ll61870. 

1997 K. W. Lewis, O. Aharonson, J. P. Grotzinger, A. S. McEwen, and R. L. Kirk. 

1998 Global significance of cyclic sedimentary deposits on Mars. In Lunar and 

1999 Planetary Institute Science Conference Abstracts, volume 41 of Lunar and 

2000 Planetary Institute Science Conference Abstracts, page 2648, 2010. 

2001 G. E. Liston and J.-G. Winther. Antarctic Surface and Subsurface Snow 

2002 and Ice Melt Fluxes. Journal of Climate, 18:1469-1481, May 2005. doi: 

2003 10.1175/JCLI3344.1. 

2004 D. Loizeau, N. Mangold, F. Poulet, J.-P. Bibring, A. Gendrin, V. Ansan, 

2005 C. Gomez, B. Gondet, Y. Langevin, P. Masson, and G. Neukum. Phyllosili- 

2006 cates in the Mawrth Vallis region of Mars. Journal of Geophysical Research 

2007 (Planets), 112:E08S08, July 2007. doi: 10.1029/2006JE002877. 

2008 J.-B. Madeleine, F. Forget, J. W. Head, B. Levrard, F. Montmessin, and 

2009 E. Millour. Amazonian northern mid-latitude glaciation on Mars: A pro- 

2010 posed climate scenario. Icarus, 203:390-405, 2009. doi: 10.1016/j.icarus. 

2011 2009.04.037. 

2012 J.-B. Madeleine, F. Forget, E. Millour, L. Montabone, and M. J. Wolff. Revis- 

2013 iting the radiative impact of dust on Mars using the LMD Global Climate 

2014 Model. Journal of Geophysical Research (Planets), 116:E11010, November 

2015 2011. doi: 10.1029/2011JE003855. 

2016 M. C. Malin and K. S. Edgett. Sedimentary Rocks of Early Mars. Science, 

2017 290:1927-1937, 2000. doi: 10.1126/scicncc.290.5498.1927. 

2018 M. C. Malin and K. S. Edgctt. Mars Global Surveyor Mars Orbiter Cam- 

2019 era: Interplanetary cruise through primary mission. Journal of Geophysical 

2020 Research, 106:23429-23570, 2001. doi: 10.1029/2000JE001455. 

2021 M. C. Malin, K. S. Edgctt, B. A. Cantor, M. A. Caplinger, G. E. Danielson, 

2022 E. H. Jensen, M. A. Ravine, J. L. Sandoval, and K. D. Supulver. An overview 

2023 of the 1985-2006 Mars Orbiter Camera science investigation. International 

2024 Journal of Mars Science and Exploration, 5:1-60, 2010. doi: 10.1555/mars. 

2025 2010.0001. 

2026 M. Manga, A. Patel, J. Dufek, and E. S. Kite. Wet surface and dense 

2027 atmosphere on early Mars suggested by the bomb sag at Home Plate, 

2028 Mars. Geophysical Research Letters, 39:L01202, January 2012. doi: 

2029 10.1029/2011GL050192. 

2030 N. Mangold, L. Roach, R. MiUiken, S. Le Mouehc, V. Ansan, J. P. Bibring, 

2031 P. Masson, J. F. Mustard, S. Murchie, and G. Neukum. A Late Amazonian 

2032 alteration layer related to local volcanism on Mars. Icarus, 207:265-276, 

2033 2010. doi: 10. 1016/j.icarus.2009. 10.015. 

2034 N. Mangold, E.S. Kite, M.G. Kleinhans, H. Newsom, V. Ansan, E. Hauber, 

2035 E. Kraal, G. Quantin, and K. Tanaka. The origin and timing of fluvial 



71 



2036 activity at Eberswalde crater, Mars, accepted by Icarus. 2012. 

2037 N. Mangold, V. Ansan, P. Masson, and C. Vincendon. Estimate of aeolian dust 

2038 thickness in Arabia Terra, Mars: Implications of a thick mantle (>20 m) for 

2039 hydrogen detection. Geomorphologie - Relief Processus Environnement, (1): 

2040 23-31, 2009. ISSN 1266-5304. 

2041 C. V. Manning, C. P. McKay, and K. J. Zahnle. Thick and thin models of the 

2042 evolution of carbon dioxide on Mars. Icarus, 180:38-59, January 2006. doi: 

2043 10.1016/j.icarus.2005.08.014. 

2044 C. V. Manning, Y. Ma, D. A. Brain, C. P. McKay, and K. J. Zahnle. Paramet- 

2045 ric analysis of modeled ion escape from Mars. Icarus, 212:131-137, March 

2046 2011. doi: 10. 1016/j.icarus.2010. 11.028. 

2047 D. R. Marchant and J. W. Head. Antarctic dry valleys: Microclimate zonation, 

2048 variable geomorphic processes, and implications for assessing climate change 

2049 on Mars. Icarus, 192:187-222, December 2007. doi: 10.1016/j.icarus.2007. 

2050 06.018. 

2051 P. Marsh and M.-K. Woo. Wetting Front Advance and Freezing of Meltwa- 

2052 ter Within a Snow Cover 1. Observations in the Canadian Arctic. Water 

2053 Resources Research, 20:1853-1864, 1984. doi: 10.1029/WR020i012p01853. 

2054 M. Masse, O. Bourgeois, S. Le Mouelic, C. Verpoorter, A. Spiga, and L. Le 

2055 Deit. Wide distribution and glacial origin of polar gypsum on Mars. Earth 

2056 and Planetary Science Letters, 317:44-55, February 2012. doi: 10.1016/j. 

2057 epsl.2011. 11.035. 

2058 I. Matsuyama and M. Manga. Mars without the equilibrium rotational figure, 

2059 Tharsis, and the remnant rotational figure,. J. Geophys. Res., 115:E1020, 

2060 2010. doi: doi:10.1029/2010JE003686. 

2061 I. Matsuyama, J. X. Mitrovica, M. Manga, J. T. Perron, and M. A. Richards. 

2062 Rotational stability of dynamic planets with clastic lithospheres. Journal of 

2063 Geophysical Research (Planets), 111:E02003, February 2006. doi: 10.1029/ 

2064 2005JE002447. 

2065 C. P. McKay and W. L. Davis. Duration of liquid water habitats on early 

2066 Mars. Icarus, 90:214-221, April 1991. doi: 10.1016/0019-1035(91)90102-Y. 

2067 C. P. McKay, R. A. Wharton, Jr., S. W. Squyres, and G. D. Clow. Thickness 

2068 of ice on perennially frozen lakes. Nature, 313:561, February 1985. doi: 

2069 10.1038/313561a0. 

2070 B.C. McKelvey. The Lithologic Logs of DVDP Cores 10 and 11, Eastern Taylor 

2071 Valley. In McGinnis, L.D., editor. Dry Valley Drilling Project, volume 33 

2072 of Antarctic Research Series, page 63. American Geophysical Union, 1981. 

2073 N. K. McKcown, J. L. Bishop, E. Z. Noc Dobrea, B. L. Ehlmann, M. Parente, 

2074 J. F. Mustard, S. L. Murchie, G. A. Swayze, J. -P. Bibring, and E. A. Silver. 

2075 Characterization of phyllosilicates observed in the central Mawrth Vallis 

2076 region. Mars, their potential formational processes, and implications for past 

2077 climate. Journal of Geophysical Research (Planets), 114:E00D10, November 

2078 2009. doi: 10.1029/2008JE003301. 

2079 D. McKnight. Onyx River at Vanda gauge measurements. National Science 

2080 Foundation McMurdo Dry Valleys Long Term Ecological Research Signa- 



72 



2081 ture Datasets knb-lter-mcm.9022.4, 2011. 

2082 S. M. McLennan and J. P. Grotzinger. The sedimentary rock cycle of Mars. 

2083 In Bell, J., Ill, editor. The Martian Surface - Composition, Mineralogy, and 

2084 Physical Properties, volume 9 of Cambridge Planetary Science Series, page 

2085 541. Cambridge University Press, 2008. 

2086 S. M. McLennan, J. F. Bell, W. M. Calvin, P. R. Christensen, B. C. Clark, P. A. 

2087 de Souza, J. Farmer, W. H. Farrand, D. A. Fike, R. Gellert, A. Ghosh, T. D. 

2088 Glotch, J. P. Grotzinger, B. Hahn, K. E. Herkenhoff, J. A. Hurowitz, J. R. 

2089 Johnson, S. S. Johnson, B. JoUiff, G. Klingelhofer, A. H. Knoll, Z. Learner, 

2090 M. C. Mahn, H. Y. MeSween, J. Pocock, S. W. Ruff, L. A. Soderblom, S. W. 

2091 Squyres, N. J. Tosca, W. A. Watters, M. B. Wyatt, and A. Yen. Provenance 

2092 and diagenesis of the evaporite-bearing Burns formation, Meridiani Planum, 

2093 Mars. Earth and Planetary Science Letters, 240:95-121, November 2005. 

2094 doi: 10.1016/j.epsl.2005.09.041. 

2095 M. T. Mellon and B. M. Jakosky. The distribution and behavior of Martian 

2096 ground ice during past and present epochs. Journal of Geophysical Research, 

2097 1001:11781-11799, 1995. doi: 10.1029/95JE01027. 

2098 J. R. Michalski and E. Z. Noe Dobrea. Evidence for a sedimentary origin of 

2099 clay minerals in the Mawrth Vallis region. Mars. Geology, 35:951, October 

2100 2007. doi: 10.11302FG23854A.1. 

2101 R. E. Milliken and D. L. Bish. Sources and sinks of clay minerals on 

2102 Mars. Philosophical Magazine, 90:2293-2308, June 2010. doi: 10.1080/ 

2103 14786430903575132. 

2104 R. E. Milliken, J. P. Grotzinger, and B. J. Thomson. Paleoclimate of Mars as 

2105 captured by the stratigraphic record in Gale Crater. Geophysical Research 

2106 Letters, 370:L04201, 2010. doi: 10.1029/2009GL041870. 

2107 E. Millour, F. Forget, and S. R. Lewis. European Mars Climate Database 

2108 v4.3 Detailed Design Document. Technical report, Lab. de Meteorol. Dyn., 

2109 Paris. (Available at http://www-mars.lmd.jussieu.fr/), 2008. 

2110 D. W. Ming, R. V. Morris, and B. C. Clark. Aqueous alteration on Mars. In 

2111 Bell, J., Ill, editor. The Martian Surface - Composition, Mineralogy, and 

2112 Physical Properties, volume 9 of Cambridge Planetary Science Series, page 

2113 519. Cambridge University Press, 2008. 

2114 D. A. Minton and R. Malhotra. Assessing the Massive Young Sun Hypothesis 

2115 to Solve the Warm Young Earth Puzzle. Astrophysical Journal, 660:1700- 

2116 1706, 2007. doi: 10.1086/514331. 

2117 M. A. Mischna and M. I. Richardson. Climate simulation of recent climate 

2118 changes on Mars. In F. Forget, M. A. Lopez- Valverde, M. C. Desjean, 

2119 J. P. Huot, F. Lefevre, S. Lebonnois, S. R. Lewis, E. Millour, P. L. Read, 

2120 & R. J. Wilson, editor. Mars Atmosphere Modelling and Observations, page 

2121 312, 2006. 

2122 M. A. Mischna, M. I. Richardson, R. J. Wilson, and D. J. McCleese. On the 

2123 orbital forcing of Martian water and CO2 cycles: A general circulation model 

2124 study with simplified volatile schemes. Journal of Geophysical Research 

2125 (Planets), 108:5062, 2003. doi: 10.1029/2003JE002051. 



73 



2126 F. Montmessin, R. M. Haberle, F. Forget, Y. Langevin, R. T. Clancy, and 

2127 J. -P. Bibring. On the origin of perennial water ice at the south pole of 

2128 Mars: A precession-controlled mechanism? Journal of Geophysical Research 

2129 (Planets), 112:E08S17, August 2007. doi: 10.1029/2007JE002902. 

2130 G. A. Morgan, J. W. Head, F. Forget, J.-B. Madeleine, and A. Spiga. Gully 

2131 formation on Mars: Two recent phases of formation suggested by links bc- 

2132 tween morphology, slope orientation and insolation history. Icarus, 208: 

2133 658-666, 2010. doi: 10.1016/j.icarus.2010.02.019. 

2134 S. Murchie, L. Roach, F. Seelos, R. Milliken, J. Mustard, R. Arvidson, 

2135 S. Wiseman, K. Lichtenberg, J. Andrews- Hanna, J. Bishop, J.-P. Bibring, 

2136 M. Parcnte, and R. Morris. Evidence for the origin of layered deposits in 

2137 Candor Chasma, Mars, from mineral composition and hydrologic model- 

2138 ing. Journal of Geophysical Research (Planets), 114:E00D05, 2009a. doi: 

2139 10.1029/2009JE003343. 

2140 S. L. Murchie, J. F. Mustard, B. L. Ehlmann, R. E. Milliken, J. L. Bishop, 

2141 N. K. McKeown, E. Z. Noe Dobrea, F. P. Seelos, D. L. Buczkowski, S. M. 

2142 Wiseman, R. E. Arvidson, J. J. Wray, C. Swayze, R. N. Clark, D. J. Des 

2143 Marais, A. S. McEwen, and J.-P. Bibring. A synthesis of Martian aqueous 

2144 mineralogy after 1 Mars year of observations from the Mars Reconnaissance 

2145 Orbiter. Journal of Geophysical Research (Planets), 114:E00D06, 2009b. 

2146 doi: 10.1029/2009JE003342. 

2147 C. D. Murray and S. F. Dermott. Solar System Dynamics. February 2000. 

2148 J. F. Mustard, B. L. Ehlmann, S. L. Murchie, F. Poulet, N. Mangold, J. W. 

2149 Head, J.-P. Bibring, and L. H. Roach. Composition, Morphology, and 

2150 Stratigraphy of Noachian Crust around the Isidis basin. Journal of Geophys- 

2151 ical Research (Planets), 114:E00D12, 2009. doi: 10.1029/2009JE003349. 

2152 G. A. Neumann, J. B. Abshire, O. Aharonson, J. B. Garvin, X. Sun, and 

2153 M. T. Zubcr. Mars Orbiter Laser Altimeter pulse width measurements and 

2154 footprint-scale roughness. Geophysical Research Letters, 30 (11): 1561, June 

2155 2003. doi: 10.1029/2003GL017048. 

2156 C. E. Newman, S. R. Lewis, and P. L. Read. The atmospheric circulation and 

2157 dust activity in different orbital epochs on Mars. Icarus, 174:135-160, 2005. 

2158 doi: 10.1016/j.icarus.2004.10.023. 

2159 C.A. Nezat, W.B. Lyons, and K.A. Welch. Chemical weathering in streams 

2160 of a polar desert, Taylor Valley, Antarctica. Geological Society of American 

2161 Bulletin, 113(11):1401-1408, NOV 2001. 

2162 P. B. Nilcs and J. Michalski. Meridiani Planum sediments on Mars formed 

2163 through weathering in massive ice deposits. Nature Geoscience, 2:215-220, 

2164 2009. doi: 10.1038/ngeo438. 

2165 P. B. Niles and J. Michalski. Origin and Evolution of Sediments in Gale Crater 

2166 Through Ice-Hosted Processes. In Lunar and Planetary Institute Science 

2167 Conference Abstracts, volume 43 of Lunar and Planetary Institute Science 

2168 Conference Abstracts, page 2575, March 2012. 

2169 F. Nimmo and K. Tanaka. Early Crustal Evolution of Mars. Annual Review 

2170 of Earth and Planetary Sciences, 33:133-161, 2005. doi: 10.1146/annurev. 



74 



2171 carth.33.092203. 122637. 

2172 E. Z. Noc Dobrca, J. L. Bishop, N. K. McKcown, R. Fu, C. M. Rossi, J. R. 

2173 Michalski, C. Heinlein, V. Hanus, F. Poulet, R. J. F. Mustard, S. Murchie, 

2174 A. S. McEwen, G. Swayze, J.-P. Bibring, E. Malaret, and C. Hash. Miner- 

2175 alogy and stratigraphy of phyllosihcate-bearing and dark manthng units in 

2176 the greater Mawrth Valhs/west Arabia Terra area: Constraints on geological 

2177 origin. Journal of Geophysical Research (Planets), 115:E00D19, July 2010. 

2178 doi: 10.1029/2009JE003351. 

2179 C. H. Okubo. Strength and deformabihty of light-toned layered deposits ob- 

2180 served by MER Opportunity: Eagle to Erebus craters. Mars. Geophysical 

2181 Research Letters, 34:120205, October 2007. doi: 10.1029/2007GL031327. 

2182 A. A. Olsen and J.D. Rimstidt. Using a mineral lifetime diagram to evaluate 

2183 the persistence of olivine on Mars. American Mineralogist, 92(4):598-602, 

2184 2007. doi: {10.2138/am.2007.2462}. 

2185 M. M. Osterloo, F. S. Anderson, V. E. Hamilton, and B. M. Hynek. Geologic 

2186 context of proposed chloride-bearing materials on Mars. Journal of Geo- 

2187 physical Research (Planets), 115:E10012, 2010. doi: 10.1029/2010JE003613. 

2188 H. Pahke, R. D. Norris, J. O. Herrle, P. A. Wilson, H. K. Coxall, C. H. 

2189 Lear, N. J. Shackleton, A. K. Tripati, and B. S. Wade. The Heartbeat of 

2190 the Ohgocene Chmate System. Science, 314:1894-, December 2006. doi: 

2191 10.1126/science.ll33822. 

2192 J. T. Perron, J. X. Mitrovica, M. Manga, I. Matsuyama, and M. A. Richards. 

2193 Evidence for an ancient martian ocean in the topography of deformed shore- 

2194 lines. Nature, 447:840-843, 2007. doi: 10.1038/Nature05873. 

2195 R. J. Phillips, M. T. Zuber, S. C. Solomon, M. P. Golombek, B. M. Jakosky, 

2196 W. B. Banerdt, D. E. Smith, R. M. E. WiUiams, B. M. Hynek, O. Aharon- 

2197 son, and S. A. Hauck. Ancient Geodynamics and Global-Scale Hydrology on 

2198 Mars. Sctence, 291:2587-2591, March 2001. doi: 10.1126/science.l058701. 

2199 R. J. Phillips, B. J. Davis, K. L. Tanaka, S. Byrne, M. T. Mellon, N. E. Putzig, 

2200 R. M. Haberle, M. A. Kahre, B. A. Campbell, L. M. Carter, I. B. Smith, 

2201 J. W. Holt, S. E. Smrekar, D. C. Nunes, J. J. Plant, A. F. Egan, T. N. 

2202 Titus, and R. Seu. Massive CO2 Ice Deposits Sequestered in the South 

2203 Polar Layered Deposits of Mars. Science, 332:838-, 2011. doi: 10.1126/ 

2204 science. 1203091. 

2205 J. J. Plant, G. Picardi, A. Safaeinih, A. B. Ivanov, S. M. Milkovich, A. Ci- 

2206 cchetti, W. Kofman, J. Mouginot, W. M. Farrell, R. J. PhiUips, S. M. 

2207 Clifford, A. Frigeri, R. Orosei, C. Federico, I. P. Williams, D. A. Gurnett, 

2208 E. Nielsen, T. Hagfors, E. Heggy, E. R. Stofan, D. Plettemeier, T. R. Wat- 

2209 ters, C. J. Leuschen, and P. Edenhofer. Subsurface Radar Sounding of the 

2210 South Polar Layered Deposits of Mars. Science, 316:92-, April 2007. doi: 

2211 10. 1126/science. 1139672. 

2212 W. Pollard, C. Omelon, D. Andersen, and C. McKay. Perennial spring occur- 

2213 rence in the Expedition Fiord area of western Axel Heiberg Island, Canadian 

2214 High Arctic. Ganadian Journal of Earth Sciences, 36:105-120, 1999. doi: 

2215 10.1139/cjes-36-l-105. 



75 



2216 W.H. Press, S.A. Teukolsky, W.T. Vettcrling, and B.P. Flanncry. Numerical 

2217 Recipes 3rd Edition: The Art of Scientific Computing. Cambridge University 

2218 Press, 2007. 

2219 N. E. Putzig and M. T. Mellon. Apparent thermal inertia and the surface 

2220 heterogeneity of Mars. Icarus, 191:68-94, 2007. doi: 10.1016/j.icarus.2007. 

2221 05.013. 

2222 N. E. Putzig, M. T. Mellon, K. A. Krctkc, and R. E. Arvidson. Global thermal 

2223 inertia and surface properties of Mars from the MGS mapping mission. 

2224 Icarus, 173:325-341, February 2005. doi: 10.1016/j.icarus.2004.08.017. 

2225 M. I. Richardson and M. A. Mischna. Long-term evolution of transient liquid 

2226 water on Mars. Journal of Geophysical Research (Planets), 110:E03003, 

2227 2005. doi: 10.1029/2004JE002367. 

2228 M. I. Richardson and A. Soto. Controls on precipitation and aridity for an- 

2229 cient Mars. In Second Workshop on Mars Valley Networks, Smithson. Inst., 

2230 Moab, Utah., 2008a. 

2231 M. I. Richardson and A. Soto. Water at the surface of ancient Mars. In Second 

2232 Workshop on Mars Valley Networks, Smithson. Inst., Moab, Utah., 2008b. 

2233 L. H. Roach, J. F. Mustard, M. D. Lane, J. L. Bishop, and S. L. Murchie. 

2234 Diagenetic haematite and sulfate assemblages in Valles Marineris. Icarus, 

2235 207:659-674, June 2010. doi: 10. 1016/j.icarus.2009. 11.029. 

2236 M. R. Salvatore, J. F. Mustard, M. B. Wyatt, and S. L. Murchie. Definitive ev- 

2237 idence of Hesperian basalt in Acidalia and Chryse planitiae. Journal of Geo- 

2238 physical Research (Planets), 115:E07005, 2010. doi: 10. 1029/2009 JE003519. 

2239 N. Schorghofer. Dynamics of ice ages on Mars. Nature, 449:192-194, Septem- 

2240 ber 2007a. doi: 10.1038/nature06082. 

2241 N. Schorghofer. Theory of ground ice stability in sublimation environments. 

2242 Physical Review E, 75(4):041201, April 2007b. doi: 10.1103/PhysRevE.75. 

2243 041201. 

2244 N. Schorghofer. Fast numerical method for growth and retreat of subsurface 

2245 ice on Mars. Icarus, 208:598-607, 2010. doi: 10.1016/j.icarus.2010.03.022. 

2246 N. Schorghofer and O. Aharonson. Stability and exchange of subsurface ice on 

2247 Mars. Journal of Geophysical Research (Planets), 110:E05003, May 2005. 

2248 doi: 10.1029/2004JE002350. 

2249 D.H. Scott and K.L. Tanaka. 1-1802A: Geological map of the western equato- 

2250 rial regions of Mars. Technical report, U.S. Geological Survey, 1986. 

2251 M. M. Selvans, J. J. Plant, O. Aharonson, and A. Safaeinili. Internal structure 

2252 of Planum Boreum, from Mars advanced radar for subsurface and iono- 

2253 spheric sounding data. Journal of Geophysical Research (Planets), 115: 

2254 E09003, September 2010. doi: 10. 1029/2009 JE003537. 

2255 J. Shaw and T.R. Healy. Morphology of the Onyx River system, McMurdo 

2256 sound region, Antarctica. New Zealand Journal of Geology and Geophysics, 

2257 23(2):223-238, 1980. 

2258 D. E. Shean. Candidate ice-rich material within equatorial craters on 

2259 Mars. Geophysical Research Letters, 372:L24202, 2010. doi: 10.1029/ 

2260 2010GL045181. 



76 



2261 J. A. Skinner, Jr., T. M. Hare, and K. L. Tanaka. Digital Renovation of the At- 

2262 las of Mars 1:15, 000, 000-Scale Global Geologic Series Maps. In S. Mackwell 

2263 & E. Stansbery, editor, 37th Annual Lunar and Planetary Science Confer- 

2264 ence, volume 37 of Lunar and Planetary Inst. Technical Report, page 2331, 

2265 March 2006. 

2266 R. J. Soare, J. S. Kargel, G. R. Osinski, and F. Costard. Thcrmokarst processes 

2267 and the origin of crater-rim gullies in Utopia and western Elysium Planitia. 

2268 Icarus, 191:95-112, November 2007. doi: 10.1016/j.icarus.2007.04.018. 

2269 S. W. Squyres and Athena Science Team. Initial Opportunity Rover Results 

2270 at Endeavour Crater, Mars. In Lunar and Planetary Institute Science Con- 

2271 ference Abstracts, volume 43 of Lunar and Planetary Inst. Technical Report, 

2272 page 1892, March 2012. 

2273 S. W. Squyres and J. F. Kasting. Early Mars: How Warm and How Wet? 

2274 Science, 265:744-749, 1994. doi: 10.1126/science.265.5173.744. 

2275 K.L. Tanaka and D.H. Scott. I-1802C: Geological map of the polar regions of 

2276 Mars. Technical report, U.S. Geological Survey, 1987. 

2277 P. Thollot, N. Mangold, V. Ansan, S. Le Mouelic, R. E. Milliken, J. L. Bishop, 

2278 C. M. Weitz, L. H. Roach, J. F. Mustard, and S. L. Murchie. Most Mars 

2279 minerals in a nutshell: Various alteration phases formed in a single envi- 

2280 ronment in Noctis Labyrinthus. Journal of Geophysical Research (Planets), 

2281 117:E00J06, April 2012. doi: 10.1029/2011JE004028. 

2282 B.J. Thomson, N.T. Bridges, R. Milliken, A. Baldridge, S.J. Hook, J.K. Crow- 

2283 ley, G.M. Marion, C.R. de Souza Filho, A.J. Brown, and C.M. Weitz. Con- 

2284 straints on the origin and evolution of the layered mound in Gale Crater, 

2285 Mars using Mars Reconnaissance Orbiter data. Icarus, 214(2):413 - 432, 

2286 2011. doi: 10.1016/j.icarus.2011.05.002. 

2287 F. Tian, M. W. Claire, J. D. Haqq-Misra, M. Smith, D. C. Crisp, D. Catling, 

2288 K. Zahnle, and J. F. Kasting. Photochemical and climate consequences of 

2289 sulfur outgassing on early Mars. Earth and Planetary Science Letters, 295: 

2290 412-418, 2010. doi: 10.1016/j.epsl.2010.04.016. 

2291 T. N. Titus, H. H. Kieffer, and P. R. Christensen. Exposed Water Ice Dis- 

2292 covered near the South Pole of Mars. Science, 299:1048-1051, 2003. doi: 

2293 10. 1126/scicncc. 1080497. 

2294 O.B. Toon, J.B. Pollack, W. Ward, J. A. Burns, and K. Bilski. The astronom- 

2295 ical theory of climatic change on mars. Icarus, 44(3) :552 - 607, 1980. ISSN 

2296 0019-1035. doi: 10.1016/0019-1035(80)90130-X. 

2297 N. J. Tosca and A. H. Knoll. Juvenile chemical sediments and the long term 

2298 persistence of water at the surface of Mars. Earth and Planetary Science 

2299 Letters, 286:379-386, 2009. doi: 10.1016/j.epsl.2009.07.004. 

2300 J. Touma and J. Wisdom. The chaotic obhquity of Mars. Science, 259:1294- 

2301 1297, 1993. doi: 10.1126/science.259.5099.1294. 

2302 S. Turck-Chieze, L. Piau, and S. Couvidat. The Solar Energetic Balance Revis- 

2303 itcd by Young Solar Analogs, Helioseismology, and Neutrinos. Astrophysical 

2304 Journal Letters, 731:L29+, 2011. doi: 10.1088/2041-8205/731/2/L29. 

2305 D. T. Vaniman. Deriving a History of Deposition and Alteration from 



77 



2306 Phyllosilicate-Sulfatc Associations at Gale Crater. AGU Fall Meeting Ab- 

2307 stracts, pages abstract 7^P32A-04, December 2011. 

2308 M. Vincendon, F. Forget, and J. Mustard. Water ice at low to midlatitudes 

2309 on Mars. Journal of Geophysical Research (Planets), 115:E10001, 2010. doi: 

2310 10.1029/2010JE003584. 

2311 N. Warner. S. Gupta, J.-P. MuUer, J.-R. Kim, and S.-Y. Lin. A refined chronol- 

2312 ogy of catastrophic outflow events in Ares Vallis, Mars. Earth and Planetary 

2313 Science Letters, 288:58-69, October 2009. doi: 10.1016/j.epsl.2009.09.008. 

2314 N. H. Warner, S. Gupta, J.-R. Kim, J.-P. MuUer, L. Le Corre, J. Morley, S.- 

2315 Y. Lin, and C. McGonigle. Constraints on the origin and evolution of lani 

2316 Chaos, Mars. Journal of Geophysical Research (Planets), 116:E06003, June 

2317 2011. doi: 10.1029/2010JE003787. 

2318 J.K. Warren. Evaporites: sediments, resources and hydrocarbons. Springer, 

2319 2006. 

2320 S. G. Warren. Optical Properties of Snow. Reviews of Geophysics and Space 

2321 Physics, 20:67, February 1982. 

2322 S. G. Warren. Review. Impurities in snow: effects on albedo and snowmelt. 

2323 Annals of Glaciology, 5:177-179, 1984. 

2324 S. G. Warren and R. E. Brandt. Optical constants of ice from the ultraviolet 

2325 to the microwave: A revised compilation. Journal of Geophysical Research 

2326 (Atmospheres), 113:D14220, July 2008. doi: 10.1029/2007JD009744. 

2327 S. G. Warren and W. J. Wiscombe. A Model for the Spectral Albedo of 

2328 Snow. II: Snow Containing Atmospheric Aerosols. Journal of Atmospheric 

2329 Sciences, 37:2734-2745, 1980. doi: 10.1175/1520-0469(1980)037$($2734: 

2330 AMFTSA$)$2.0.CO;2. 

2331 C. M. Weitz, M. D. Lane, M. Staid, and E. N. Dobrea. Gray hematite 

2332 distribution and formation in Ophir and Candor chasmata. Journal of 

2333 Geophysical Research (Planets), 113:E02016, February 2008. doi: 10.1029/ 

2334 2007JE002930. 

2335 L. Wendt, C. Gross, T. Kneissl, M. Sowe, J.-P. Combe, L. Ledeit, P. C. 

2336 McGuire, and G. Neukum. Sulfates and iron oxides in Ophir Chasma, Mars, 

2337 based on OMEGA and CRISM observations. Icarus, 213:86-103, May 2011. 

2338 doi: 10.1016/j.icarus.2011.02.013. 

2339 J. Whiteway, L. Komguem, and C. Dickinson. Observations of Mars atmo- 

2340 spheric dust and clouds with the lidar instrument on the Phoenix mission. 

2341 In F. Forget & E. Millour, editor. Mars Atmosphere: Modelling and obser- 

2342 vation, pages 171-174, February 2011. 

2343 K. E. Williams, O. B. Toon, J. L. Heldmann, C. McKay, and M. T. Mellon. 

2344 Stability of mid-latitude snowpacks on Mars. Icarus, 196:565-577, 2008. 

2345 doi: 10.1016/j.icarus.2008.03.017. 

2346 K. E. WiUiams, O. B. Toon, J. L. Heldmann, and M. T. Mellon. Ancient 

2347 melting of mid-latitude snowpacks on Mars as a water source for gullies. 

2348 Icarus, 200:418-425, 2009. doi: 10. 1016/j.icarus.2008. 12.013. 

2349 R. M. E. Williams. Global Spatial Distribution of Raised Curvilinear Features 

2350 on Mars. In Lunar and Planetary Institute Science Gonference Abstracts, 



78 



2351 volume 38 of Lunar and Planetary Institute Science Conference Abstracts, 

2352 page 1821, 2007. 

2353 S. A. Wilson, A. D. Howard, J. M. Moore, and J. A. Grant. Geomorphic 

2354 and stratigraphic analysis of Crater Terby and layered deposits north of 

2355 Hellas basin. Mars. Journal of Geophysical Research (Planets), 112:E08009, 

2356 August 2007. doi: 10.1029/2006JE002830. 

2357 S. A. Wilson, J. M. Moore, A. D. Howard, and D. E. Wilhelms. Evidence for 

2358 ancient lakes in the Hellas region, pages 195-222. 2010. 

2359 S. A. Wilson, J. A. Grant, and A. D. Howard. Distribution of Intracrater 

2360 Alluvial Fans and Deltaic Deposits in the Southern Highlands of Mars. In 

2361 Lunar and Planetary Institute Science Conference Abstracts, volume 43 of 

2362 Lunar and Planetary Inst. Technical Report, page 2462, March 2012. 

2363 S. M. Wiseman, J. C. Andrews-Hanna, R. E. Arvidson, J. F. Mustard, and 

2364 K. J. Zabrusky. Distribution of Hydrated Sulfates Across Arabia Terra 

2365 Using CRISM Data: Implications for Martian Hydrology. In Lunar and 

2366 Planetary Institute Science Conference Abstracts, volume 42 of Lunar and 

2367 Planetary Institute Science Conference Abstracts, page 2133, March 2011. 

2368 M. J. Wolff and R. T. Clancy. Constraints on the size of Martian aerosols 

2369 from Thermal Emission Spectrometer observations. Journal of Geophysical 

2370 Research (Planets), 108:5097, 2003. doi: 10.1029/2003JE002057. 

2371 M. J. Wolff, M. D. Smith, R. T. Clancy, N. Spanovich, B. A. Whitney, M. T. 

2372 Lemmon, J. L. Bandfield, D. Banfield, A. Ghosh, G. Landis, P. R. Chris- 

2373 tcnsen, J. F. Bell, and S. W. Squyres. Constraints on dust aerosols from 

2374 the Mars Exploration Rovers using MGS overflights and Mini-TES. Jour- 

2375 nal of Geophysical Research (Planets), 111:E12S17, December 2006. doi: 

2376 10. 1029/2006 JE002786. 

2377 M. J. Wolff, M. D. Smith, R. T. Clancy, R. Arvidson, M. Kahre, F. See- 

2378 los, S. Murchic, and H. Savijarvi. Wavelength dependence of dust aerosol 

2379 single scattering albedo as observed by the Compact Reconnaissance Iniag- 

2380 ing Spectrometer. Journal of Geophysical Research (Planets), 114:E00D04, 

2381 June 2009. doi: 10.1029/2009JE003350. 

2382 B. E. Wood, H.-R. Miiller, G. R Zank, J. L. Linsky, and S. Redfield. New 

2383 Mass-Loss Measurements from Astrospheric Lja Absorption. Astrophysical 

2384 Journal Letters, 628:L143-L146, August 2005. doi: 10.1086/432716. 

2385 R. Wordsworth, F. Forget, and V. Eymet. Infrared coUision-induced and far- 

2386 line absorption in dense CO2 atmospheres. Icarus, 210:992-997, 2010. doi: 

2387 10.1016/j.icarus.2010.06.010. 

2388 R. Wordsworth, F. Forget, F. Sclsis, and R. Pierrehumbert. How 

2389 important is CO2 to planetary habitability? In Exoclimes 

2390 2: The Diversity of Planetary Atmospheres, page downloaded 

2391 from http://www.astro.ex.ac.uk/exoclimes/2012/pdf/talks/ 

2392 Day02_Wordsworth.pdf, 2012. 

2393 J. J. Wray, B. L. Ehlmann, S. W. Squyres, J. F. Mustard, and R. L. Kirk. 

2394 Compositional stratigraphy of clay-bearing layered deposits at Mawrth 

2395 VaUis, Mars. Geophysical Research Letters, 35:L12202, June 2008. doi: 



79 



2396 10.1029/2008GL034385. 

2397 J. J. Wray, R. E. Milliken, C. M. Dundas, G. A. Swayze, J. C. Andrews- 

2398 Hanna, A. M. Baldridge, M. Chojnacki, J. L. Bishop, B. L. Ehlmann, S. L. 

2399 Murchie, R. N. Clark, F. P. Seelos, L. L. Tornabene, and S. W. Squyres. 

2400 Columbus crater and other possible groundwater-fed paleolakes of Terra 

2401 Sircnum, Mars. Journal of Geophysical Research (Planets), 116:E01001, 

2402 2011. doi: 10.1029/2010JE003694. 

2403 R Yang, B.-C. Gao, W. J. Wiscombe, M. I. Mishchenko, S. E. Platnick, H.-L. 

2404 Huang, B. A. Baum, Y. X. Hu, D. M. Winker, and S.-G. Tsay. Inherent and 

2405 apparent scattering properties of coated or uncoated spheres embedded in 

2406 an absorbing host medium. Applied Optics, 41:2740-2759, May 2002. doi: 

2407 10. 1364/A0.41. 002740. 

2408 A. S. Yen, R. Gellert, C. Schr5der, R. V. Morris, J. F. Bell, A. T. Knud- 

2409 son, B. C. Clark, D. W. Ming, J. A. Crisp, R. E. Arvidson, D. Blaney, 

2410 J. Briickner, P. R. Christensen, D. J. Desmarais, P. A. de Souza, T. E. 

2411 Economou, A. Ghosh, B. C. Hahn, K. E. Herkenhoff, L. A. Haskin, J. A. 

2412 Hurowitz, B. L. Joliff, J. R. Johnson, G. Klingclhofcr, M. B. Madsen, 

2413 S. M. McLennan, H. Y. McSween, L. Richter, R. Rieder, D. Rodionov, 

2414 L. Soderblom, S. W. Squyres, N. J. Tosca, A. Wang, M. Wyatt, and J. Zipfel. 

2415 An integrated view of the chemistry and mineralogy of martian soils. Na- 

2416 ture, 436:49-54, July 2005. doi: 10.1038/nature03637. 

2417 K. Zabrusky, J.C. Andrews-Hanna, and S.M. Wiseman. Reconstructing the 

2418 Distribution and Depositional History of the Sedimentary Deposits of Ara- 

2419 bia Terra, Mars. Icarus, (0):-, 2012. 

2420 A. M. Zalucha, R. A. Plumb, and R. J. Wilson. An Analysis of the Effect of 

2421 Topography on the Martian Hadlcy Cells. Journal of Atmospheric Sciences, 

2422 67:673-693, March 2010. doi: 10.1175/2009JAS3130.1. 

2423 M. Zanetti, H. Hicsingcr, D. Reiss, E. Haubcr, and G. Ncukum. Distribution 

2424 and evolution of scalloped terrain in the southern hemisphere. Mars. Icarus, 

2425 206:691-706, April 2010. doi: 10.1016/j.icarus.2009.09.010. 

2426 J.R. Zimbelman and S.P. Scheldt. Hesperian Age for Western Medusae Fossae 

2427 Formation, Mars. Science, 2012. doi: 10. 1126/scicnce. 1221094. 

2428 M. T. Zuber, R. J. Phillips, J. C. Andrews-Hanna, S. W. Asmar, A. S. Kono- 

2429 pliv, F. G. Lemoine, J. J. Plant, D. E. Smith, and S. E. Smrekar. Density of 

2430 Mars' South Polar Layered Deposits. Science, 317:1718-, September 2007. 

2431 doi: 10.1126/science.ll46995. 



80 



